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LIST OF SYMBOLS 


A Area of a structural member cross section (sq. in.) 

Tributary area of a structure for p\irpose of confuting load 
(sq. ft) 







Area of steel beam cross section, coii5)Osite beam construction only 
(sq. in.) 

Net area of the rear face of a rectangular structure with openings 
(sq. ft) 

Area of concrete in cross section (sq. in.) 

Area of flange of structural steel beam or girder (sq. in.) 

Net area of the exterior front wa3J. of a rectangular structvire 
with openings (sq. ft) 

Cross section area of concrete column (sq. in.) 

Gross area of tne exterior front face of a rectangular structure 
(sq. ft) 

Area of each portion of the subdivided face of a rectangular struc¬ 
ture with openings (sq. ft) 

Area of openings in the front wall of a rectargular structure 
(sq. ft) 

Net area of the back wall of a rectangular structure with openings 
(sq. ft) 


A Area of reinforcing steel bars peuallel to steel beams in composite 

^ construction (sq. in.) 

A Area of tension steel in reinforced concrete memiber (sq. in.) 

s 

Ay Area of one stirrup of the web reinforcment of a reinforced .con¬ 

crete menfoer (sq. in. ) 


A Area of web of steel beam (sq. in.) 

w 

A* Area of conq)ression steel in reinforced concrete member (sq. in.) 

s 

a Clear distance between flanges of a steel beam (in.) 

Depth of compression area of concrete (in.) 

Short side dimension of two-way slabs (ft) 

Column spacing in flat slab construction (ft) 

Acceleration of mass (ft/sec^) 


*g,l^ \,n-l» ®g,n*l 


Acceleration of the g^ floor mass at times t^, 

t^, t ,, t T, respectively 
2 * n-1' ml^ ^ 


»g(t) 

* 1 ' * 2 ** 


Acceleration of the 



floor mass as a function of time 


.a 


h-l» n" “nfl 


Acceleration at time t^, t^, "^n+l^ 

respectively 
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Distance from edge "a" to centroid of the load on area "A” 

Distance from edge "a" to the centroid of t)he inertial force 
Peak value of externally applied loadL 

Elastic stiffness coefficient for soil under rocking foundation 

Equivalent peak value of the externally (applied load 

till 

Peak value of the effective load on the/g floor mass of a 
multi-story building 

"til 

Peak value of the external load, f (t), on the g floor mass 

S ! 

Width of flange of steel beam (in.) J 

Width of compression flange of reinforced concrete beam (in.) 

Long side dimension of two-way slahs (ft) 

Flange width of channel section () 

Width of longitudinal web stiffer/er for. steel beams (in.) 

j 

Width of stem of reinforced concrete tee beam (in.) 

Width of flange of steel beam u$ed in composite construction (in.) 
Distance from edge "b" to centrpid of load on area "B" 

Distance from edge "b" to the qentroid of the inertial force 

i 

Compression force developed in i reinforced concrete meniber (kips) 
Static compression colinnn influence factor for shear vails 

Crater depth factor for underj^ound burst 

Dynamic compression column i^luence factor for shear walls 
Average drag coefficient / 

Local drag coefficient 

Maximum ACI Code moment coefficient obtainable for any point on 
a slab (equation 6 * 52 ) 

Ratio of maximum resistance to peak load, C = R /b 

R 

Parameter defined by equation 5 . 9 U 
Seismic wave velocity (fps) 

Ratio of load duration to natural period of oscillation, =* t/T 

Approximate coefficient for determining elastic truss deflection 
(varies with type of truss) 

W 

Ratio of maximum work done to absolute maximum work done, C « — 

W W 

Local drag coefficient for a cylindrical segment of infinite len^h 

Value of determinate formed by omitting the row and column con¬ 
taining the coefficient from the matrix [cj 
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c Distance from neutral axis to the extreme fiber (in.) 

Cohesive strength of the soil as determined by conventional 
laboratory testing procedures (psi) 

The clear distance between compression flange and stiffener 
c^ Velocity of sound in undisturbed air (fps) 

Crefi Velocity of sound in the region of reflected overpressure (fps) 

D Diameter of a cylindrical segment or a spherical dome (ft) 

Diameter of the smallest longitudinal reinforcing bar (in.) 

Diameter of spiral core of a concrete column 

D.L.F Dynamic Load Factor = x /x 

m s 

D-I.F Dynamic Increase Factor for strength of materials 

d Depth of burst of bomb below ground surface (ft) 

Height of burst of bomb above ground surface (ft) 

Depth from the compression face of a beam or slab to the centroid 
of longitudinal tensile reinforcement (in.) 

Least lateral dimension of rectangular column (in.) 

Total depth of steel beam (in.) 

Over-all depth of truss (ft) 

Diameter or width of column capital in flat slab constiuction (ft) 
Diameter of circular reinforced concrete column (in.) 

d* Distance between centroids of compression and tensile steel in 

doubly reinforced concrete member 

E Modulus of elasticity, including Young*s Modulus for steel and 

soil, of any structural material except concrete 
Energy absorbed by the equivalent system 
Explosive factor for cratering 

Conpressive modulus of elasticity of concrete, kips/sq. in. 

Initial tangent modulus of elasticity for concrete 
Explosive factor for underground burst 
Maximum allowable energy absorption for the g^^ story 
Base of natural logarithms (2.718) 

Eccentricity of load, distance from gravity axis to point of 
application of load (in.) 

Center to center spacing of two trusses 

Strain at failure (‘^) 

Strain at initiation of strain hardening (^) 

Strain at yield {^o) 

Horizontal force 

"th 

Sum of all external loads above the g story, f(t) + f -(t)+.... 
4 f (t) ® 

Design lateral load on frame at time t 
n ^ n 



VII 




EM IIIO-345-IH5 
15 Mar 57 


Symbols 


Summation of all external horizontal forces applied to structure*, 
including foundation reactions 

F Pressure factor for underground burst 

Normal component of total passive resistance force (kips) 

F(t) Load applied to a structural element or system as function of time 

f Fiber stress (psi) 

f Column critical buckling stress (psi) 

fg^^ Average stress at design level due to axial load (psi) 

fg Dynamic yield strength for structural steel in composite beam 

construction only (psi) 

f^ Elastic buckling stress due to bending (psi) 

f^^ Maximum stress due to bending at design level for steel members 
subjected to axial load and bending (psi) 

f, Dynamic yield strength of steel (psi) 

"till 

f Effective load on the g floor of a multi-story structure as a 

function of time 

f^ Initial static stress (psi) 

Lower static stress (psi) 

f^ Dynamic yield strength for reinforcing steel for con 5 )Osite beam 

construction only (psi) 

f^ Tensile strength of concrete (psi) 

f(t) Load applied to a structural element of system as a fxmetion of 

time 

f Upper static yield point stress (psi) 

f static yield point stress 

y 

f^, fp, f_ Factors used in the design of tee beams 

^ j "bh 

f^(t), f2(t), f (t) Load applied to the first, second, and g floors 

° of a multi-story building as a function of time 

f^ Static ultimate compressive strength of concrete 

Dynamic ultimate compressive strength of concrete 

G Solidity ratio of truss 

Shear modulus of elasticity (psi) 

G.Z. Ground zero (point on ground below point of detonation of air 
burst) 

g Acceleration of gravity 

H Impulse per unit area, in^iulse per foot of length, or total in^Dulso 

(kip-sec) 

Depth of footing or wall below the surface of the ground 
Height of shear wall to center of roof beam 
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Equivalent ln5)ulse acting on the equivalent system (kip-sec) 

"til 

Total impulse of the external load on the g floor mass (kip-sec) 

Total inpulse of the external load (kip-sec) 

Height (ft) 

Story height (ft) 

Column height (ft) 

Height of a structure above ground (ft) 

Maximum thickness of channel flange (in.) 

Clear height of column (ft) 

"th 

Height of the g story of a multi-story building (ft) 


Clearing dimension for each portion of the subdivided front face 
of rectangular structure with openings (ft) 

Clearing height of a closed rectangular structure (ft) 

Vertical distance from roof to point under consideration for bturied 
structure (ft) 

Weighted average buildup height of the exterior rear face of a 
rectangular structure with openings (ft) 

Weighted average clearing height for the exterior front face of 
a rectangular structure with openings (ft) 

Weighted average clearing height for the interior rear face of a 
rectangular structure with the openings (ft) 

Weighted average buildup height for the interior front face of a 

rectangular struct!ire with openings (ft) 

( h k\ 

Moment of inertia I in. or ft ) 

Inertial force of "A" portion of two-way slab 


Average of the gross and transformed moments of inertia 
(in.^ or ft^) 

Inertial force of "B" portion of two-way slab 

( k 

Moment of Inertia of the gross section or ft j 

( ' 

Moment of inertia of the transformed section I in. or ft 


Moment of inertia of rectangular section in concrete tee beam 

I 

design ^ in. j 

Moment of. inertia of tee section in concrete tee beam design (in. 
l(t^j Velocity inqiulse at time, t^ 

Nteiss moment of inertia of structure about axis of rotation "O" 


Ratio of distance between centroid of compression and centroid 
of tension to the depth, d 

Modulus of strain hardening 


DC 
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Symbols 


KE 

(KE) 

(KE) 




So 


S0 


S 

K 


KT 


K’ 


K" 


k 


k 

c 

k 

e 

k 

ep 


k 

g 


Kinetic energy 

Kinetic energy of the actual system 

i 

Kinetic energy of the equivalent system 

Either load, mass, or load-mass factor for slab fixed on fotir sides 

(i = 1, 2, 3,..-; j 2, 3—•) Stiffness influence coefficient 

(kips/ft) 

Load factor 
Load-mass factor 
Mass factor 

Normal component of passive pressiire coefficient accounting for 
the cohesive effect of any soil having an internal friction equal 
to 0 

Normal component of passive pressure coefficient for any soil with 
internal friction angle equal to 0 and with zero wall friction 
developed 

Normal conponent of passive pressure coefficient for any soil with 
internal friction angle equal to 0 and with maximum wall friction 
developed 

Resistance factor 

Either load, mass, or load-mass factor for slab sinply supported 
on four sides 

Kiloton, 1000 tons 

Either load, mass, or load-mass factor for special edge conditions 

Ratio of the maximum average overpressure to the reflected over¬ 
pressure existing on an inclined roof 

Beam equivalent length coefficient used in the design of steel 
beam columns 

Column length factor used in the design of steel columns and beam 
columns 

Spring constant, force required to cause unit deflection of spring 
(kips/ft) 

Soil pressure factor (psi) 

Soil constant for underground explosion (psi) 

Effective spring constant (kips/ft) 

Equivalent spring constant (kips/ft) 

Spring constant in the elasto-plastic range 

th 

k^j-'-k^-"- C'pring conatants for the first, second, and the g 
stories of a multi-story building (kips/ft) 

Simplified notation of k 

g,g-l 

X 


^ynibols 





gi 


^'4 

gi 


gl' 


L' 


M 


M 


M 


cn 


M 


cp 


“d 


M 


dy 


M 


ed 


M 


m 


M 


fa 


M 


tb 


M 


tnn 


M 


nqp 
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^rlng constant of the coupling spring between the g^^ and the 
floors of a multi-story building (kips/ft) 


k* = 
gi 


gi 




Spring constants of the springs 
the masses 


connecting m 


g 


with 


Spein length of beam or truss (ft) 

Unsupixjrted length of beam (ft) 

Length of a rectangular structure in the direction of propagation 
of the blast wave (ft) 

Length of a cylindrical segment along the cylinder axis (ft) 

Length of shear wall center to center of column steel 
Spacing of columns in each direction 

Distance from the outside of the front face to the inside rear 
face of a rectangular structure with openings (ft) 

Distance from the front face of a rectangular structure to a point 
under consideration on the roof or sides in the direction of pro¬ 
pagation of the blast wave (ft) 

Length between sections of zero and maximum moment being considered 
(composite beams) 

Bending moment applied to a section 

Moment of forces on piles about their centroidal axis 
Resisting moment of soil on the footing per unit width 

Bending moment at center line of a beam or slab 

Negative resisting moment in colxnnn per foot, flat slabs (kip-ft) 

Positive resisting moment in column strip per foot, flat slabs 
(kip-ft) 

Maximum design moment in a member \mder axial load, Pj^ 

Dynamic bending moment 
Elastic dynamic buckling moment 


Maximum bending moment 

Component in a plane perpendicular to edge "a" of the total re¬ 
sisting moment along the fracture lines bounding area "A," two- 
way slabs (kip-ft) 

Con5)onent in a plane perpendicular to edge "b," of the total re¬ 
sisting moment along the fracture lines bounding area B, two-way 
slabs (kip-ft) 

Negative resisting moment in middle strip per fpot, flat slabs 
(kip-ft) 

Positive resisting moment in middle strip per foot, flat slabs 
(kip-ft) 


I 


XI 
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^■^Pcn 

^cp 

^fa 


«Pfb 


*S>nip 

M 

pos 

»^sa 

^sb 

^b 


Negative bending moment at support 

Moment of all external forces on structure about axis of rotation 
" 0 " 

oummation of external moments about the point of rotation ex¬ 
cluding footing projection moments 

Plastic resisting moment under bending orG.y 

Negative plastic resisting moment in coliamn strip per foot, flat 
slabs (kip-ft) 

Positive plastic resisting moment in colxunn strip per foot, flat 
slabs (kip-ft) 

Component of the total plastic bending moment capacity along the 
fracture line boundary of area "A" which is in a plane perpendic¬ 
ular to edge "a," or the total positive plastic bending moment 
capacity for a section parallel to edge "a," two-way slabs (kip-ft) 

Con5>onent of the total plastic bending moment capacity along the 
fracture line boundary of area "B" which is in a plane perpendic¬ 
ular to edge "b," or the total positive plastic bending moment 
capacity for a section parallel to edge "b,” two-way slabs (kip/ft) 

Plastic resisting moment at centerline of beam or slab 

Negative plastic resisting moment in middle strip per foot, flat 
slabs (kip-ft) 

Positive plastic resisting moment in middle strip per foot, flat 
slabs (kip-ft) 

Maximum positive bending moment (kip-ft) 

Plastic resisting moment at support (kip-ft) 

Total negative plastic bending moment capacity along edge "a,” 
two-way slabs (kip-ft) 

Total negative plastic bending moment capacity along edge "b," 
two-way slabs (kip-ft) 

Yield resisting moment (kip-ft) 

Plastic resisting moment per unit of width of slab (kip-ft/ft) 

Plastic positive bending moment capacity per unit width for short 
span, two-way slabs (kip-ft/ft) 




Plastic negative bending moment capacity per unit width at center 
of edge "a" for long span, two-way slabs (kip-ft/ft) 

Plastic negative bending moment capacity per unit width at center 
of edge ''b'' for short span, two-way slabs (kip-ft/ft) 
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Moment at the intersection of the two-linear portions of the P vs ! 
curve of a steel beam section (kip-ft) 

Theoretical plastic resisting moment (kip-ft) 

Mass per unit length (kip-sec^/ft^) 

Point mass (kip-sec' ^/ft) 

In concrete design: f^^O.85 

Total moving mass of structure and earth included betveen footings 
Number of fundamental dimensional q[uantities 


Mass of the equivalent system I kip-sec 


Vft") 


1 ^ 2 


back 


Mass of structure considered to rotate as well as translate 
Total mass of the element or structural system under consideration 
(kip-sec ^/ft) 

Total mass of the equivalent system (kip-sec' -/ft) 

th 

.m Mass of the firsts second, and g floors of a multi-story 
® building 

Number of non-dimensional parameters 

Total number of shear connectors required between the points of 
zero and maximum moment of a composite beam 
Weighted number of piles in group 

Ratio of the length of fixed-edge perimeter to the total perimeter. 

Ratio of the length of simply-supported perimeter to the total 
perimeter 

Number of columns in a story 

Ratio of modulus of elasticity of steel to modulus of elasticity 
of concrete 

Number of dimensional variables 
Number of stories in a multi-story frame 

Load or force (kips) 

Total dynamic load on slab (kips) 

Panel influence factor for shear walls 

Total load on portion of two-way slab (kips) 

Average axial load acting on each of several columns of a frame 
(kips) 

Local overpressure on the back face of a buried rectangular 
structure (psi) 

Total load on ’’B" portion of two-way slab (kips) 

Compression mode overpressure (psi) 

Uniform compressive pressure applied radially on arch 
Uniform radial pressure that produces buckling of arch. 
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Symbols 


'cyl 




dome 


^ed 

p^(t) 


■front 


g 


3 

-P 

^refl 

p 

'r-a 


Local overpressxare normal to the exterior surface of a cylindrical 
segment (psi) 

Maximum axial load on colimm with given (kips) 

Reduced maximum axial load for long colvunns 

Deflection mode overpressiure (psi) 

Total dynamic axial load 

Local overpressure normal to the exterior face of a spherical 
dome (psi) 

■s^iuivalent concentrated load for equivalent system (kips) 

Elastic dynamic buckling load (kips) 

Equivalent load on an element as a function of time 

Average value of P(t) for the far (or leeward) side of the arch 
(psi) 

Local overpressure on the front face of a buried rectangular 
structure (psi) 

Peak underground overpressure resulting at a given location from 
an underground burst (psi) 

Total vertical load (blast plus static) on column 
Force acting at any time, t^ 

Average value of P(t) for the near (or windward) side of the arch 
(psi) 

Dynamic plastic axial load capacity of coliunn 

Reflected shockwave overpressixre for angle of incidence of zero 
degrees (psi) 

Reflected shock wave overpressui'e for angle of incidence other 
than zero degrees (psi) 


P Overpressure existing in the incident shock wave for any value of 

" t-t^ (psi) 

^sb Overpressure existing in the incident shock wave when t-t^ = t^ 

^ (psi) 


P Local overpressure on the exterior side walls of a rectangular 

structure 


Initial peak incident overpressure (psi) 

Peak overpressure of the shock wave formed in the interior of a 
structure with openings (psi) 

P . Stagnation overpressure; the overpressure existing in a region in 
Vhich the moving air has been brought con^l^tely to rest (psi) 

P^ Axial load on the columns of a frame due to the vertical live and 
dead loads (kips) 
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Syndbols 

[r.P J T^ie time average of the total vertical load on the column in the 

' time interval t < t < t^ 

ge gm 

P- . Average overpressure on the exterior hack wall of a rectangular 
^ structure (psi) 

P, . Average net oveirpressure acting on back wall of a rectangular 
structure (psi) 

l^backlmax value of the average overpressure on the exterior back 

' c /max „3^21 of a rectangular structure (psi) 

P , Average overpressxire on the closed ends of a cylindrical segment 
(psi) 

P_ . Average net overpressure acting on front wall of a rectangular 
structure (psi) 

P Average overpressure on the exterior of the front wall of a 

rectangular structure (psi) 

^ + Average overpressure on the interior front wall of a rectangular 

structiire with openings (psi) 

P. -- Reflected shock wave average overpressure in the interior of a 
rectangular structure with openings (psi) 

F. - Average overpressure on the interior roof of a rectangular struc- 
ture with openings (psi) 

F. .. Average oveipress\ire on the interior side walls of a rectangular 
structure with openings (psi) 

P . Net average horizontal overpressure exerted on a rectangular 
“ structure (psi) 

P Average net overpressure acting on the roof of a rectangular 

structure (psi) 

P _ Average overpressure on the front slope of a gable roof (psi) 

Average overpress\ire on the exterior of a rectangular structvire 
(psi) 

F' - Peak average overpressure on the front slope of the first gable 
of a multi-gabled roof (psi) 

? . Average net overpressure acting inward on sidewall of rectangular 

structure (psi) 

F .. Average overpressvire on the exterior side walls of a rectangular 
structure (psi) 

P^ Axial load determined by the intersection of the two linear 

portions of the P vs M curve of a steel beam section (kips) 

Actual load on a structural element as a fiinction of time 

Ratio of tensile reinforcement in reinforced concrete members to 
concrete area, 

Uniformly-distribute'd load intensity 

XV 
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Symbols 


Critical steel ratio for reinforced concrete member 
Equivalent static load for an arch (psi) 

Ratio of volume of spiral reinforcement to the volume, of the 
concrete core (out-to-out of spirals) of a spirally reinforced 
concrete colimm 

Ratio of compressive reinforcement in beams to the concrete area, 
A^/bd 

Statical moment of the section about the centroidal axis 
Strength of a shear connector in coH^xjsite construction 

Radiant energy on a unit area (cal/ cm^) 

Unit drag pressure produced by the incident shock wave (psi) 
Dynamic design bearing pressure on soil ^ kips/ft] 

Maocimum bearing pressure on soil (kips/ft^'j 

Maximum unit drag pressure produced by the Incident shock wave 
(psi) 

Radius of a spherical dome 

Total resistance of structural element or structural system (kips) 
Dosage of gamma radiation (roentgens) 

Crater radius (ft) 

Total resistance of "A" portion of two-way slabs (kips) 

Crater radius (ft) 

Horizontal static load resistance of shear wall at Tirst cracking 

Horizontal dynamic load resistance of shear wall at first cracking 

Reynolds number of the high velocity wind, in the incident shock 
wave 

Resistance of the equivalent system (kips) 

Simplified notation of R ,, the resistance function developed 

gj g“-i- 

between the and (g-l)^^ floors of a multi-story building 
The average value of R in the time interval from 0 to t]^ 

Resistance function developed between the and the i^^ floors 
of a multi-story building (kips) 

th 

Msiximum resistance of the g story of a multi-story building 
(kips) 

Resistance developed in the spring which connects m to the ground 

"til ^ 

Total resistance acting on the g floor mass (kips) 

Maximum resistance developed by a structural system (kips) 

Plastic resistance of the "A" portion of two-way slabs where the 
edge is fixed (kips) 
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nB 


R 


me 


R 


mf 


R 


u 


R 

3 

E, 


R 


Im 


1 mA. 


^ InrR 

RCt) 

H(t) 

cv*o 

*0 

R(x) 


R 


g 


Plastic resistance of the "B" portion of two-way slabs where the 
edge is fixed (kips) 

Maximum resistajice of the equivalent system (kips) 

Fictitious maximum resistance (kips) 

Resistance of a structural element or system at time, t^ (kips)' 

Horizontal ultimate load resistance of a shear wall (kips) 

Yield resistance of the struct\ire (kips) 

Resistance in the elastic range (kips) 

Maximum resistance in the elastic range (kips) 

Maximum total resistance in the elastic range - "A" portion of 
two-way slabs (kips) 

Maximum total resistance in the elastic range - "B" portion of 
two-way slabs (kips) 

Resistance as a function of time 

Resistance of the g story columns as a function of time 


Time variation of the resistance of the equivalent slngle- 
degree-of-freedom dynamic syst^ for the g"!”^ story 

Dosage of nuclear radiation without shielding (roentgens) 

Resistance ais a function of displacement 

Radius of gyration of section (in.) 

Resistance per unit length of a beam or per vmit area of slab 
Ratio of web reinforcement » A^^/bs 

Distance from point of explosion (ft) 

Roentgens 

Ratio of steel reinforcing placed perpendicular to the steel beam 
(in composite construction) in excess of that required to carry 
the slab bending stresses 

bb 

Average value of the resistance of the g story when the relative 
displacement between the g^^ and the g-1^^ story is negative 


S Section modulus 

SE Strain energy absorption (kip-ft) 

(SE) Strain energy absorption of actual structirre (kip-ft) 

(SE)g Strain energy absorption of equivalent structvure (kip-ft) 

S' Section modulus about the weak axis 

s Coordinate axis at 45° angle to x axis in x-y plane 

Spacing of stirrups and spacing of ties in reinforced concrete 
members (in.) 
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Symbols 


I 

T 




Distance along s axis to centroid of loading (ft) 

Distance along s axis to centroid of inertia force (ft) 

Duration of the external load (sec) 

Resultant tension force (kips) 

Fxmdaniental period of vibration for complete spherical shell of 
centroid thickness (sec) 

Natiural period of oscillation or fundamental period (sec) 

Lowest natiiral period of circular arch with pinned or fixed ends 
Natural period of the i mode of oscillation (sec) 

Time or rise of a step-load (sec) 

Tp..., T Time duration of the external loads on the first, second, 

and g floors of a multi-story building (sec) 

T2 jj...T^ Time duration when the effective loaids on the first, 
second, and g^^ stories are negative (sec) 

Period of oscillation of a fictitious system representing the 
g'th story (sec) 

Thickness of concrete slabs (in.) 

Thickness of flange of beams (in.) 

Time measured after the arrival of the incident shbck wave (sec) 
Time variable 

Thickness of deep beam web 

Time of arrival, or time required for the shock wave to travel 
from the point of explosion to the chosen location (sec) 

Average flange thickness of standard steel beam (in.) 

Time required for overpressure on the rear face of a closed 
rectangular structure to rise from zero to its maximum value (sec) 

Time required to clear the front face of a structixre from the 
reflection effects (sec) 

Time displacement factor; the time required for the shock front 
to travel from the frontmost element of a stnictvure to the point 
or siirface under consideration (sec) 

Time at which the limiting elastic deflection is reached 
Flange thickness of a WF steel beam (in.) 

Time at which the limiting elastic deflection of the g^^ story 
is reached 

Time required for the vortex generated at the front face of a 
structure to travel a distance L' across the structvure (sec) 

Time required for maximum displacement of element or structure 
to occur (sec) 

Duration of the positive phase of the incident shock wave (sec) 


XVIII 


Symbols 


EM IIIO- 3 U 5 -I+I 5 
15 Mar 57 


Time of rise; time required for the overpressure in the incident 
shock wave to rise from zero to its maximum value (sec) 

t Thickness of stiffener (in.) 

D 

t^ Time required for the shock to pass over the length of the 

structiire (sec) 

t^ Time of rise of an underground overpressure pulse (sec) 
t^ Web thickness of a steel beam or channel (in.) 

t Time required to reach yield point of material (sec) 

^1' ^2*‘'^n^ ^n+l sequence 

^Im' ..t Time at which the maximum absolute displacements of the 
first, second and floor masses are reached 

Time required for the overpressure on the exterior rear face of 
a rectangular structure with openings to rise from zero to its 
maximum value (sec) 

t^ Time required to clear the front face of a structizre with openings 

from reflection effects (sec) 

t^, t* Time at which the relative displacements and X , respectively, 

® are equal to zero ® 



Time immediately before t 

n 

Time immediately after t^ 

Time lag 

Time interval used in numerical analysis 


At » t - - t A time interval 
n ru-l n 

Velocity of the incident shock front (fps) 

U.^ Velocity of the shock front of the shock imve formed in the 
interior of a rectangular structxire with openings (fps) 

u Particle velocity in the incident shock wave (fps) 

Bond stress per unit of surface area of bar in reinforced con¬ 
crete design (psi) 

V Dynamic reaction (kips) 

Total shear (kips) 

Total vertical force on piles 

V^ Total dynamic reaction along one edge "a," two-xra.y slabs (kips) 

V Average vertical shear in length, L’ (kips) 

V_ Total dynamic reaction along one edge "b," two-way slabs (kips) 

V Total column load in flat slab design (kips) 

c 

V Maximum shear cai>acity of web of deep beam section 

m 

V Maximum vertical shear force (kips) 

max 
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Symbols 


V Vortex velocity (fps) 

Shear stress (psl) 

V, Dynamic shear yield strength 

u.y 4*1^ 

V = y (t \ Velocity of the g floor at time, t 

g,n n; , n 

V Maximum shearing stress (psi) 
zu 

Vn * Velocity at time, t^ 

v+ Velocity at t+ 

n •' n 


V- 

n 


Velocity at t^ 



Horizontal velocity of axis of rotation "0" 

Velocity of mass at time t^ 

Static shear yield strength of steel (psi) 
Average unit shear stress in concrete (psi) 


v(t) 

W 


W 

a 

W 

e 

W 


W 


gm 


W 


m 


Velocity as a function of time 

Total energy yield of an atomic homib expressed KT of TNT 
required for an equivalent total energy yield 
Weight (lbs) 

Total load on element of structure (kips) 

Work done (ft-lbs) 

Dynamic peak blast load 

Work done on the actual system (ft-kips) 

Work done on the equivalent system (ft-kips) 

Maximum vork done on the g floor mass (ft-kips) 

"til 

Absolute maximum work done on the g floor mass (ft-lbs) 

Maximum vork done on the equivalent system by the equivalent load 


Fictitious maximum work done on the equivalent system 

W(t) Work done as a function of time ^ 

w Uniformly distributed load (kips/ft) 

Length of channel shear connector (in.) 

Width of front face of rectangular building (ft) 

Wj^ Weighting factor for vertical and batter piles 

X Relative displacement in a story of a multi-story building (ft) 

4“ Vi 

Maximum relative displacement in the g story (ft) 

X Relative displacement of the mass g when the average external 

loads are applied statically 

X^, Relative displacement in the first, second, and g^^ stories 

of a multi-story building (ft), X = x - x , 

g g g-1 
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X Limiting elastic relative displacement for tlte g story column 

^ (ft) 

X Distance the incident shock wave travels after impinging on the 

frontmost element of a cylindrical segment or spherical dome to 
any point being considered 

Deflection of a structural element or system (ft) 

Distance of .any pile from the centroidal aocis 

Collapse deflection 

Xg Limiting elastic deflection (ft) 

X- Forced solution of a dynamic system for a given external load 

X Absolute displacement of the g floor at time t 

th 

X Maximum absolute displacement of the g floor (ft; 

th 

X X , X ^ Absolute displacements of the g floor at time, 

g,n-± g,n g,nfi t , and t ^ 

i)-l' n^ ml 

Xj^ Ifexifflum displacement (ft) 

X T, X , X - Displacement at time t t , and t ,, resnectively 
n-1 n' m-1 n-l-' n' ml' 

Xp Initial displacement (ft) 

Horizontal displacement of axis of rotation "0" 

.Xg Displacement of an elastic system subjected to the peak load B 

acting statically (ft) 

x(t ) Displacement of a mass at time t 

' ®® th 

X-, Xp, X Absolute displacement of the first, second, and g floors 

® of a multi-story building (ft) 

X Distance to centroid parallel to x axis (ft) 

Horizontal velocity of axis of rotation "0" 

jd. Horizontal acceleration of axis of rotation "0" 


Acceleration in the x direction (ft/ sec^) 

Displacement or deflection of a stinictural system (ft) 

Deflection of act\aal element (ft) 

Midspem deflection of actual element (ft) 

Midspan deflection (ft) 

Limiting elastic deflection (ft) 

Deflection of equivalent system (ft) 

Limiting deflection in the elasto-plastic range (ft) 

MAxlraum displacement (ft) 

Displacement at time, t^ (ft) 
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Symbols 




y 

y 

^cl 

y 

z 

Z' 

z 

a 


a 


n 


u 


Vertical distance from axis of rotation "0" to mass centroid of 
rotating mass, of structure 

Deflection at midspan of column strip - flat slab (ft) 

Deflection at midspan of middle strip - flat slab (ft)' 

Maximum elastic range deflection (ft) 

Maximum elasto-plastic range deflection (ft) 

Distance to the centroid parallel to y axis (ft) 

Vertical distance from axis of rotation "0" to mass centroid of 
total moving mass, m 

Acceleration in the y direction (ft/ sec^ ) 

Velocity in the y direction (ft/sec) 

Velocity of displacement of the center of the beam (ft/sec) 

Increment ^of displacement (ft) 

Plastic modulus of the cross section (in.^) 

Plastic modulus of the cross section in the weak direction 

Depth of the resultant compressive force in concrete tee beam 

Angle of incidence between the normal to the surface and the 
direction of propagation of the blast wave (degrees) 

Central angle of the arch 

Deflection coefficient for two-way slabs 

Design load ductility reduction factor 

Angular acceleration of structure about axis of rotation ’*0" 

Dimensionless ratio •- 0.5P- /l^*7 
Ductility ratio 

Maximum bending deflection due to the application of a unit load 
(ft/kip) 

Shear wall lateral deflection at first cracking 

Stib-ai-ea clearing factor for rectangular structures with openings 

Maximum shear deflection due to the application of a unit load 
(ft/kip) 

Shear wall lateral deflection at ultimate resistance 
Deflection at theoretical yield (ft) 

Angle of Inclination to the horizontal of a gabled roof (degrees) 
Parameter used to define the location of a point on the surface 
cf a cylindrical segment or n spherical dome (degrees) 

End rotation of structural element 
Angular displacement 
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Joint rotation at bottom of colxmm 
Rotation of beam at midspan (radians) 

Angular displacement of structure about 8 ucls of rotation " 0 " 
Joint rotation at top of colxann 

Viscosity of air in the incident shock wave (ib-sec/ft^) 
C!oefflcient of friction 

Poisson's ratio 

Mass density of air in the incident shock wave (Ib-sec^/in.^) 

/ 2/4 2/' 4^ 

Mass density of soil (lb-sec /in. , or lb-sec /ft j 

Critical buckling stress (psl) 

Parameter used to.define the location of a point on the surface 
of a spherical dome 

Phase angle in the transient response (radians) 

Internal friction angle (degrees) 

‘t/ll 

Phase angle of the 1 mode (radisuis) 

Unit weight of soil (ibs/cu. ft) 

Angular velocity of structure about axis of rotation "0" 

(radlana/sec) 
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ENGINSESINC AND DESIGN 

DBSIGH OP SIROCTDKES TO BESIST THE EFFECTS OF ATOMIC WEAPONS 
PRINCIPLES OF DYNAMIC ANALYSIS AND DESIGN 

raTRODacraON 

5-01 PURPOSE AND SCOPE. This manual Is one in a series issued for the 
guidance of engineers engaged in the design of permanent type military 
structures required to resist the effects of atomic weapons. It is appli¬ 
cable to all Corps of Engineers activities emd Installations responsible 
for the design of military construction. 

The material is based on the results of full scale atomic tests and 
analytical studies. The problem of designing structures to resist the ef¬ 
fects of atomic weapons is new and the methods of solution are still in the 
development stage. Continuing studies are in progress and supplemental 
material will be published as it is developed. 

The methods and procedures were developed throxigh the collaboration 
of many consultants and specialists. Much of the basic analytical work was 
done by the engineering firm of Ammann and Whitney, New York City, under 
contract with the Chief of Engineers. !nie Massachusetts Institute of Tech¬ 
nology was responsible, under another contract with the Chief of Engineers, 
for the conpllatlon of material and for the further study and development 
of design methods and procedures. 

It is requested that any errors.-and deficiencies noted and any sug- 
* gestions for improvement be transmitted^o HQDA (DAEN-MCE-D) ViASH DC 20314* 

5-02 REFERENCES. Manuals - Corps of Engineers - Engineering and Design, 
containing interrelated subject matter are listed as follows: 

DESIGN OF STRUCTURES TO RESIST THE EFFECTS 
OF ATOMIC WEAPONS 

IM 1110-345-413 Weapons Effects Data 

EM 1110-345-414 Strength of Materials and Structural Elements 
EM 1110-345-415 Brinciples of I)ynaraic Analysis and Design 
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EM 1110-345.416 
EM 1110-345-417 
EM 1110-345-418 
EM lllC-345-419 
EM 1110 - 345-420 
EM 1110-345-421 


Structural Elements Subjected to Dynamic Loads 

Single-Story Frame Buildings 

Multi-Story Frame Buildings 

Shear Wall Structures 

Arches and Domes 

Underground Structures 


a. References to Material in Other Manuals of This Series* In the 
text of this manual references are made to paragraphs, figures, equations 
and tables in the other manuals of this series in accordance with the 
number designations as they appear in these manuals. The first part of the 
designation which precedes either a dash, or a decimal point, identifies a 
particular manual in the series as shown in the table following. 


EM 

paragraph 

figure 

equation 

table 

1110-31+5-^13 

3- 

3. 

(3. ) 

3. 

1110-31+5-1+14 

4- 

4. 

) 

4. 

1110-31+5-1+15 

5- 

5. 

( 5 - ) 

5. 

1110-31+5-1+16 

6 - 

6 . 

( 6 . ) 

6 . 

1110-345-1+17 

7 - 

7. 

(7. ) 

7 - 

1110-345-418 

8 - 

3 . 

( 8 . ) 

8 . 

1110-345-419 

9 - 

9 - 

(9. ) 

9. 

1110-345-420 

10 - 

10 . 

( 10 . ) 

10 . 

1110-345-421 

11 - 

11 . 

( 11 . ) 

11 . 


b. Bibliography. A bibliography is given at the end of the text. 

Items in the bibliography are referenced in the text by numbers inclosed in 
brackets- 

c. List of Symbols. Definitions of the symbols used throughout this 
manual series are given in a lisi^ following the table of contents. 

5-03 RESCISSIONS. (Draft) EM 1110 - 345-415 (Part XXIII - The Design of 
Structures To Resist the Effects of Atomic Weapons, Chapter 5 - Principles 
of Dynamic Analysis.) 

5-04 BASIC PRINCIPLES USED IN DYNAMIC ANALYSIS AND DESIGN. Before discuss¬ 
ing the fundamental principles of dynamic analysis, it will be helpful to 
review briefly the principles used in the analysis of structures under static 
load. Two different methods are used either separately or concurrently in 
static analysis; one is based on the principle of equilibrium and the other 
is based on work done and internal energy considerations. 

Under the application of external loads, a given structure will be 
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deformed and internal forces will be set up in its members. In order to 
satisfy static equilibrium, the vector sum of all the external and internal 
forces acting on any free body portion of the structure must be equal to 
zero. For the equilibrium of the structure as a whole, the vector sum of 
all the external forces and the reaction of the foundation must also be 
equal to zero. This principle is used in the analysis of statically de¬ 
terminate structures. 

The method based on work done and energy considerations is sometimes 
used when it is necessary to determine the deformation of a structure. In 
this method, use is made of the fact that the deformation of the structijre 
causes the point of application of the external load to be displaced. The 
force then does work on the struct\ire. Meanwhile, because of the structural 
deformations, potential energy is stored in the structure in the form of 
strain energy. By the principle of consemration of energy, the work done 
by the external force and the energy stored in the members must be equal. 

In static analysis, siirplified methods such as the method of virtuial work 
and the method of the unit load are derived from the general principle of 
energy conservation. 

In the analysis of statically indeterminate structures, in addition 
to satisfying the equations of equilibrium, it is necessary to include a 
calculation of the deformation of the structiu'e in order to arrive at a 
complete solution of the internal forces in the structure. The methods 
based on energy considerations, such as the method of least work and the 
method based on Castigliano's theorems, are generally used. 

In the analysis of structures under dynamic loading, basically the 
same two methods are used, but the load changes ra.pidly with time and the 
acceleration, velocity, and hence, the inertial force and kinetic energy 
are of magnitudes requiring consideration. Thus, in addition to the in¬ 
ternal and external forces, the equation of equilibrium includes the in¬ 
ertial force, and the equation of dynamic equilibrium takes the form of 
Newton's equation of motion, which is given by equation (5-l)- 

(mass) (acceleration) = external force - internal force (5-l) 
As for the principle of conservation of energy, the work done must be equal 
to the sum of the kinetic energy and the strain energy, namely: 
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work done = kinetic energy + strain energy (5-2) 

The strain energy in equation (5*2) includes both the reversible elastic 
strain energy and the irreversible plastic strain energy. Thus the differ¬ 
ence between structures voider static and dynamic loads is the presence of 
inertial force [|(mass) (acceleration )3 in the equation of dynamic equi¬ 
librium and of kinetic energy in the equation of energy conservation. Both 
terms are related to the mass of the structure. Hence, the mass of the 
structure becomes a very in^rtant consideration in dynamic analysis. 

In the above discussion, the two methods of analysis are described 
as if they were independent of each other. This is not the case. The force 
and energy equations are inter-dependent and axe derived from each other. 
Although these two methods are not independent, the convenience of applying 
them to a given problem varies greatly depending on the particular problem 
vmder consideration. Generally speaking, the force equation is convenient 
for analysis of stnictures, and the energy equation is convenient for design. 
5-05 DCTAMICALLr EQUIVALEMT SYSTEMS. In dynamic analysis, there are only 
three quantities to be considered. These are, (l) the work done, (2) the 
strain energy, and ( 3 ) the kinetic energy. To evaluate the work done, the 
displacement at any point on the structure under external distributed or 
concentrated load must be known. The strain energy is eqvsal to the summa¬ 
tion of strain energies in all the structural elements, which may be in 
bending, conqjression, shear, or torsion. The kinetic energy Involves the 
energy of translation and rotation of all the masses of the structvcre. The 
actvial evaluation of these quantities for a given structure under dynamic 
load would be complicated. However, for practical problems this can be 
avoided by using appropriate assumptions. 

In order to simplify the problem, a giyen structure is replaced by 
a dynamically equivalent system. The distributed masses of the given 
structure are Ivimped together into a number of concentrated masses. The 
strain energy is assvimed to be stored in a number of weightless springs 
which do not have to behave elastically. Similarly, the distributed load 
is replaced by a number of concentrated loads acting on the concentrated 
msses. Hence, the equivalent system consists merely of a number of con¬ 
centrated nBsses Joined together by weightless springs and subjected to 
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concentrated loads which vary with time. This concentrated mass-spring- 
load system is defined as an equivalent dynamic system. 

The reduction of a given structxire to an equivalent dynamic system 
involves the principle of dyimmic similarity which is merely the require¬ 
ment that the work done, strain energy, and kinetic energy of the equiva¬ 
lent system must he identical respectively with those of the given struc¬ 
ture [ 2 , 7 ] • The details of this principle and illustrative examples are 
given in paragraphs 6 .I 3 throxigh 6.2k in EM 1110-3^5-^16• 

For the time being it is svtfficlent to state that any given struc- 
tijre, be it a multi-story building or a structural member such as a beam or 
column, can be approximated by an equivalent dynamic system for the purpose 
of dynamic analysis. The discussion in this manual is limited to the 
analysis and design of equivalent dynamic systems consisting of concentrated 
masses, springs, and loads. 

a. Basic lynamic System. The siiiq)lest dynamic system consists of a 
concentrated mass supported by a weightless spring and subjected to a con¬ 
centrated load. This is shown in flgvire 5*1* 



In figure 5*1> the following notation is used: 
f(t) m external applied load which varies with time 
in m IDSLSS 
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R a Internal resisting force which depends on the displacement 
X a displacement which also varies with time. 

The directions in which the variables f(t), R and x are considered to 
he positive are indicated by arrows in figure 5.1. 

Suppose at a given time^ t , the displacement, velocity, and ac¬ 
celeration of the mass are x , dx/dt, and d^x/dt^ respectively) then, 
the following expressions may be written: 


Internal resisting force a R « kx when the spring is elastic 


Inertial force a m 


Work 


done by the load a J* f(t) dx a f(t) ^ dt. 


Kinetic energy 




c 12 

Strain energy « J Rdx a — kx when the spring is elastic 


Newton’s equation of mption is: 


f(t) - R 


(5.3) 


and the equation of conservation of energy is 


jf f(t) ■J* = I m (g) + / Bax 


(S.lt) 


These are the two basic equations used in dynamic analysis. Before attempt¬ 
ing to evaluate the displacement, x , from these equations for' a given ex¬ 
ternal load, f(t) , a few concepts about the external load and the resist- 
emce function need clarification. 

( 1 ) External Load and Work Done. Load Duration and Load Impulse. 

A typical load applied to a dynamic system is shown in figure 5-2. This 
load acts over an interval of time, T , which is defined as the duration 
of the load. The laqjulse of the external load, denoted by H , is equal 
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to the area under the load-time c\jrve; namely; 
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H - I f(t)dt 

Jq 


(5.5) 


Expression for Work Done. The work 
done by the external load up to any time, t , 
is given by: 

DC 

( 5 - 6 : 



Time 

Figure o.2. A typical 
dynamic load 


In order to evaluate the work done, it is neces¬ 
sary first to determine the velocity, dx/dt-, which can be obtained from 
equation (5*3) ty integration. This is given in symbolic form as: 


^ / [f (t) - e] dt 

e/n 


( 5 . 7 ) 


Substituting equation (5*7) into equation (5*6), the expression for work 
done then is. 




(5.8) 


Equation (5"8) indicates that the work done by the load varies with time. 

A typical case is shown in figure 5-3a'‘ At time, t^^^ , when the displace¬ 
ment of the mass reaches its maximum value, the work done is also a maximum. 
Ihe maximum work done W by the external load is obtained when the in- 
tegral limits in equation (5-8) are from 0 to t : 


- jT “ f(t) j I [f (t) - r] dt dt 


( 5 . 9 ) 


Because of the presence of m and R in the Integrand, it is evident that 
t he n»t v<tnnm work done depends not only on the external load but also on the 
mass and resistance function, R , of the dynamic system* 

Equation ( 5 . 9 ) is the general expression for the maximum work done. 
For the cases shown in figures 5*313 and 5*3c, the external load is zero 
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(0) ♦m<T 


over the time interval from T to 

W is a constant 
T . Hence, for the case where 
, the maximum work done is 
given hy: 

.t 


t . The value of 
m 

after 

t > T 
m 




W t„>T 


! 


[f(t) - r] dtidt ( 5 - 10 ) 


(c) t„»T 


Pure Pulse Load. For the 

case shown in figure 5 -Sc, the value 

of t is much greater than T and 
in 

the value of internal resisting 
force, R in equation ( 5 .IO) is 
negligible compeured with the value 
of the external load over the entire 
load dxiration, T . Under this con¬ 
dition, the external load is defined 
as a pure pulse load. In equation 

( 5 . 10 ), if it is assumed that the resistance R is zero for the full load 
dTiration, the maximum work done by a pure pulse load is given by; 



Figure 5.3. External load, displacement function 
and work done curves 


“"-f 


f(t) dt 


dt * 


£ 

2 m 


( 5 . 11 ) 


where 


H 


I f(t)dt = the impulse of the external load. 

Jci 


A pure pulse load corresponds to an impact load. During impact, the 
internal resisting force and the strain energy are both ass-umed to be zero. 
After impact, the mass, m , acquires an initial velocity equal to H/m . 
The kinetic energy of the mass after impact is given by equation (5-11) 
which is also the maximum work done by the load. This equation indicates 
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that for a pure pulse load, the maoclmum work done depends on the load Im¬ 
pulse, hut is independent of the detailed shape of the external load. 

Fictitious Maximum Work Done. In many actual problems, the in¬ 
ternal resisting force, R , in the time interval from 0 to T is not 
snail enough to be.neglected. From equations (5-3) and (5.10), it can be 
seen that if the internal force in the time interval from 0 to T is 
neglected, the actioal acceleration, velocity, and displacement will always 
be larger than the corresponding values when the resistance is considered. 
Likewise, the internal force in equation (5-9) always acts to reduce the 
work done by the external load. In other words, the actual maximum work 
dope on a dynamic system will always be smaller than the maximum work done 
on the same system whe^e the load is considered to be a pure pulse load. 
Hence, the work done as computed from equation (5-11) is defined as the 
fictitious maximum work done, namely. 



(5.12) 


The maximum work done, , as computed from equation (5*9) is al¬ 
ways sms Her than . The ratio depends on m , R , and on the 

detailed shape of the external load-time curve. The work done ratio is an 
ingport'""^ non-dimensional parameter to be considered in the design of dy¬ 
namic systems. 

( 2 ) Spring Constant, Resistance and Strain Energy. From para¬ 
graph 5 - 05 a, it is seen that the function of the internal force is to re¬ 
sist the movement of the mass, hence the internal force is defined as the 
resistance of the structure* The variation of the resistance with dis¬ 
placement of the mass is expressed by a resistance function. 

In figure 5.1, the strain energy is assumed to be stored in a spring. 

For a linearly elastic spring, the resistance function is eqioal to kx , 

1 2 

and, the strain energy is ^ hx . For the general case, the strain energy 
is equal to Rdx . Then the spring constant, resistance fxinction, and 
strain energy are related to one another by the following two equations; 

aifl.R ( 5 . 13 ) 
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dx 


k 


(5.14) 


where 

SE » internal strain energies which include both the reversible 

elastic strain energy and the irreversible plastic strain energy 

H = resistance 

k = spring constant 

Equations (5*13) (5.l4) are obtained from the basic dynamic system shown 

in figure 5 . 1 . complicated cases when more than one displacement is 

necessary to specify the configuration of the dynamic system, these two 
equations still hold true except partial derivatives are used. 

The relationship between resistance and strain energy as expressed by 
equation ( 5 * 13 ) indicates that two different methods may.be used to evaluate 
the resistance function. One method is based on the direct evaltiation of 
internal resisting force. The other method is based on the evaluation of 
strain energy, and the resistance is then derived from eqviation ( 5 . 13 ). 

The resistance wh^n plotted against the displacement, x , usually 
takes one of the forms shown in figure 5.4a. For structures made of brittle 
materials, the resistance f-unction is generally as indicated by curve A. 

For structures made of ductile materials with marked yielding, the 



R 



Displocement 


( 0 ) Typicol Resistonce Functions (b) Ideolized Resistonce 

Functions 

Figure 5 . 4 . Resistance functions 
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resistance function is given by curve B. For curve C, the slope of the re¬ 
sistance function is negative for large values of x . This is generally 
not a property of the material used, but is due to: (l) the combination of 
axial and lateral load on a given member of a given structiare, and/or (2) 
the effect of elastic or inelastic instability. In evalmting the response 
of structures under dynamic load, if better accuracy in the result is de¬ 
sired, the resistance function as shown in figure 5.^a should be used. But 
the method of preliminao:y design is usmlly based on a number of design 
charts which are evaluated on the basis of idealized resistance functions 
such as shown in figure 5-^h where one or two straight lines are used to 
approximate the actxoal resistance function. The type of idealized resist¬ 
ance to be used is discussed in paragraph 5-0?. 

b. Degree of Freedom of Dynamic Systems. The basic dynamic system 
considered in the preceding peLragraph consists of one mass, one spring and 
a concentrated load. Hence, a single displacement variable, x , is suf¬ 
ficient to describe its motion. This basic system is called a single- 
degree-of-freedom dynamic system, or a single-degree dynamic system. 

For the case of a multi-story building, more than one displacement 
variable is needed to describe its motion. The equivalent system consists 
of several concentrated masses connected by springs. 


The degree of freedom of a dy¬ 
namic system is defined as the num¬ 
ber of independent displacement vari¬ 
ables needed to specify completely 
the configuration of the system. 

Based on this definition, the dy¬ 
namic systems shown in figures 5•5a 
and b have single, and two degrees 
of freedom, respectively. 

For ordinary building struc- 
tvires under lateral loads, the mass 
of the structure can be ass\amed to 
be concentrated at the floor levels 



(o) Single Degree of Freedom 


(b> Two Degrees 
of Freedom 


and at the roof. The framework be¬ 


tween floors can be assumed to be 
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Figure 5 . 5 . Single-and two-degree-of-freedom 
dynamic systems 
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weightless splines or flat springs. The degree of freedom of a huilding is 
therefore equal to the number of stories of the building. In certain spe¬ 
cial cases, the displacements of any two floors can be assumed to always 
bear a constant ratio- The degree of freedom is then reduced. A multi¬ 
story building may even be approximated by a single-degree dynamic system. 
However, the error involved in this simplification may be very large. It 
is not recoinmended unless it is only for a first approximatioh in analysis 
or design. 

Beams, columns, slabs, etc., which have a continuous distribution of 
mass as shown in figure 5•5a are actually infinite-degree dynamic systems. 
However, if it is assumed that the deflected shapes of the member at any 
two instants of time are geometrically similar to each other, only one dis¬ 
placement is necessary (for example, the deflection at midspan). Thus, an 
infinite-degree-of-freedom system is reduced to a sirgle-degree dynamic 
system. The method of evaluating the parameters of the equivalent dynamic 
system for these members is discussed in paragraphs 6-13 through 6-24, 

EM 1110-345-416. 

c. Difficulties in Application of Dynamic Systems. A basic dif¬ 
ficulty in dynamic analysis results from the fact that the method of super¬ 
position is not directly applicable. If there are several loads applied 
simultaneously on a dynamic system with elastic springs, the maximum values 
of the displacements produced by individual loads do not generally occur at 
the same time. The maximum value of the displacement produced by all the 
loads is not equal to the s-ummation of the maximum values produced by in- 
dividml loads. Hence, in order to apply the method of superposition, the 
so-called phase differences in the dynamic responses to individual loads 
must be considered. Moreover, in some cases, the equivalent dynamic system 
of a structure may be non-linear, for example when a certain member is 
stressed beyond the yield point or when there is a vertical load as well as 
a lateral load acting on the structure. For non-linear dynamic systems, 
the method of superposition is not applicable at all. Hence, in general, 
all the loads acting on a dynamic system have to be considered 
simultaneously. 

In dynamic analysis, a given structure is approximated by its 
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equivalent dynamic system. Any structure is made up of a number of ele¬ 
ments such as beams, columns, girders, and structural connections. Each 
member is a distributed mass and possesses an infinite degree of freedom. 
Rigorously speaking, any structure is a multi-infinite degree dynamic 
system. However, by making certain assumptions and by applying the prin¬ 
ciples of dynamic similarity, a complicated structure is simplified to an 
equivalent dynamic system with a finite degree of freedom. Essentially, a 
distributed mass-load-spring system is replaced by a lijmped mass-load- 
spring system. The accuracy of the results obtained from the equivalent 
system as compared with those obtained from the actual structure depends 
entirely on the accToracy of the assumption used in the derivation of the 
equivalent system. In some cases, it may involve large error. It should 
be understood that the replacement of an actml structure by an equivalent 
dynamic system is one of the basic approximations used in dynamic analysis. 

Even for t]X)se cases when the equivalent system is acceptable as a 
good approximation for the given struct\jre, there are difficulties in the 
evaluation of its parameters. One typical example is the evaluation of the 
resistance of a one-story frame building. The Internal force which resists 
the lateral movement of the roof is equal to the total shear force at the 
top of the colxunns and this shear force is developed because of the bending 
and shear deformations of the columns. The bending deformation, hence the 
resistance, depends on the degree of fixity of the ends as well as on the 
size of the columns. The determination of the end fixity involves the size 
of girders, foundations, and the type of end connections. In practical 
problems, for the sake of simplicity, the two extreme conditions of fixity 
are generally used, that is, the end is ass\imed either hinged or fixed. 

For the determination of the resistance function of a member with stresses 
exceeding the elastic limit, it is necessary to consider the variation of 
the elastic limit with the rate of strain which is discussed in paragraph 
k-03&, EM 1110-3U5-kl3. There are still other effects to be considered 
such as the effect of dead weight and vertical load on the roof. Because 
of these factors, the actual resistance function is a complicated f-unction 
of the displacement. Although in the analysis of structxxres, these factors 
nay be partially taken into account in order to obtain a resistance function 
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whlcli is a closer approximation to the actml resistance function, idealized 
resistance fvmctions are always used in design- 

The same difficulty is encountered in the determination of the load as 
in the determination of the resistance function. The load is determined 
either directly or indirectly from the blast pressure which is generally a 
complicated function of time. Here again, the actual variation of load can 
be taken into account in analysis, but simplified loads are required for 
design. 

From the above discussion, it is seen that the parameters of the dy¬ 
namic system atre subject to iincertainty in-dileir evaluation. The analysis 
and design of a structure based on its equivalent dynamic system can only 
be carried out approximately. This is especially true in design since 
idealized resistances and sin^lified loads are used. Although rules are 
given at appropriate places which will help the designer to make reasonable 
assumptions and approximations, a basic understanding of the effect of dy- 
tiamic loads is necessary for the intelligent application of these rules to 
actml problems. 

In paragraphs 5-06 throiigh 5-08, methods of aralysis of single¬ 
degree dynamic systems are discussed- The results of a systematic anal¬ 
ysis of single-degree dynamic systems for several simplified loadings are 
plotted in non-dimensional design charts in paragraph 5-10. The applica¬ 
tion of these charts to the design of a single-degree dynamic system is 
discussed in paragraphs 5-12 throiigh 5-l4. The analysis and design of 
multi-degree dynamic systems are discussed in the remaining portion of 
the manual. 


Amnsis OF SmOLE-DECmEE-OF-FREEDOM SYSTEMS 

5-06 IMTOODUCTIOK. The analysis of single-degree dynamic systems consists 
of the evalmtlon of the displacement using Newton's equation of motion 
(equation 5-3). Two methods are generally used in solving this type of 
differential equation. One is the rigorous method in which the solution 
is obtained directly from the differential equation. This method is prac¬ 
tical only when both the load and the resistance function can be expressed 
in simple mathematical forms. The other method is the numerical method in 
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which the solution is obtained by a step-by-step integration. The nvonerical 
method is generally applicable for any type of load and resistance function. 
The two methods of analysis are discussed in paragraphs 5-0? and 5-08 re¬ 
spectively. A numerical example is given in paragraph 5-09 to illustrate 
the different methods of analysis. 

5-07 ANADfSIS BY RIGOROUS METHOD, a. General. The problem considered in 
this paragraph is the evaluation of the displacement using Newton’s equation 
of motion for a single-degree dynamic system as shown in figure 5*1* The 
eqviatlon of motion is: 

,2 

m ^ = f(t) - R ( 5 . 15 ) 

dt 

For a given load and a given resistance function, the procedure for 
evaluating the displacement is the same as that for solving a second order 
differential equation. First, the complementary function of the differen¬ 
tial equation is obtained by assuming the load to be zero. There are two 
constants of integration in the complementary function. Next the partic¬ 
ular Integral for the given load is evaluated- The complete solution for 
the displacement is the sum of the complementary function and particular 
integral which are defined respectively as the transient and the forced 
solutions in dynamic problems [ 2 , tI- From the known initial condi¬ 
tions, that is, the displacement and the velocity of the mass at the 
initiation of load, the constants of integration are evaluated. Thus an 
expression of the displacement as a function of time is obtained. It is 
obvious that the displacement function depends on the mass, m , the dy¬ 
namic load, f(t) , and the resistance function, R . 

The rigorous methdd of analysis is generally impractical for problems 
when either the dynamic load or the resistance 'is a complicated fimction. 
However, many important concepts related to the basic principle of dynamic 
analysis, and design can be illustrated when this method is applied to 
Idealized cases in which both the load and the resistance are of simple 
mathematical form. 

Three types of dynamic load of simple form are shown in figure 5*8 
together with the expressions for them. Each type of load is con^letely 
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apeclfled by two pexameters; the peak load B and the load dviration T . 



(0) Rectangular (b) Triangular (c) Step Load with 

Load Load Finite Rise Time, 

Tr 

Figure 5.6. Three idealized dynamic loads 

Three types of idealized resistance fvmctions as shown in figure 5*7 
are chosen for study. The dynamic systems characterized by these resistance 


R R R 



Displacement Displocement Displacement 


(o)Elostic (b) Completely Plastic (c) Elosto-Plostic 

Figure 5.7. Three types of resistance function 

functions are defined as: (a) linearly elastic, (b) completely plastic, 
and (c) eldsto-plastic. The expressions of these resistance functions are 
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also given in figure 5-7- It is noted that the linearly elastic and com¬ 
pletely plastic systems are only special cases of the elasto-plastic system. 
Ihe rigorous method of analysis of these basic types of single-degree dy¬ 
namic systems are given separately in the following three paragraphs. 

b. Linearly Elastic System. The resistance function of a linearly 
elastic system is always equal to kx . Substituting kx for R in 
equation ( 5 . 15 )» the equation of motion becomes: 

2 

m—^ = f(t) - kx ( 5 . 16 ) 

dt 

Referring to a textbook on Differential Equations, the complete solution of 
equation ( 5 .I 6 ) is given by: 

-jZ5^ + x^ (5*17) 

where 

= 2 nVm 7 ^ = natural period of oscillation 
x^ ® forced solution 
C s= constant of integration 

= phase angle (constant of integration) 

The sinusoidal term of eqmtion (5-1?) is the transient solution which is 

associated with the oscillatory nature of the dynamic system. The two con¬ 

stants of integration, .C and 0 , are determined from the requirement 
that the displacement and velocity of the mass must satisfy certain pre¬ 
scribed Initial conditions. For struct\u-es under dynamic load, the usual 
initial conditions are: 

at t = 0 , X « 0 , and dx/dt =0 ( 5 .I 8 ) 

Once C and 0 are determined, the displacement at any time is obtained. 

Equation (5.17) shows that, for any given dynamic load, the speed of 
response of the linearly elastic system depends on the natural period of 
oscillation, T^^ , in the transient solution. On the other hand, the 
rapidity of variation of the applied load depends on the load duration, T 


X = C sin 


("’'t 
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Thus the non-dimensional ratio, t/t , is expected to be an important 
parameter in dynamic analysis and design. 

As an illustrative exaii 5 )le, the case of a triangular external load 
given in figure 5* 6 b is considered. The expressions for the load are dif¬ 
ferent depending on whether t is greater or smaller than T . The cor¬ 
responding displacement functions are also different. 

For t < T , the expression for the load is: 


f(t) = B (1 - I) 

The forced solution for this load is: 



(5.19) 


(5.20) 


Substituting equation (5-20) 
conditions given by equation 
C and 0 can be evalviated. 


into equation ( 5 . 1 ?) and taking the initial 
( 5 - 18 ) into acco\int, the arbitrary constants 
These are given by: 



(5.21) 


f » tan"^ ^2it 


(5.22) 


Thus, in the time interval from 0 to T , the displacement is given by: 

where 

B 

Xg « « static deflection produced by the peak load, B . 

After t « T , the external load is zero. The forced solution is also zero. 
From the displacement and velocity at time, T , determined from equation 
( 5 * 23 )# the constants of integration of the transient solution after t > T 
can be determined. The displacement fxmctlon after t > T is given by: 
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Equations (5*23) and (5-24) shov that the non-dimensional displace¬ 
ment variable x/x at a given time, depends only on T/t . In figure 

S XI 

5«8a, the displacement of a linearly elastic system subjected to a tri¬ 


angular load is plotted versus time 
for three values of , namely, 

T/T^ = 0.2, 0 . 5 , and 2.0. The cor¬ 
responding cases of rectangular load 
and step load with rise time, , 
are plotted in figxjres 5* 8b and c. 
These curves show: (l) the maximoam 
displacement, (2) the time when the 
maximum displacement occttrs, and 

( 3 ) that the time for one coii^jlete 
oscillation is equal to T^ , and 

(4) the displacement at any time. 



The resistance of a completely 


plastic system is a constant equal 


to R . No practical system ever 
m 

behaves like a plastic system. But 


:^or a given system, if the maximum 


displacement x^^ is much larger 
■ hh«.n the niftTfinniin elastic deflection 


X , the elastic portion of the 

e 

resistance function can be neglected 



O 


(q) Triongular Lood 



(b| Rectongutcif Lood 



TImt Rotio 


(c) Stop Lood with Ri«# Timo 

figwrc 5 . 8 . The displccement functions of elnstic 
systems subjected to three types of external load 
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and the system treated as a completely plastic system. 

Substituting for R in equation (5.15)> the eqviation of motion 
of a completely plastic system is: 


mi^ = f(t) - R 


(5-25) 


Once f(t) is given, the displacement function can be evaluated directly 
by integration. 

For the case of a rectangular load given in figirre 5.6a, the expres¬ 
sions of the load are: 

f(t) = B for t < T ( 5 . 26 ) 

f(t) = 0 for t > T ( 5 . 27 ) 

Substituting eqmtions ( 5 . 26 ) and ( 5 . 27 ) into equation (5*25), integrating 
and taking equation ( 5 - 18 ) into account, the displacement functions of the 
mass for a rectangular load are given by: 


for t < T 




(5.28) 


for T < t < t 




where 


1 BT 
■ s “ 2 m 


t = the time when the maximum displacement x is reached, 
m m 


Substituting equation (5-28) into equation (5.6) and integrating, 
the work done by the external load up to any time, t , which is less 
than T can be obtained. This is given by equation (5*30). 


where 




(5.30) 


W - o — (see equation 5-12) 
•p 2 m 

_ 1 (BT)^ 

~ 2 m 
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Ite maximum work done by the load occurs at time, T , and Is given by; 


W 


m 


W 



(5-31) 


The preceding four equations show that the non-dimensional q.uantity 

is an in^xartant parameter in the analysis and design of completely 

plastic systems. In figure 5-9a-» the displacement and the work done by the 

load are plotted versus time for a rectangular load with B *= . The 

corresponding functions for a triangular load are plotted in figure 5*9b* 

These curves show: (l) the displacement function for a completely plastic 

system is non-oscillatory; (2) the maximum displacement; (3) the maximum 

work done; and (4) the time t when the maximum displacement x_ occurs. 
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Time Rotio 

(a) Rectangulai:- Lood 



Time Ratio 

(b) Triangular Load 


Figure 5.9. The displacement and work done on a completely plastic system B >■ 2R^ 


d. Elasto-Plastic System. This is tlie dynamic system which is com¬ 
monly used to represent a structure as a whole, as well as its members, for 
the purpose of analysis and design of blast resistant structures. Under 
blast loads, it is often uneconomical to design a structure to behave 
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entirely in the elastic range. Many ductile construction materials are 
characterized by a large yield range which enables them to absorb a large 
amount of energy. Use of this property in design results in much lighter 
nemberB than wo\ild be required should the energy absorption ability of the 
yielding region not be utilized. 

For an elasto-plastic system, the resistance function is given by: 

R = kx for X < Xg ( 5 - 32 a) 

R » R^ forx>Xe (5.32b) 

After the maximum displacement is reached, the dynamic system rebounds and 
the resistance function is given by; 


R - k [x^ - (*m - ^)] (^m - 

R = - R for (x - x\ > 2 x 

. m \ m / e 


( 5 . 32 c) 

(5.32d) 


For the rectangular and triangular loads shown in figure 5 . 6 , the expres- 
slonB of the load are different depending on whether t is greater or 
amaller than T . Hence, the six differential equations given below are 
needed to express the equation of motion. 


t < T, X < X 


t < T, X > 


■(&) ■ 


f(t) - kx 


= f(t) - R^ 


m ^ I d^X I 

T, X < x^, m 1—5 I = - kx 

' \atV 

/(i^x\ 

t>T, x>x^ ml — 7s I = - R 

' VatV 


T P' t > t , m 
m^ 


d^x 

t > t and t > T, m ^ = - k 
in c. 

dt 


( 5 . 33 a) 


( 5 . 33 b) 


(5.33c) 


( 5 . 33 d) 


(5.33e) 


h - (=V - *)] 


The first foTu: equations express the motion before the maximum displacement 
is reached and the last two equations express the motion when the dynamic 
system is in elastic rebound. 

The sequence of using these equations for the evaluation of the 
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displacement function of a given problem depends on: (l) the time, t^ , 
when the naximum elastic deflection is reached, and ( 2 ) the time, t^ , 
when the maximum displacement is reached. 

For the linearly elastic system, T/T^ is an important parameter 
while Rjjj/B is inqxsrtant for a plastic system. Both parametCTS are im¬ 
portant for an elasto-plastic system. As an example, the case of a tri¬ 
angular load as shown in figure 5-6b is considered. The ntimerical values 
of the non-dimensional parameters are assumed for illustrative purposes as: 

T 

= 0.3 

n 


Equation (5-338') is first used. Following the method described in para- 
graph 5-07b, the displacement function is given by: 


B 

* * k 




T 


=— = 


T 


n 



(5.31^) 


where 

» tan"^ 2Jt ^ » 62 °h' 
■‘‘n 


Equation (5-3^) is true so long as t < T and x < x 


Substituting 


_E _ c - 
k * S k 


into equation (5-3^)» the time, t^ , when the maximum elastic deflection, 

x^ , is reached can be evaluated. This is given by: 

t /T = 0.64 
e 

The displacement and velocity of the mass at time, t^ , acre respectively 
given by: 


"(‘e) - I 

''(s) ” ra 
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In the time interval, < T , equation (5-33h) should he used- 

The displacement function is obtained, by direct integration. This is given 
by: 

B 

*-k 

) [(ij - 3 5 (t) ♦ 2 (t) ] (5.35) 

At tlmej t as T y the displacement and velocity are respectively given by: 

x(T) = 0.831 2 

v(T) = 0.702 

In the time interval, T < t < , equation (5 •33d) should he used. The 

displacement■ in this time interval is: 

X -1 j^o.831 + 0.702(1 - 1^ _ (i - i^^j (5.36) 

From equation (5*3-)^ 'the maximum displacement, , is given hy: 

and the maximum displacement occurs at 

^ = 1.397 

After t > t^ , the dynamic system is in elastic rehoimd, and equa- 
ion (5-33f) should he used. The displacement after t = t^^^ is given hy; 

0.471 + I cos ^2« —( 5 . 37 ) 

e displacement function as expressed by equations (5.3^)> (5-35)^ (5-36)» 
( 5 . 37 ) is plotted in figure 5-10 in non-dimensional form. 
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B 


Figure S.IO^ Displacement and work done by a triangultu load on an elasto^plastic system 


Cf^^~0.3, C^-~-0.S 

^n 


B 


Once the displacement function is known, the work done by the ex¬ 
ternal load at any time, as expressed by equation (5.6), can be evaluated 
and this is also plotted in figure 5*10. The maximum work done by the ex¬ 
ternal load, which is eq.ual to the total strain energy at ^he msocimum dis¬ 
placement, is given by eq.uation ( 5 * 38 ). 

t 

"m - X \ fm - t) (5.33) 

As shown in figure 5*10, the ratio of ^or the given problem is 

equal to 0.8l4. 

5-08 AHAUSIS BY NUMERICAL METHODS, a. Need for Numerical Method. The 
rigorous method described in paragraph 5-0? is practical only for the cases 
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when both the load amd the resist¬ 
ance are simple mathematical expres¬ 
sions. However, for an actual prob¬ 
lem, the load may be quite compli¬ 
cated. A typical case is shown in 
figure 5-11* Fovir separate equations 
are needed to express f(t) . For 
this load applied to an elasto- 
t plastic system, the equations of mo¬ 
tion are expressed by eight separate 
differential equations. The solution 
Figure 5.11. A typical dynamic load of these equations by rigorous method 

is generally impractical. 

The rigorous method is fvirther complicated if the resistance function 
is not a siii5»le shape as shown in figure 5*7. Two practical cases are 
shown in figure 5*12. One is the idealized resistance function of a fixed- 
end beam, and the other is of the columns of a single story frame building 
when the vertical load is taken into consideration in the evaluation of . 
the resistance. 

For practical applications, when the load and resistance fiinctions 
are conq)lex, because of the impracticability of solving the rigorous 

R R 

* 

u 
c 
o 

« 

q: 


(a) Resistance Function (b) Resistance Function 

of 0 Fixed-end Beom of o One-Story Building 

Figure 5.12. Two practical resistance functions 
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differential equations, it is necessary to have other methods of analysis. 
Many methods have been developed. One is the method of numerical integra¬ 
tion which is discussed in this section. 

Perhaps it should he mentioned that these methods are by no means in¬ 
exact. Any degree of accuracy desired in the end result can be obtained. 
But, for the convenience of conqputation, these methods always embody sim¬ 
plifying assusqrtlons which affect the accxjracy of the results. The inac¬ 
curacies are not due to Inherent properties of the numerical methods, but 
rather due to the simplifying assumptions used in carrying them out. 


b. Basic Principles of Numerical Analysis. The differential equa¬ 
tion of motion (equation 5-15) expresses implicitly the acceleration of a 
dynamic system at any time, a » . it is known that the velocity 


and displacement can be obtained 
from the acceleration by inte¬ 
gration. The niamierlcal method 
of evaluating the displacement 
from equation (5.15) is called 
the method of numerical Inte¬ 
gration. 

The load, resistance, ac- 
cel^atlon, velocity, and dis¬ 
placement are plotted versus 
time to figure 5*13. for a typ¬ 
ical case. 



Suppose t^, tp 

t... ..t t.t .T*** is 

2 n - 1' n' n + 1 

a time sequence. The time in¬ 


terval from t^ to "tn + 1 1® 

denoted by dynamic 

is assumed to be initiated 


at t ■ t . The acceleration, 
o 

velocity, and displacement at 



'^0 denoted by a^, V^, figure 5.13. Load, resistance, acceleration, velocity 

respeOtively. If the and displacement versus time 
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acceleration in the time intearval, At^ , is represented by a(t) , the 
velocity and displacement at ^ are given reapectively by the follow¬ 
ing equations: 

% + 1 - ''n + X'“ " ^ “ (5-39) 

’‘n + 1 ■ *n ■'n (‘“n) * [/»(*) 4t] « (S-kO) 

These two equations indicate that the velocity and displacement at ^ 3 ^ 

can be obtained by extrapolation from the corresponding values at t^ 

once the acceleration in the time interval. At , is known. 

n 

In the analysis of structtnres under dynamic load, the velocity and 

displacement at the initiation of load are given as zero, that is: 

v^ ■ 0 and ■ 0 ( 5 .U 1 ) 

In applying equations (5*39) Mi<i (5«^0), the values of v^ and x^^ can be 

obtained provided that a(t) is known in the time interval from t^ to 

t^ . From the values of v^ and , the values of v^^ and Xg can be 

obtained. Ohis process can be continued until the values of v and x 

n n 

for any values of n are obtained. Hence, in the above numerical proce¬ 
dure the displacement, x^ Xg, ^c^..., • • ®re evaluated by a step- 

by-step extrapolation method starting from the known initial conditions. 

Here the difference between the rigorous method and the nvunerical 
method can be perceived. In the rigorous method, the value of x is ob¬ 
tained as a continuous function of time while In the numerical method, the 
value of X is obtained at discrete values of time. In other words, the 
actual relationship between x emd t is given by a differential equa¬ 
tion, but in the numerical method, this relationship is replaced by a 
finite difference equation. As the time interval in the sequence of time 
is reduced, the finite difference equation approaches the actual differ¬ 
ential equation as a limit. It can, therefore, be concluded that the ac¬ 
curacy of the numerical method based on the finite difference equation for 
most problems depends almost entirely on the magnitude of the successive 
time intervals. 
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Many extrapolation forzmilas have been derived for the solution of the 

finite difference equation j]5] • For engineering application, sla^jlifled 

extrapolation for mul as are often used which are obtained by assuming a 

slnple acceleration-time relationship in any time interval. At . IVo 

n 

procedures are developed in paragraphs 5-08c and 5-08d. In one case the 
acceleration is assumed to be varying linearly with time. In the other 
ease, the acceleration is assumed to consist of a series of pure pulses. 

c. Linear Acceleration Extrapolation ttethod. In this method, the 
acceleration, a(t) , in any time interval, At^ , as shown in figure 5* 14a, 
is assumed to vary linearly with time as shown in figure 5»l4b. 



Time Time 


(a) Actual Acceleration Curve (b) Approxlmote Acceleration 

Curve 

Figure 5.14. Comparison of actual and approximate acceleration curve 
used in linear acceleration extrapolation method 

In the time interval, At^ , a(t) is then given by: 

“(‘i - \ (* - "n) (5.42) 

irtiere At ■ t^ . , - t^ 
n n + 1 n 

Substituting equation (5-42) into equations (5-39) and (5-40), the recur¬ 
rence formulas for velocity and displacement are given respectively by: 
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( 5 .^^ 3 ) 


\^) (5.Wt) 

Once and ^ are known, these two equations are applicable to the 

evaluation of the velocity and displacement at t . from the known values 

n + 1 

at t^ 

According to definition, a(t) is given by: 

a(t) [f(t) -R] 

dt 

and ^ T is: 
n + i 


V * 1 - S ['(‘b + 1 ) - H. l)_ 


{5M) 


In equation (5*^5)> + l) known, but the resistance, + l) ^ 
may depend on ^ ^ which is still an unlsnown. Hence, in the applica¬ 
tion of equations (5*^3) (5«^), a trial and error procedure has to be 

r» 


adopted. A value of x 


h + 1 


be assumed to be 


X + 

n 


obtained. The value of x 


(At \ 

from which R(t 

\ W. 

V n + l| 


and a 


n + 1 


can 

are 


h + 1 


conqjuted value is compared with the assumed value, nils procedure is re¬ 
peated until the assumed and computed values of x ^ agree with each 

n JL 


other. 


Of course, in the plastic region, the resistance is a constant; 

^ can be obtained directly from equation (5.^^) and the trial and 
error procedure is not necessary. 


There are two basic characteristics of this method. First, it is, 
in general, a trial and error method, hence it may be inconvenient to carry 
out, and secondly, the acceleration is assumed to vary linearly with time, 
thus the variation of acceleration in the time interval, At^ , is par¬ 
tially taken into account. So this method is called a trial and error 



5-06d 


EM 1110-345-415 
15 Max 57 


linear acceleration extrapolation method. 

In the application of the above method, the successive time intervals, 
Atg,.. ‘At^ do not have to be eq,ual although they are chosen to 
be e<],ual for the convenience of conqputatlon. To determine the magnitude of 
tbe time Interval to be used for a given problem, a number of factors have 
to be considered. First, the time interval must be chosen such that within 
It the acceleration can be reasonably approximated by a straifdit line. 
Secondly, the convergence of the trial and error procedure must be rapid. 
Thirdly, the nuniber of steps necessary for the evaluation of the maximum 
displacement must be reasonable. For the first and second considerations, 
tbe time Interval should be small while for the third consideration, the 
time Interval should be as large as tolerable. Because of these conflict¬ 
ing requirements, it is found that a time Interval equal to about T^/6 is 
a good cdspromlse to be used in the elastic region, but a larger time in¬ 
terval can be tised in the plastic region provided that the external load in 
tbe time Interval is a reeisonably straight line. 

d. Acceleration Impulse Extrapolation Method. In this method, the 
actual acceleration curve shown in figure 5*15® is replaced by a train of 


a 



(o) Givtfi Acceleration 
Curve 





(t>) Acceierotion Pulses 
to Reploce the Given 
Acceleration Curve 



(c) Velocity 



Tim* 

(d) Displacement 


Figure 5.15, AcceleraUon impulse extrapo/aiion method 
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equally spaced impulses occurring at t^, t^, ^2 '’**^n * magnitude 

of tlie acceleration ingpulse at t^ is given by: 

l^t^^ • (At) (5«li6, 



QMs is shovn in figure 5'15t>« Since an ioipulse is applied at t^ , there 

is a discontinuity in the value of velocity at t^ . In the tine Interval 

from t to t , . the velocity is constant and the dlsplacenent varlee 
n n *r i. 

linearly with tine. The velocity and displacement thus obtained ere shoini 
in figures 5 '150 and 5 » 15 d. 

Suppose t^ and t^ Indicate the time Inmediateiy before and after 
the application of the impulse at t , and let v' . and v'*' indicate 
respectively the velocity at t“ and t^ j these two velocities are relatei 
by the following equation: 




(5A7) 


The relationship between 
are given by: 


n - 1 


and x^ , and between x_ and x_ ^ , 
n ' n n 4 1 


‘n - ’‘n - 1 - % > 

‘n + 1 - 


(5.W) 


ombinlng equations (5*^7) and (5-^8), the three successive diaplacenente 
re related by: 

Phis is the basic recurrence formula for the acceleration inpulae extrapo* 

Lation method. Once the values of x at t . and t are known, the 

n - 1 21 ' 

rabie at ^ ^ can be directly conqputed without resorting to a trial aad 
HrrOr procedure. 
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The evaluation of by the recurrence formu l a needs special con¬ 
sideration. In equation (5.*<'9), when n * 0, is given by: 

^1 “ ^o ■ ^-1 ®o 

Both x^ and x_^ are equal to zero, and the expression of x^ is sim¬ 
plified to: 

Xi * (At)^ ( 5 . 50 ) 


The correct value of a used in the computation is not equal to a'' as 

a'' o o 

indicated in figure 5 - 16 , curve a or curve b, but is given by: 



Time Time 

(a) (b) 

Figure 5 . 16 . Discontinuities in the acceleration curve 

When equation ( 5 . 51 ) is used in equation (5-50) for the evaluation of 
the value thus obtained is that given by equation (5*^) which is based on 
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the Linear Extrapolation Method. In order to avoid a trial an(^ error pro¬ 
cedure for the evaluation of , the value of a^ in equation ( 5 - 51 ) is 


taken to he 


instead of its exact value 


f (t J . R (t,) 


In figure 


5 . 16 a, if a^ = 0, -a^ is equal to . In figure 5 .l 6 h if a^ and 


are approximately equal to each other, the value 


may be used for a 


Equation (5*^9) is directly applicable for the evaluation of x^, x^... 

^ whenever the acceleration, a(t) , is a continuous curve. When* 
there is a discontinuity in the acceleration this equation is still appli¬ 
cable provided a modified value of the acceleration is used. This is il¬ 
lustrated in figures 5*l6a and b. In figure 5-l6a, there is a discontinu¬ 
ity at tj^ which occurs at the end of the time interval. At . Under this 
condition, the value of a^ used in the numerical procedure is the average 
value at discontinuity, n amely •+ for aj^ . In figure 5 . 16 b, the 

discontinuity occttrs within the time interval. At . The correct value for 


in this procedure is: 




( 5 . 52 ) 


For a given problem if there are one or two discontinuities, the 
method described above can be applied conveniently, fewever, if the number 
of discontinuities is large, it is more convenient to replace the given ac¬ 
celeration by a smooth cvnrve. 

The main advantage of the acceleration impulse extrapolation method 
is that the numerical computation is straightforward and involves no trial 
error procedures. Its main disadvantage is that the result is less 
accurate than the linear acceleration method. It has been proved that if 
the acceleration curve is a straight line or a series of straight lines 
with the same slope, the result obtained from the acceleration inqpulse 
method is exact. However, the actual acceleration curve may be rather 
complicated so that the result obtained by this method is usually only ap¬ 
proximate. A number of examples have been carried out, and the results 
indicate that the accviracy in the evaluation of the maximum displacement is 
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within 3^ provided the time interval used in the con^nitation is not greater 

than T /lO . 
n' 

Another disadvantage of this method is that any nianerical mistake in 
the conqputation is increased hy geometrical progression in subsequent 
steps. Hence, one must be extremely careful to avoid any numerical mistake 
especially in the early stages of computation. 

5 -O 9 NUMERICAL EXAMPLE. An exeimple is given to illustrate the actual ap¬ 
plication of the different methods of analysis. The same exan^jle is eval¬ 
uated by three different methods: (l) the rigorous method; ( 2 ) the linear 
acceleration extrapolation method; and ( 3 ) the acceleration impulse extrap¬ 
olation method. 

The external load and the resistance fimction for the example are 
shown in figures 5 * 17 a' and 5* 17 b respectively. 




Displacement 
(b) Resistance Function 


Figure 5.17. External load and resistance function for the illustrative example 

atsstoDed dimensional parameters of the dynamic system are 

2 

m « 2.5 kips-sec /ft 

k » 9860 kips/ft 

R « 750 kips 
in 

X - 0.0760 ft 
e 

T ■ 0.10 sec 
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B = 1000 kips 

=* 2ir-^= 0.10 sec 


The non-dimensional parameters are: 


R 

T_m 

"R ” B 


= 0.75, 




_T_ 

T 


= 1.00 


n 


In the rigorous method, the numerical values of the parameters are 
substituted into equation (5.23). The displacement in the time interval, 
0 < t < t^ is given by: 

^ - 60°57’') + 
n / 


X = 0.01012 


1.012 sin f 2 tc 






(5.53a) 


The time, t^ , when the displacement x is equal to the maximum elastic 
deflection, x^ , can be solved from equation (5.53a) and is equal to 0.022 
sec. The displacement and velocity of the mass at t^ are; 



x|t^^ = 0-0760 ft 

= 0.01012 (1.012) cos ^2jt- 80°57' 


5 .Ml- ft/sec 


Substitutii 3 g the numerical values of the parameters into equation (5.33b), 

and integrating, the displacement function in the time interval, 

t < t < t , is given by: 
e m 

X = 0.0760 + ^.kk (t - 0 . 022 ) -f 6 (t - 0 . 022 )^ - 666 (t - 0 . 022 )^ ( 5 . 53 ^) 


From equation (5.53b), it is found that the maximum displacement is 

X = 0.2825 ft 

m 

which occurs at t^ given by 

t = 0.077 sec 

HI 

The displacement computed from either equation (5.53a) or equation (5.531') 
is tabulated in column ^4- of table 5.1. 

The same example is evaluated by applying the two methods of numerical 
nalysls described in paragraph 5-08. The detailed comgputations of the 
.inear acceleration and acceleration impulse extrapolation methods are 
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given in tables 5>2 and 5*3 respectively. The displacements obtained by 
these three methods are tabulated in table 5 . 1 . It is seen that the error 
in using the numerical methods in the evaluation of the maximum displace¬ 
ment is less than 1 . 5 ^ 6 . 


Table 5.1. Comparison of the Results Obtained by 
Different Methods of Analysis 


Time 

Displacement Function x 

t 

(sec) 

Linear 

Acceleration 

Extr apo lat ion 
Method 

Acceleration 

Impulse 

Ebctr apo lat ion 

Method 

-1 

Rigorous 

Method 

0 

0 

0 

0 

0.01 

0.0181 

0.0200 

0.0186 

0.02 

0.0633 

0.0683 

0.0651 

0.03 

0.1182 

0.1216 

0.1203 

O.Ol)- 

0.1713 

0.1729 

0.1718 

0.05 

0.2184 

0.2182 

0.2184 

0.06 

0.2556 

0.2535 

0.2548 

0.07 

0.2787 

0.2748 

0.2785 

0-08 

0.2839 

0.2781 

t = 0.077 sec 

m 

0.09 


0.2594 

= 0.2825 

Maxiiiium 

Displace¬ 

ment 

0.2839 

0.2781 

0.2825 
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Table 5.2. Details of Computation by the Linear Acceleration Extrapolation Method 


1 

Trial 




a - /p - R V® 

Assumed 


Ml 

*n + 1 “ 


(kMqputed 


t 

n 




X . , 

^n + 1 



^n + 1 

-n.l 

(ft) 

1 

Ifo. 

(sec) 

(ft) 

(ft/sec) 

^ft/aec'^/ 

n + 1 
(ft) 

(klpa) 


' “ / “ t\^" 

(ft/B«‘) 

(ft/scc) 

B 

let 

0 

0 

0 

-jliOO 

0.0150 

900 


4282 

3.41 

0.0180 


2Dd 





0.0180 

900 

■SB 

4289 

3.44 

0.0181 

■ 

3rd 





0.0181 

900 


4286 

3-44 

0.0181 

B 

1st 

0.01 

0.0181 


•♦288 

0.0525 

800 

52B 

+113 

5.44 

0.0640 


2Dd 




0.0^0 

800 

631 

468 

5.22 

0.0632 


3rd 





0.0633 

800 

624 

+70 

5.23 

0.0633 

2 

0.02 

0.0632 

+5.23 

+70 

* 

700 

750 

-20 

5.48 

0,1182 

3 


0.03 

0.1182 

+5.W 

-20 

* 

600 

750 

-60 

45.06 

0.1713 

k 


0.04 

0.1713 

+5.08 

-60 

* 

500 

750 

-100 

44.28 

0.2184 

5 


0.05 

0.21B4 

Ak.2Q 

-100 

* 

400 

750 

-140 

+3.06 

0.2556 

6 


0.06 

0.2556 

+3-00 

! -lUO 

* 

300 

750 

-180 

+1.48 

0.2767 

7 


0.07 

0.2787 

-I.I48 

-180 

* 

200 

750 

-220 

-0.52 

0.2839 

« 

The dlsplaceoent Is larger than the maximum elastic deflection. H 

Is a ocmstant equal to 750 kips. 

The 


cosQutatlon of ^ and 

n X n 

^ becomes stral^tforward vltbout resorting to a trlaJL and error procedure. 

fcte: 

n 

- -1«» 

[1 - J t . 0.01 sec 










- 9660 X 

^ f or X < 
a 

m 0.0760 ft; R 

^ » 750 kips for x^ > 0.0760 ft 







* ‘x + iN 

\*x-\ 



)(i»)® 





+ 1* 

r 

a J 

-r 


' <> , 




Table 5.3. 


Details of Computation by Acceleration Impulse Extrapolation Method 



t 

P 

n 

R 

n 

P - R 
n n 

a m (t - R \/“ 

a (Afc)^ 

2x 

n 




(sec) 

(klpa) 

(klpa) 

(kips) 


(ft) 

(ft) 



0 

0 

+1000/2 

0 

+500 

4200 

40.02 

0 

0 

40.0200 

1 

40.01 

4900 

+197 

+703 

+283 

40.0203 

40.0400 

0 

40.0^ 

2 

+0.02 

4800 

4674 

+126 

+50.4 

40.0050 

40.1366 

40.0200 

40.1216 

3 

40. 03 

+700 

+750 

>50 

-20 

-0.0020 

40.2432 

40.0683 

40.1729 

4 

40,04 

4600 

+750 

-150 

-60 

-0.0060 

40.3458 . 

40.1216 

40.2i82 

5 

40.05 

+500 

+750 

-250 

-100 

-0.0100 

40.43^ 

40.1729 

40.2535 

6 

40.06 

4400 

+750 

-350 

-140 

-0.0140 

40.5070 

40.2182 

40.2748 

7 

40.07 

+300 

+750 

-450 

-180 

-O.OlflO 

40.5496 

••6.2535 

40.2^1 

8 

40.08 

4200 

+750 

-550 

-220 

-0.0220 

40.5562 

40.2748 

9 

40.09 









Note 

: P„.f(t 

.J.iooo[ 

‘-■(ei)] 








1 kx for ■ 

In n e 

0.0760 







B - 










1 R - 750 for X > X 









m 

n e 







X 

n + 1 * 

- X , + 

n - 1 

a„(At)® 








(At) n 0.01 sec 
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5-10 EESIGW CHARTS FOR SIMPLIFIED IQADINGg. a. General. Paragraph 5-10 
presents, in non-dimensional graphs, the results of a systematic analysis of 
single-degree dynamic systems for simplified loads. In order to these 
graphs simple and convenient for design use, only those cases vhich Involve 
no more than three independent non-dimensional parameters are considered. 

For these cases, the complete results of analysis can he clearly presented 
on a simple graph with one parameter as the abscissa, a second j^ameter as 
the ordinate, and the third as a running parameter. 

IVom the theory of dimensional analysis f l], the number of Independent 
non-dimensional parameters needed to express the functional relationship of 
all the variables of a given problem is given by equation ( 5 . 5 ^). 

N « n - m ( 5 * 5 ^) 

where 

N = number of non-dimensional parameters 
n m number of dimensional parameters 
m » number of fundamental dimensional quantities 
For a slngle-de^ee dy nami c system, there are three fundamental dimensional 
quantities: namely, mass, length, and time. If consideration is limited 
to problems with three or less non-dimensional parameters, the total number 
of dimensional variables can be no more than six. Of the six variables, 
one is the maximum displacement and one is the mass of the system. Hence, 
the load and the resistance must be conpletely specified by no more than 
four variables. If two variables, the spring constant, k , euid the plastic 
resistance, , are needed for, the resistance function, the dynamic load 
must be conpletely specified by no more than two variables. These are the 
peedc load, B , and the load duration, T '. It is for this reason that 
only very simple types of loads as shown in figure 5-6 are chosen as the 
basis for the preparation of design charts. Althou^ only three simplified 
inadH are considered, most actual loads can be replaced by one of these 
with sufficient accxuracy for preliminary design prxrposes. 

The non-dimensional graphs for the linearly elastic, completely 
plastic ftTid elasto-plastic systems sub^jected to the simplified loads shown 
in figure 5*6 are given separately in the following three sections. The 
non-dimensional parsuneters chosen for each case are also discussed. For 
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the linearly elastic system, all three types of loads are used hut only the 
rectangular and triangular loads are used for the contpletely plastic and 
elasto-plastic systems. 

h. Linearly Elastic System. For* a linearly elastic system subjected 
to simplified loads shown in figure 5 - 6 , the problem Involves five dimen¬ 
sional variables; m, k, B, and T Hence, only two non- 

dimensional parameters are needed in the presentation of the design charts 
and these parameters are chosen to be: 

X 

D.L.F. « — ■ Dynamic Load Factor 

C 

n 

^ere x is eg.ml to B/k which is the displacement produced in the 

s 

system when the peak load B is applied statically. 

The method described in paragraph 5 “ 07 b can be used for the evalua¬ 
tion of the dynamic load factor. Ihe displacement function for a given 
load is derived, from which the jnexlmum displacements, , for different 
values of C_ and the time, t , at which x 's occur, are obtained. 

The dynamic load factors for the rectangular load, the triangular 
load and the step load with rise time, , are plotted respectively in 
figures 5.I8, 5*19> and 5-21. The time ratios t^^ are also plotted in 
these figures for each loading. Ihe conposlte of figures 5*18 and 5*19 is 
given in figure 5*20 for the purpose of comparing the effects of the 
rectangular and triangular loads. ^ 

The dynamic load factor in this manual is defined as =2 vhere both 

^ s 

Xj^ and Xg are the displacements of the equivalent dynamic system. In 
the evaluation of the equivalent dynamic system, it is assumed that the 
deflected shapes of the struct\ire, for exaB5)le a simply-supported beam, at 
any two instants of time are geometrically similar to each other. Under 
this assumption, the slope and coivatiue at a given point of the beam de¬ 
pend only on the mid-span deflection of the beam Thus the dynamic load 
factors also give the ratio of the maximum Internal force developed at any 
point under dynamic load to the internal force developed at the same 
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point when the peak load., B , is applied statically. 

It should he noted that ordinary structures such as beams or frame 
buildings, involve a continuous distribution of mass. High modes of oscil¬ 
lation as well as the fundamental mode which is used in the evaluation of 
the single-degree dynamic system are present. As shown in appendix A, 

EM 1110-3^5-^16, although the higher modes have very little effect as far 
as.the mid-span deflection of beams is concerned, their effects on the 
bending moment and shear are quite large. Strictly speaking, the dynamic 
load factor based on a single-degree dynamic system, can only be referred 
to as the ratio of displacement but not as the ratio of internal forces. 

For the correct evaluation of internal forces such as bending moment and 
shear, high modes have to be considered and the structure must be replaced 
by a-multi-degree dynamic system. In practical application, for the sake 
of simplicity a single-degree equivalent system is still generally used to 
approximate the given structin'e. The same dynamic load factor is applied 
to the displacement as well as to the bending moment. The high modes of 
oscillation are sometimes partially taken into account in the evaluation 
of the shear and this is discussed in appendix A, EM 1110-3U5-4l6. 

Another point worth mentioning is that the dynamic load factor for 
the equivalent system of a stinictural member and the dynamic load factor 
for the -quivalent system of a structure as a whole have slightly different 
meanings. For the equivalent system of a member, an elastic beam for 
example, the elastic strain energy is stored all along the beam. If the 
shape of deflection curve remains the same, the internal force at any 
point of the beam is proportional to the mid-span deflection. Hence, the 
dynamic load factor gives the ratio of internal forces under dynamic and 
static loads at any point of the beam. But for the equivalent system of a 
structxare as a whole, the strain energy is assumed to be stored in a few 
of its members, for example, the colvmins of a one-story frame building with 
very stiff girders. For this case, the dynamic load factors give the ratio 
of internal forces Tinder dynamic and static loads at any point of those 
members in which the strain energy of deformation is stored. 

In the application of the dynamic load factor curves of figures 5*l8» 
5 . 19 ^ 8uid 5*21, the value of 0,^, is determined from the mass, m , the 
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spring constant, k , and the load dixration, T , of the equivalent dynamic 
system. Using , the D.L.F. and t^T are determined from figures 
5.18, 5*19> S’Did 5-21 for rectangular load, triangular load, and step load 
with rise time, respectively. 

The loading cxarves determined from the data in EM 1110-345-413 gen¬ 
erally do not conform to any of the simplified load shapes for which figures 
5.18, 5 - 19 > and 5-21 have been derived. However, with reasonable accuracy, 
these figures may be used for analysis and design where the given load 
curves have been approx im ated by one of these simplified load curves. The 
method of performing this approximation is explained in paragraph 5 - 13 - 

c. Completely Plastic System. For a completely plastic system sub¬ 
jected to the simplified loads shown in figure 5-6, the problem involves 
five dimensional variables: x , R , m, B, and T . Only two non- 

IQ in 

di^nsional parameters are needed to express the relationship of- these 
variables in chart form, and these parameters are: 

W 

^ = work done ratio 

P 

R 

^R “ ~ resistance-load ratio 


where 


W » the maximum work done by the external load and 

ss the fictitious maximum work done defined in equation (5-12) 

The work done ratios for the rectangular and triangular loads, eval¬ 
uated by the method described in paragraph 5-07c, are plotted in figures 
5.22 and 5-23 respectively. The time ratio, also plotted in 

these figures for each type of loading. These charts indicate that the 
work done ratio is independent of the duration of the load, but depends on 
the detail shape of the load; that is, whether the load is rectangular or 
triangular. 


In the application of figures 5*22 and 5*23, the value of W and 

Jr ■ 

Cjy are first determined from the peak load, B , the load duration, T , 
R 

the mass, m , and the plastic resistance, . Using is de¬ 

termined from these figures. The maximum displacement, x^^^ , is equal 
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wbere W is the product of C„. and W 
nr m m W p 

d. Elasto-Plaatlc System. There are six dimensional variables, x^, 
Mp k, B, and T, involved in the problem of an elasto-plastic system 

subjected to a rectangular or a triangular load. Hence, three non- 
dimensional parameters are needed and these are chosen to be 

T ^R “ B 

n 

For given values- of and , applying the method described in para¬ 

graph 5-07d, the displacement function for a given load can be obtained. 
Sid)stltutlng the maximum displacement into equation (5>38}> the maximum 
vork done by the external load can be evaluated.- This procedure vas applied 
to an elasto-plastic sys-tem subjected to rectangular and triangular loads. 
The results thus ob-talned are plotted in figures 5*24 and 5* 27* Because of 
the importance of the elasto-plastic system in actual application, four 
more graphs are plotted in figures 5-25, 5*26, 5*28, and 5-29* The first 
two are for rectangular loads and the last -two are for triangular loads. 

In figures 5.25 and 5*28, the ordinates are where x^ is the 

■nviiimni elastic deflection. These two graphs are useful in determining the 
■atr-tiwini displacement of a given structure for a given load. 

In both figures 5*26 and 5*29# the ordinate is t^^T where t^^ is 
the time when the maximum displacement is reached. 

In the application of -these graphs to an elasto-plastic system, the 
non-dimensional parameters, Cp and C^, are determined from the variables 
of the system. Using Cj^ and , the values of 

are determined from figures 5.24, 5 - 25 , and 5.26 respectively for rectan¬ 
gular load, for figures 5-27, 5-28, and 5-29 for triangular load. The 
maximum -work done, , the maximum displacement, x^^ , and the time, t^ , 

when X occurs can then be determined, 
m 

Althoiigh these graphs are derived for rectangular or triangular 
loads, they are generally applicable to loads which do not conform to 
these 8 lii 5 >lifled shapes. For this condition, it is necessary to idealize 
the given load and this is discussed in paragraph 5-13* 


W 


m 
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Ct = t/t„ 

Figure 5.19. Dynamic load factor and tJT cwrves for linearly 
elastic system, triangular load 
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Ct=T/T„ 

Figure 5.20. Composite dynamic load factor and tJT curves for 
linearly elastic system, rectangfilar and triangular load 
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Figure 5.21. Dynamic load factor and tJT curves for lihearly 
elastic system, step load with finite rise time, 
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Figure 5.22. Work done ratio and t^^/T curves for completely 
plastic system, rectangular load 
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0.1 1.0 10 40 

Cj = T/Tn 

Figure 5.24, lH/ork done ratio for elastO’-plastic system,, 
rectangular load 
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Figure 5.25. 



curves for elasto^plastic system, 
rectangular load 
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Ct = T/Tn 



0.1 1.0 10 40 

Ct = T/Tn 

Figure 5.26* j curves for elasto-plastic system, 
rectan^lar load 
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Figure 5 . 27 . ^ork done ratio for elas to-plastic system, 
triangular load 


53 








EM 1110-345 4 .5 
15 Mar 57 


5-lOd 



Fiffire 5 . 28 . 


■^j^i-curves for elasto-plastic system, 
* triangular load 
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e. Discussion of the Design Charts. In the preceding paragraphs (b, 
c and d) the response of single-degree dynamic systems subjected to simpli¬ 
fied loads are presented in non-dimensional graphs. In actual application, 
it would be helpful to know how the response of dynamic systems varies with 
the peak value, shape and dijration of the dynamic load and how the response 
is affected by the mass, spring constant and plastic resistance of the dy¬ 
namic system. In appendix A, following text, a quantitative discussion of 
the effects of small increments of the variables, namely: T, m, k, B, 
and , on the response are presented in detail and these effects are ex¬ 
pressed in terms of percentage increment ratios. In this paragraph, only a 
stumnary of the results derived in appendix A is described. 

First, the response of a linearly elastic system is considered. The- 
dynamic load factor curves for both rectangular and triangular loads are 
shown in figure 5-20. Depending on the value of , these two curvea are 
separated into three regions. In the first region, where > 3*0 , the 
D.L.F. for the rectangular load is constant; and the D.L.F. for the tri¬ 
angular load is practically constant. At a given , the D.L.F. ffiir both 
types of load are approximately equal to each other. Hence, in this region 
the maximum displacement, x^ , is: (l) directly proportional tua tiBe peak 
load, B; (2) inversely proportional to the spring constant, k , and 
(3) practically independent of the mass, m , the load duraticm, T , and 
the detail shape of the initially-peaked load. In the second region where 
< 0.3, the D.L.F. for the rectangular load at a given is about 

twice the corresjxmding value for a triangular load. This indicates that 
the area under the load-time curve (the load impulse defined by equation 
( 5 * 5 )) is important while the shape of the load is not important. In this 
region, the maximum displacement, > is: (l) directly proportional to 
the peak load, 3; (2) directly proportional to the load duration, T; 

(3) inversely proportional to the square root of the mass, m; and ( 4 ) in¬ 
versely proportional to the square root of the spring constant, k . In 
the third region, where 0.3 < < 3 * 0 ^ sill the variables m, B, T, k, 

X 

and the shape of the load are important. It can therefore be concluded 
that, so far as the maximum displacement of a linearly elastic system sub¬ 
jected to rectangular or triangulai- load is concerned, the peak load, B , 
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is always an important variable. The spring constant, k , is also ingwr- 
tant especially when is large. The mass, m , and the load duration, 

T , are important only when is small. The detail shape of the load 

is Important when 0,^, is between 0.3 and 3'0- 

Second, the response of a completely plastic system subjected to a 
rectangular or triangular load is considered. The maximum work done, , 
is: (1) inversely proportional to the mass, mj ( 2 ) proportional to the 
square of the load duration, Tj ( 3 ) approximately proportional to the 
squar^ of the peak load, Bj and ( 4 ) inversely proportional to the plastic 
resistance, B . The above conclusions are also true for the maximum dis- 
placement, , with the exception that the maximum displacement, x^ , 
is imrersely proportional to the square of the plastic resistance, . 

TOiearefare it can be concluded that all the variables B, T, m and R 
are always in^portant in the design of a completely plastic system. 

Third, the response of an elasto-plastic system subjected to a tri¬ 
angular load as shown in figures 5*27 and 5-28 is considered. Rererring 
to figure 5 - 27 , the work done ratio, , is practically constant when 
is smaller than 0 . 2 . In this region, the maximum work done is: pro¬ 
portional to the sqviare of the peak load; proportional to the square of 
the load duration; and inversely proportional to the mass. Hence, the peak 
load, the load duration and the mass are important variables for short 
duration load. For a given value of , the work done ratio is also 
practically a constant for large values of when is small. Thus 

the peak load, the load duration and the mass are also important variables 
for long duration loads provided is small. On the other hand for 

large values of and , the values of decrease with the in¬ 

crease of , hence the spring constant is important under these condi¬ 
tions. The value of Cy also decreases with the increase of , hence 
the plastic resistance, , is always important especially at large 
values of . 

Similarly, referring to figure 5 ' 28 , it can be deduced that, in the 
determination of the maximum displacement, the peak load is by far the most 
important variable. The plastic resistance is important for large values 

of C_, provided 0.5 < C„ < 1.25. The mass and the load duration are 

T R 
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important for short duration loads and those two variables are also impor¬ 
tant for long duration loads provided is small. The spring constant 

is important for large values of provided is also large. 

DESIGN METHODS FOR SINGIE-DEGREE-OF-FREEDOM D'lNAMIC SYSTEMS 

5 "H IRTROPICTION. The following paragraphs deal with the design of stru 
tures or their members wh^ch can be represented by a single-degree dynamic 
system. 

For a given design problem, the dynamic load is always given. Of tt 
three variables of the dynamic system, namely, m, k, R , the mass can t 

'XU' 

estimated and hence is considered to be known. There are also the inde¬ 
pendent variable, time, and the dependent variable, displacement of the dj 
namic system, the magnitude of which may not exceed specified limits. 
Hence, the designer's problem is to design a single-degree dynamic system 
with a known mass subjected to a given load such that the maximum deflec¬ 
tion is equal to or smaller than a certain amount. 

There are two variables, k and R^ , of the dynamic system which 
must be determined to conq)lete the design. Tiieoretically, there should b€ 
an infinite number of combinations of k and R^ to satisfy the given d€ 
sign specifications. However, considerations of economy, type of con¬ 
struction, and standard structural practice limit the relative variation c 
the spring constant, k , and the ultimate resistance, R^ . Thus, once 
the type of construction has been established, the size or composition of 
the stiructural members may generally be \miquely determined. 

In paragraph 5 - 10 , non-dimensional design charts are presented for 
idealized loeids. Ihe methods of design, using these charts, are first die 
cussed in paragraph 5 - 12 . When the given load is different from the 
idealized loadings, the method of approximating the given load so that the 
design charts can be used is discussed in paragraph 5-13- 
5-12 DESIGN METHODS USCTG CHARTS - IDEALIZED LOADINGS. Three methods 
using the design charts are given in this paragraph for the design of 
single-degree dynamic systems. These are the energy method and the de¬ 
flection method for cases where plastic deformation is allowed, and the 
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method based upon the dynamic load factor for cases where only elastic de¬ 
formation is permitted. 

a. Design for Plastic Deformation - Energy Method. The energy 
method for the design of elasto-plastic systems is based upon the fact that 
the work done on the system up to any time is equal to the sum of the ki¬ 
netic and strain energies. At the time of maximum displacement, the ki¬ 
netic energy is zero and the work done by the external load is equal to the 
strain energies. 

When the dynamic load and the mass of the system are given, the fic¬ 
titious maximum work done on the system by the external load is given by 
equation (5-55)- 


W 

p “ 2m 


f(t)dt 


(5.55) 


where 


H as area \mder the load-time curve 

as BT for rectangular load, or BT/2 for trisingular load. 

Suppose, using the method which is described later in this section, 

a trial size of standard structural section is chosen; the value of k and 

R can be determined from the geometrical dimensions of the structure and 
m 

the data in EM 1110-345-414; then the following non-dimensional parameters 


are computed: 



(5.56a) 

(5.56b) 


For the values of and , the work done ratio, 

then obtained from the design charts. Figure 5.24 is used for rectangular 
load, and figure 5.27 is used for triangular load. The actual maximum work 
done on the system, , is then computed from the relation • 

The nwYlnnnn allowable strain energy to be stored.in the system is 
equal to the area under the resistance-displacement c\arve up to the maxi¬ 
mum displacement, x^^^ . For an elasto-plastic system as shown in 
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figure 5* 7c if the maximum allowable deflection is x , the maaciimun allow- 

m 

able strain energy storage in the system is 

\ 

Comparison of the maximum work done and the maximum allowable energy 
storage determines the suitability of the' selected size. 

Based upon the above discussion, a trial and error design procedure 
is set up for an elasto-plastic system as follows: 

Step 1. From the parameters B and T of the given load and the 
estimated mass, the fictitious maximum work done, , is computed from 
equation (5.55). 

Step 2. A value of Wis first assumed and an approximate 
value for calculated using from step 1. A first approx im ation of 

the required maximum plastic resistance, , is determined from equation 

( 5 . 57 ) using this approximate value for and neglecting x^ as 

follows: 


or 


where 


R 


m 



R 


m 


W + E 
m 



( 5 . 58 a) 

( 5 . 58 b) 


X is maximum allowable deflection, 
m 

Step 3- A. trial section is then chosen having a value of Rj^^ ap¬ 
proximately equal to the value given by equation ( 5 * 58 ). 

Step 4. The actual value of R^ , the spring constant, k , and the 
maximum elastic deflection, x^ , are then computed from the selected sec¬ 
tion and the dimensions of the men±)er. 

Step 5- Factors 0,^, and Cj^ as expressed in equations ( 5 . 56 a) and 
( 5 . 56 b) are computed. 

Step 6. Using the values of C^, and is determined from 

the given design chart. Figure 5-24 is used for rectangular load and 
figure 5.27 is used for triangular load. 

Step 7- Hie actual maximum external work done, , is obtained'as 
the product of (calculated from equation (5*5^)) sJid obtained 

from the chart. 
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B and T , similar charts can be prepared. The above procedure would 
still give the exact solution provided the value of were obtained from 

the design chart corresponding to the given type of loading. However, in 
actual problems, the loading is rather complicated and generally cannot be 
specified by two parameters. It is impractical to prepare design charts 
which will give the exact solution for any type of loading. Hence, the de¬ 
sign charts are only prepared for two basic types of loading. For an actual 
problem, the loading must be approximated by either a rectangular or a tri¬ 
angular shape. The method of approximating the given load so that the de¬ 
sign chart for rectangular or triangular load can be used is discussed in 
paragraph 5 - 13 * 

The application of the energy method of design to actual problems is 
given in EM 1110-345-416 and 4l7. In paragraph 6-11, EM 1110-345-416, a 
numerical exanple is given to illustrate the design of beams. 

b. Design for Plastic Deformation - Deflection Method. The deflec- 

.. .—. — . . . r,..., 

tion method of design is practically the same as the energy method with the 
exception that the non-dimensional parameter, * is used instead of 

the work done ratio, . The design procedxxre for the deflection method 
is as follows: 

Step 1. A value of » Ris assumed, and a first approximate 
value of the resistance, R , determined. 

HI 

\-CrB (5.59) 

Step 2. A trial section is chosen having a plastic resistance, R^ , 
approximately equal to the value given by eqmtion (5-59)' 

Step 3- The actual values of R^^ , the spring constant, *k , and the 
mRylTmim elastic deflection, x^ , are computed from the selected section 
and the dimensions of the member. The mass, m , of the system is also 
determined. 

Step 4. Factors and as expressed in equations ( 5 . 56 a) and 

( 5 . 56 b) are computed. 

Step 5 . Xj^Xg is determined from the given design charts. Figure 
5.25 is used for rectangular loads, and fig\xre 5-28 is used for triangular 




loads. 


5-12C EM 1110 - 345-^15 

15 Mar 57 

Step 6. The maximum deflection, , is obtained as the product of 

x^ and X /x . 
e nr e 

Step ?• Con^ffltrison of the computed maximum deflection, x , with 

ID. 

the design specification for the maximum allowable x determines the 

m 

suitability of the selected trial section. By assuming a new value of , 

the design proced\ire from steps 1 to 7 is repeated until a section which 

gives the closest agreement between the computed x and the allowable 

m 

X is obtained, 
m 

c. Design for Elastic Deformation - Dynamic Load Factor Method. 

This design procedure is for the case in which the external load can be ap¬ 
proximated by either a rectangular or triangular shape or a step load with 
a finite rise time as covered by the design charts in figures 5-l8» 5 - 19 > 
and. 5*21. Methods for approximating a given load to conform to one of 
these sinqple types are covered in paragraph 5 - 13 * 

Step 1. A value of dynamic load factor is first assumed. 

Step 2. A trial section is selected having a maximum resistance, 

, approximately equal to 

R (D.L.P. )B 

where 

D.L.F. is the assxmied dynamic load factor. 

Step 3 . Thf» spring constant, k , is determined as described in 

paragraph b above, and the maximum elastic deflection is determined for 

various type elements as described in paragraphs 6 -I 3 through 6-24 in 

EM 1110 - 345 - 416 . ^ 

Step 4. The time ratio, « T/T^ “ ^ VI determined. 

Step 5 . Using , a revised value fop the dynamic load factor is 

determined from figures 5 .I 8 or 5.19 for rectangular or triangular loads 

respectively, or from figure 5*21 for step loads with a finite rise time. 

Step 6. The static deflection, x^ , is determined by x^ = B/k . 

Step 7* xnaximum deflection, , is obtained as the product of 

X and the revised dynamic load factor. 

8 

Step 8. Comparison of x^ and x^ determines the suitability of 
the selected trial section. The design procedure from steps 1 to 8 is re¬ 
peated \intil a section which gives the closest agreement between x^ 
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and X is ©"btained. If it is desired to limit the naximum deflection to 
e 

less than the maximum elastic deflection, a fictitious equal to the 

desired maximum deflection can be used. The maximum resistance of the sec¬ 
tion, R , vould he reduced in the same ratio as x to correspond with 
the smaller design deflection. 

An alternate dynamic load factor method may he used if desired, in 
which stresses are used rather than deflections. To use this approach, 
proceed as above but omit computation of x , x , and x • For step 8 
compare R^^^ of section with (D.L.F. )B to determine suitability of the 
selected trial section. By either approach, it is necessary to determine 
tjj^ and check the load shape as discussed in paragraph 5-13 • 

The application of this method to the design of elastic beams and 
one-story buildings is illustrated by numerical examples in EM 1110-345-^17- 
5 -I 3 METHOD OF ARPROXIMATIIC A GIVEH IDAD. The design methods given in 
the previous paragraph are based upon an idealized load acting on a dynamic 
system with an idealized resistance function. For an acttial problem, the 
resistance fxonction is very close to one of the idealized forms for an 
elastic or an elasto-plastic system; however, the given load may be radi¬ 
cally different from either a rectangular load or a triangular load or a 
step load with a finite rise time. Hence, the first step in the design is 
the replacement of the given load by one of the idealized loads for which 
the design charts are pxrovided. 

The basic consideration in replacing the given loading by an ideal¬ 
ized loeiding shape is that the maximum displacement produced by both must 
be equal. This replacement of load lasually introduces a certain amount of 
error and is one of the more difficult parts of the whole design problem. 

The amo\mt of error depends largely upon how well the loading is approxi¬ 
mated which depends entirely upon the judgement of the designer. Mb gen¬ 
eral rules can be given for replacing a given load by a siiiqE>lifled load; 
however, the following is set up as a guide. 

Suppose for a given problem, the load is shown in figure 5.30a. 

In the first approximation a rectemgular load as shown by the dotted curve 
my be used. Following the design procedures set up in paragraph 5-12, the 
ilues of and are obtained. Then the value of is obtained 
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from figure 5*l8» 5-22, or 5*26. This parameter is very in^rartajit in 

Judging how good the load approximation is. Since t is the time at 

m 

which occurs, the approximate load and the given load should (l) be 

very close to each other in the time interval of time from 0 to t , and 

in 

( 2 ) have about the same area in this time interval. 

Depending upon the value of t^ obtained from the first approxima¬ 
tion, the given load should'be replaced by one of the rectangular loads 
shown in figure 5* 30b and figure 5- 300 . However, if t^ is much larger 
than Tg as shown in figure 5 - 30 a, it is expected that the internal re¬ 
sisting force in the time interval from 0 to Tg is small. Based on the 
discussion in paragraph 5 - 05 a, for this case, the detail shape of the load 
is not important, and the given load can be replaced by a pure pulse load 
with equal total impulse. In this case, the value of is approximately 

equal to unity, hence 


First 



T T 2 t^T T 2 T T 2 

(a) (b) (c) 
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In figure 5* 30a, b, c, a rectangular load is used to replace the 
given load, hence the value of t^ is obtained from figures 5 * 18 # 5 *22, 
or 5-26. For the case vhen the given load is replaced by a triangular load 
(figures 5 * 30 d and 5 * 30 e), t^ may be obtained from figures 5 ' 19 » 5 * 23 ^ 
or 5 • 29 . When the given load is replaced by a load with a finite rise 
time (figure 5 . 30 f) t^ may be obtained from figure 5.21. 

5-l4 mSCUSSIOIT OF DESIGN METHODS, a. Accuracy of Design Methods. From 
the above discussion, it is seen that in applying the design proced\ares 
given in paragraph 5-12 to an actual problem, idealized load and resistance 
values have to be used to replace the actual values. It is expected that 
a certain amount of error will be involved. A final check of the actual 
design using the numerical methods of analysis described in paragraph 5-08 
is necessary. Therefore, this design method can be considered as a pre¬ 
liminary design procedure. The result obtained from the preliminary design 
must be checked by a numerical method of analysis using the actual load and 
resistance function. 

The maximum displacement evaluated by the numerical method may be 
different from the design specification. There are two reasons for this 
discrepancy. One is the fact that idealized loads and resistances are used 
in the preliminary design. The error from this source is usually less than 
'5^- The second is due to the use of standard structural sections. Only a 
limited nximber of sections are available. The difference in properties be¬ 
tween two adjacent sections may be more than yyf>. For a given design prob¬ 
lem, it is sometimes found that no standard section is available in order 
to satisfy the given design specifications. If a standard section close 
to the desired section is used, the maximum displacement may be different 
from the design specifications by more than 10^. However, a number of de¬ 
sign problems have been carried out and it has been found that the accu¬ 
racy of the preliminary design method, judging by the difference between 
the actual and the specified maximum displacements, is within 10^ for the 
design of one-story buildings, and 25 '^ for the design of structuual mem¬ 
bers, such as beams or slabs. The larger error in the design of beams and 
slabs is due to their shorter periods which make them more sensitive to the 
hape of the load, and to gaps in available structural sections (see 
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appendix A following the text of this manual). 

When designing systems with an elasto-plastic behavior, a small change 
in the max intum resistance of the section will cause proportionately larger 
changes in the maximum displacement (see appendix A). Because of this 
sensitivity, a 10^ error in the displacement represents a relatively 
smaller chauoge in the required maximum resistemce. Therefore, numerical 
analysis is not always justified in checking the design of an elasto- 
plastic system. 

h. Comparison of Different Methods of Design. In paragraph 5-12, 
three different methods of design are given. The energy and deflection 
methods.are applicable to the design for plastic deformation while the dy¬ 
namic load factor method is applicable only to the design of elastic 
systems. 

In both the deflection and the dynamic load factor methods, the non- 
dimensional parameter used is a ratio of deflections, 

former, and xjj/Xg for the latter. Hence, the dynamic load factor method 
is actmlly only a special case of the deflection method. 

Basically, there are two different methods of design, the energy 
method, and the deflettion metliod. Ihe design procediires for these two 
methods as given in paragraph 5-12 are on a trial and error basis. After 
a trial section is selected, the amount of nuDierical work involved in de¬ 
termining the suitability of the selected section is about the same no 
matter which method is used. But the convergence of successive trials for 
a given problem depends on which method is used. 

In the design of one-story buildings which can be replaced by a 
single-degree dynamic system, the value of is of the order of unity; 

and the value of C_ varies from 0.25 to 1.5 for an elasto-plastic system 

X\ 

Tinder triangular load. In these ranges of* values for and , the 

iDaxiinum work done, , is not very sensitive to small variations in the 
spring constant, k , or the plastic resistance, (appendix A). In 

other words, the naximum work done is practically constant when different 
sections are used in a trial and error design procedure- Thus, in suc¬ 
cessive trials, the only modification is in the strain energy storage. If 
in one trial design the section is found to be unsatisfactory, a second 
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trial can be started by using another section with maximum strain energy 
storage approximately eq.ual to the work done obtained in the previous trial 
design. On the other hand, if the deflection method is used within the 
same ranges of values of and T/T^ as mentioned before, the maximum 
displacement is rather sensitive to any change in the value of the spring 
constant, k , or plastic resistance, . If a trial design is fo\md to 

be unsatisfactory, there is practically no guide as to what the next trial 
section should be. Hence, in the design of one-story buildings, the energy 
method of preliminary design is preferred to the deflection method. 

However, in the design of individual members such as beams or slabs 
under triangular load, the value of is usually greater than 3 the 

value of is usually greater than 1. In this range of values of 

and , the variation of W„ with any change in k or R_ is about of 
the same order of magnitude as the variation of x^ (appendix A of this 
manual). Hence, in the design of beams or slabs, the deflection method is 
as convenient as the energy method so far as the convergence of successive 
trials is concerned. 

So far the discussion has been limited to the design of a single¬ 
degree dynamic system. For a multi-degree dynamic system such as a multi¬ 
story building, the design problem is complicated by the inter-actions be¬ 
tween different stories. It is shown in paragraph 5-21c that, for design 
purposes, any story of a multi-story building can be replaced by an equiv¬ 
alent single-degree dynamic system. From this, a method of design based on 
energy considerations has been developed which gives reasonably good 
results. 

MALISIS OF MULTI-DEGREE-OF-FREEDOM ETRAMIC SYSTEMS 

5 -I 5 IMRODUCTION. In the following paragraphs, the general problem of 
analysis of miilti-degree-of-freedom dynamic systems is considered. The 
types of structures and structural members which are replaced by multi- 
degree dynamic systems for analysis are, for example: (l) beams under 
lateral load when higher modes as well as the fundamental mode are consid¬ 
ered (appendix A of EM 1110-345-4-16); (2) the effect of support movement on 
he response of a beam; and ( 3 ) multi-degree buildings under lateral loads. 
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In order to make the problems under discussion easier to be visual¬ 
ized, the terminology of multi-story buildings is used in this and the sub- 
8 e<iuent sections of this manual. It should be understood that the methods 
of analysis are applicable for any multi-degree dynamic system and are by 
no means limited to multi-story buildings. In jiaragraph 5-16, the equa¬ 
tions of motion are set up. The rigorous and numerical methods of analysis 
are given in paragraphs 5-17 and 5-18 respectively. 

5-16 equations of .MOTIOW. For the dynamic analysis of a multi-story build¬ 
ing of n stories acted upon by lateral blast loads, it is necessary, to re¬ 
place the structvure by an equivalent dynamic system. The distributed masses 
of the structure are luniped together at the floor and roof levels, giving a 
series of n lumped masses. The masses are connected to each other by 
weightless springs which simulate the lateral resistance of the structure. 
Similarly the distributed loads on the structtare are Ixiaqjed together at the 
lumped masses. Hence, the multi-story building of n stories is repre¬ 
sented by an n-th degree dynamic system with n concentrated masses con¬ 
nected by weightless springs and acted upon by concentrated loads at each 
of the masses. The variables under consideration are the displacements of 
these masses. 


The following notations are used; 

For masses: m^, m^-m^ _ m^, m^ ^ ^-m^ 

For absolute displacement of these masses: 


*1' - 1’ V 

For external load on these masses: 


X ,... .X 

g + 1 n 


^’o^Ct), f• • ’'f^ _ Q^(t), fg(t), fg ^(t).. •.f^(t) 

The equation of motion of mass, m , is given by equation ( 5 . 60 ): 


- f (t) - 


( 5 . 60 ) 


where R ^ is the total resisting force which is exerted by the columns 
gt 

and walls above and below upon the mass m . Depiending upon the type of 
construction, two cases are considered for the determination of • For 

exaniple, in the case of a shear wall structure, where the building is con¬ 
sidered as a vertical cantilever and the effects of over-all bending on the 


69 



IM U.l 0 - 3 k 3^13 
15 Mar 57 


5-l6a 


building is taken into consideration, the resistance function, R . , de- 

gt 

pends upon the displacements of all the n masses- In the case of frame 
structures where the effect of over-all bending of the building is neglected, 
is a function of only three absolute displacements, x x , and 

*g + 1 * 

a- Multi-Story Shear Wall Buildings. The resisting forces acting on 
any floor of a multi-story shear wall building depend on the displacements 
of all the floors. Figure 5-31 shows the absolute displacements of any two 

floors represented by two masses, 

, *^9 _, m^ and m^ of a multi-story build- 

I g 

I j ing. The absolute displacement of 

1 II D each mass is represented by x and 

I mg ^ - «gi ^ g 

I t—j—1 x^ respectively, their relative dis- 

I n * placement being given as x^ - x^ . 

I // Due to the difference in magnitudes 

j 1 of the displacements x^ and x^ , 

j ^ pair of internal forces, , as 

[ ^ IJ I shown in figure 5*31 are developed. 

1 Rni I ' I ! The masses m. and m are said to 

1 _£L_^ m. j i g 

I I be coupled together. If the coupling 

I -^— I 

I ^ 1 V - V < spring constant is k . • the expres- 

I ‘ *I I gl 

}-->+<— -H Sion for the resistance function, 

I , in the elastic range is given 

» by equation (5.61). 


Figure 5,3L Any two floors of a 
multistory building 


R , « k . 
gi gi 


i K - 


(5-61) 


The positive direction for R . as well as the positive directions for x. 

gl X 

an^ X are indicated by arrows in figure 5 •31* 
g 

For a shear wall building, the coupling of each floor mass to every 
other floor mass and the groimd must be considered. The general expression 
for the total resisting force, R^^ , acting on the g-th floor mass of a 
mul-'-t-story shear wall building is: 


y 






5 -l 6 a 


m 1110-3^5-^15 
15 Mar 57 


+-+ +-+ R = 5^ R 

fft go gl gi cn . - £ 


( 5 . 62 ) 


i = 0 


The design of multi-story shear vail buildings will generally be simplified 
to elastic analysis because of the small allowable lateral deflections and 
the complexity of the problem. For this case, the above expression may be 
written as; 

V * Vg v(='s - *1) vk ■ ^n) 

l5o"6l{=‘g-’'l) <5-63)* 

In equation ( 5 . 62 ), R is the resistance developed by the shear wall at 

the g-th floor due to the relative movement of m with respect to the 

g 

ground. Note that when i = g, R . * R = 0 . In equation (5*63), the 
grovind displacement is assumed to be zero. 

Substituting equation (5-62) into equation ( 5 . 6 O), the equation of 


motion for mass, m , is: 

g 


n -§ . f 


w- E 


(5.61^) 


i = 0 


In the elastic range, equation (5-64) becomes 


® dt^ 


= f (t) - 7 ^ k . /x - x.\ 

g^ ^ Z-# \ g y 


( 5 . 65 ) 


i = 0 


There is one equation of the same type as equation(5.64) for each floor 
mass (see figure 5-32). Hence, the equations of motion for an n-story 
building consist of n simultaneous second order differential equations. 

* IHie resistance function, R^^ , may also be expressed by; 

E, ‘ ■‘'gi 


i as 1 


“ V V V ^g3‘ • • • V " V 


and 
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Figure 5.32. Equivalent dynamic system for an n-story shear wall building 

For a multi-story shear wall building in the elastic range, the 
variables of the equivalent dynamic system can be evaluated based on methods 
which are given in detail in EM 1110-345-416 and EM 1110-345-419. 

Having determined the variables of the system, and the load acting on 
the structure, the analysis is carried out for the elastic range using 
eqmtions of the form 5 * 65 , as explained in paragraphs 5-17 and 5-18. 

As soon as the stress at any point of a shear wall for any story ex¬ 
ceeds the elastic limit, not only the internal resisting force in this 
story is affected, but also practically all of the coupling constants of th.e 
equivalent dynamic system are affected. Whenever the stress in any story 
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exceeds the elastic limit, the resistance fvmctions of the con^^lete struc¬ 
ture have to he re-congiuted. Hence, it would not be practical to attempt 
to analyze a multi-story shear wall building with stresses exceeding the 
elastic limit. The application of the above equation is therefore limited 
to the analysis of shear wall buildings in the elastic range. 

b. Multi-Story Frame Buildings. 

In frame buildings, most of the re- j 

slstances, R . , in equation ( 5 . 62 ) f - - . . ^ 

81 I fn(t)- 

are either zero or negligibly small ^ . 

except for the resistances, R . ^ Rn,n-i " Rn 

e> s+i, Y j 

and R ,, which are the resist- ^ Rn,n-i 

g, g-1 Tn-ilO— 

ances developed between adjacent -Rn-i,n -2 = Rn-i 

masses. In an actioal building, — ’‘g 

f =i=pRg+2,g + l =Rg+2 

these resistances are developed be- Tg+nrj ^j^g-nj 

tween the floors by the columns and ^ ^ Rg+i g- Rg+i 

walls. Because of the relatively * 

large stiffness of the frame ele- V" i 

ments in the plane of the floor in -Xg-i —*•/ j '^9,g-i ^g 

frame buildings, the lateral re- fg-i(t)— 

^ - p = p I 

stralnts are not so important be- —^ g-i,g-2 g*' 

’‘ii_^ 

tween non-adiacent floor levels. , p-i—, ' R 2.1 “ ^2 

f|(t;—|» j mi I 

Eliminating all resistances except r,^ = r. 


fg+l(t)- 


fg -1 (t) 


V .Rn,n-l = Rn 

* ^n-l,n-2 ' Rn-I 

=r**f^g+2,g + l = Rg+2 

]_ 

> Rg +l,g = Rg +1 


g,g-l - Rg 


^g-i,g-2 ■ Rg-I 


R 2 .I-R: 


Eliminating all resistances except "y r,q = r, 

those developed between adjacent Y 

masses, the equivalent system in 

figure 5*33 is obtained. This Figure 5.33. Equivalent dynamic system for 

. , . . j ® multistory frame buildine 

equivalent system, which is appli- ° 

cable for both elastic and elasto-plastic behavior, is used in the anal¬ 
ysis of multi-story frame buildings. 

For the equivalent system shown in figure 5*33> using the following 
simplified notations; 


R R ; 

g, g - 1 g 


g, g - 1 


the total resisting force, R . , on mass, , 


\ ■ ^g + 1 


( 5 . 66 ) 
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In the elastic range, equation (5-66) hecoines 


5-l6b 


^rrf- = ~ 

gfc E\ & g 


+ k _ /x - - X ^ 
g + 1^ g + 1 gj 


( 5 . 67 ) 


Substituting eqmtion (5-66) into equation ( 5 . 6 O), the equation of 


motion for mass, m , is given by 

c> 


# = f„(t) - R + R 

g ^+2 g" ' g g + 1 


( 5 . 68 ) 


In the elastic range, equation (5*68) becomes 


' if'g + 1 - ^g) <5.69) 


Rearranging terms in equation (5-69), the following recurrence formula 
expressing the relationship of any three consecutive displacement variables 


is obtained. This is given by 


-^2 1 

d X 

-k X 1 + ™ —o + (k + h . , \x - k ,x . = fjt) (5*70) 

g g - 1 g ^.2 ^g g + lyg g + lg + 1 g 


For an n-th degree dynamic system as shown in figure 5*33> there are 
n separate equations of the same type as equation (5.70). These n equa¬ 


tions are 


^ 2)^1 ■ ^^2 " 


“2 ,.2 


+ (kg + k^jxg - kgx^ - k^x^ - fg(t) 


( 5 . 71 ) 


'^gr^ - 1 - -g t 1-g -t 1 - ‘g' 


Ix-kx -,-k tX . 
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m 


d X 


n 


n 


dt 


k X 

n n 


k X ^ = f (t) 
n n - 1 ■' 


In equation (5*71), all the coupling springs are assximed to be in the 

elastic range. If a certain spring such as k has reached its plastic 

S 

resistance, R , this value should be used for R instead of 

The analysis of an n-th degree dynamic system involves the evaluation 
the n simultaneous differential equations. 

5-17 RIGOROUS METHOD OF ANAITSIS. In the rigorous method, the standard 

procedures of solving simultaneous differential equations are used. Load 

and resistance functions must be expressed in simple mathematical forms. 

The displacement of any mass, say x , is first evaluated. The remaining 

S 

x's are then determined using the recurrence formula given by equation 

( 5 - 70 ) for frame buildings, or equation ( 5 * 65 ) for shear wall buildings. 

Eliminating all the x's except x from the simultaneous equations 

6 

given by equation (5-7l) for. frame buildings or equation (5-65) for shear 
wall buildings, a linear differential equation of the (2n)th order in x 

o 

is obtained. Just like the case of single-degree dynamic systems, the 

solution of X of an n-th degree dynamic system consists of two parts. 

6 

One is the transient term and the other is the forced solution. Thus the 
complete solution of x can be expressed as 

tS 


X 


8 




sin 



+ X 


gf 


( 5 . 72 : 


where is the total displacement of the mass, m^; x^^ is the portioi 

of X which is referred to as the forced solution; and the term in 
g 

brackets is the portion of x^ referred to as the transient solution. 

There are 2n constants of integration in the complete solution. These 
constants are determined from the initial conditions that all the n massei 
are at rest at the initiation of the load. 

There are n oscillatory terms in the transient solution. The fre¬ 
quencies of oscillation are the n natural frequencies of the n-th degree 
dynamic system. The lowest frequency is defined as the fundamental fre¬ 
quency of the dynamic system. 
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Several other methods, based on the same principle as described 
aboye, but differing in detailed techniques are also commonly used in the 
analysis of multi-degree dynamic systems. In one method, the set of simul¬ 
taneous differential eqxaations (5-Tl) or ( 5 - 65 ) is first transformed into a 
set of algebraic equations by applying the principle of Laplace Transforma¬ 
tion and operational calculus The complete solution is then obtained 

by inverse transformation. In another method, the solution of the problem 
is star*ted from the recurrence formula, equation (5-TO), by applying the 
methods of solving finite difference equations [S]. 

There is another method which is radically different from the above 
methods. In this case, the principle of generalized coordinates is used 
[2, 6, 7 ]. An n-th degree dynamic system is changed to an equivalent systen 
with n uncoupled masses. Each mass is individually supported by a spring 
and acted upon by a generalized load. Thus the n-th degree dynamic system 
is changed to a system of n \mcoupled single-degree dynamic systems. 

^Biis method of generalized coordinates is extensively used in the study of 
structural members such as beams, columns, or slabs, when the high modes 
of oscillation are to be conpidered. The application of this method to the 
analysis of elastic beams is given in appendix A, EM 1110-3it-5-^l6, and to 
the analysis of multi-story frame buildings in reference 6. 

The difficulty in applying the rigorous method of analysis to multi- 
degree dynamic systems is evident. The problem involves the solution of : 
simulteuaeous differential equations of second order with 2n initial con¬ 
ditions. Suppose two separate equations are needed for each load and the 
dynamic system is elasto-plastic; then, as many as i»-n sets of equations 
are needed to express the equations of motion. In each set, there are n 
simultaneous equations. Thus, in general, the rigorous method of analysis 
is impractical for nrulti-degree dynamic systems, and numerical methods mus 
be used. The two numerical methods described in paragraph 5-08 for single 
degree dynamic systems are also applicable to multi-degree systems. These 
two methods are described separately in paragraph 5-18. 

5-18 MUMERICAL methods of AHALTSIS. a. Linear Acceleration Extrapolatic 
Method. The equation of motion for mass, m^ , is given by equation (5-6( 

.... O 

Suppose a denotes the acceleration of m : the expression for a as { 

6 00 

function of time is 
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(5.73) 


where is the total internal resisting force acting upon m at'any 

time, t . The expressions of for the shear wall building and frame 

building discussed in paragraph 5 -l 6 are given respectively by the follow¬ 
ing equations: 


Shear wall building, R 




(5.Tit) 


1=0 


Frame building, R . = R - R 

gt g g + 1 


' (5.75) 


If the acceleration a^ is assumed to be linear in the time interval 
from t^ to \ ^ (5.44), the relationship between the 

displacements of at t^ and 3 ^ is given by 

\. r (“) - 

( 5 - 76 ) 


where 


/dx \ 

(^g)t = t ^ \dF/t = t 


At = t - t 

n + 1 n 


Based upon equation (5.76), a trial and error procedure can be set up for 

the computation of x’s at t from the values of x*s at t^ . A set 

n + 1 n 

of n values of x,, x„,-x ,....x are assumed for the displacement of 

i 2 g'^ n 

each mass at time t - from which the resistance R's and the accelera- 

n + 1 

tion a*s at time t , are determined. The actual values of x’s at 
n -f 1 

t^ ^ are then computed from equation ( 5 . 76 ). Comparison of the computed 
and assumed "Values of x’s indicates the correctness of the ass\xmed values. 
This trial and error procedure is tedious to carry out since it 
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involves n variables. The computed values of all the n variables must 
agree respectively with the n assvimed values. On the other hand, this 
method usually gives accurate results. When all the coupling springs are 
in the plastic range where the resistances are constant, the computation 
becomes strai^tforward. Further details of this method are given in para¬ 
graph 5“08 where it is described for the analysis of a single-degree dy¬ 
namic system. 

b. Acceleration Impulse Extrapolation Method. The recurrence for¬ 
mula, equation (5"49), for the single-degree dynamic system is also appli¬ 
cable for multi-degree dynamic systems. For the displacement, x , of 

S 

mass, m , the recurrence formula is 
g 


where 


At 


n + 1 


- t = 
n 


n 


n 


- t 


1 



= the acceleration of mass, mg , at time 
the expression given by equation (5.73) 





The method of modifying the acceleration at any point of discontinuity 
of the acceleration cxirve is the same as in the case of single-degree dy- 
namic systems (see paragraph 5-08), The advantage and disadvantage of this 
method as described for the single-degree d^mamic system also hold for 
multi-degree systems. If the time interval (At) is chosen to be about 
T^/lO where T^ is the shortest period, of oscillation, the accuracy in 
the maxlmuin absolute displacement is about 3^* 

The application of the acceleration impulse extrapolation method to 
the analysis of multi-story frame buildings is shown in EM 1110-3^5-^18- 
The equivalent dynamic system for frame buildings is the case shown in 
figure 5-33- 

DESIGN OF MULTI-STORY FRAME BUILDINGS FOR PLASTIC DEFORMATION 

5-i9 INTRODUCTION. The problem considered in this and the following para¬ 
graph is the design of a multi-story frame building which is represented by 


78 



•5-19 


IM 1110-345-415 

15 Mar 57 

the equivalent dynamic system shown in figure 5.33- The terminologies of 
multi-story buildings are used in the discussion. Thus the concentrated 
masses are approximately the masses of different floors, the displacements 
of each mass are the displacement of each floor level, and the springs are 
the coltmms between axijacent floors. 

In the design problem, the floor masses can be estimated and the ex¬ 
ternal loads acting on the masses are given. The maximum allowable dis¬ 
placements between adjacent floors are prescribed. The problem is the se¬ 
lection of the coltunn sections for each story such that the relative dis¬ 
placements between adjacent floors do not exceed the prescribed mtoimum al¬ 
lowable values. 

One approach to the design problem is the preparation of non- 
dimensional design charts as in the case of single-degree dynamic systems. 
Once charts are available, the complete design can be carried out in a 
series of routine computations. A little reflection will indicate that, 
even for a two-story building, there are too many non-dimensional paraiii- 
eters to permit the preparation of practical design charts to cover all-pos¬ 
sible cases. Thus the design of multi-story buildings cannot be carried ' 
out in as simple and straightforward a manner as the design of single¬ 
degree dynamic systems. 

Another method of approach to the design problem is to first assiome 
the sizes of all the columns for the entire building and then check the 
adequacy of the assumed sizes by a numerical analysis of the whole struc¬ 
ture. Any readjustment of the column size is then made according to the 
results of the numerical analysis. The size adjustment is difficult be¬ 
cause the change of column size of any one-story affects the relative dis-- 
placements of the entire structure. If, in the numerical analysis, the 
relative displacements between several adjacent floors are found to be un¬ 
satisfactory, it is difficult to decide which column needs modification, 
and by what amo'unt the size of the column should be changed. The detailed 
numerical analysis may have to be carried out over and over again before a 
set of satisfactory column sizes is found. This method, therefore, is 
difficult and tedious to carry out, and is not recommended except for 
problems where a rational method of design has not been developed. 
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The design method recommended for multi-story frame buildings is based 
on a floor-by-floor iteration procedxire in which the colnmn size for the 
first, second, third, and on up to the top story is selected one after the 
other. Specifically, the outline of the design procedxire is as follows: 

Step 1 . The plastic resistance and the size of the columns for the 
first story are determined on the assumption that (l) the first floor is 
a single-degree dynamic system with a mass, m^ , subjected to an external 
load, f^(t) , and (2) all the remaining floors are absent. 

Step 2 . The plastic resistance, Rp^ , of the col-umns for the second 
story is determined by (l) incorporating the result of the preliminary de¬ 
sign of the first story, and (2) assumiiig that all those floors above the 
second floor are absent. 

Step 3 - The plastic resistance, R^j^^ , of the colvanns for the third 
story is determined by (l) incorixsrating the results of the preliminary de¬ 
sign of the first and second stories, and ( 2 ) assuming that all those 
floors above the third floor are absent. 

Step 4 . The above procediure is continued until the ultimate resist¬ 
ance, R f of the columns for the top story is determined. 

Step 5. The preliminary plastic resistances determined in step 1 
through *+ for each story are adjusted by talcing into account the effect on 
the columns of a given story by the col'uinns of the story which are directly 
above the given story. For example, the shearing force developed in the 
columns of the third story is considered in modifying the plastic resist¬ 
ance of the second story. 

Step 6 . The column sizes for different stories are then determined. 
The sizes of the selected columns should be based on the revised resistance 
values obtained in step 5 * lii "the determination of the size of the se¬ 
lected columns, the effect of vertical load on the plastic moment of the 
section should be taken into account. 

Step 7. A numerical analysis for the entire building is then car¬ 
ried out to determine the suitability of the selected columns. Any revi¬ 
sion of the column size can be made according to the result of the numer¬ 
ical analysis. 

There are two basic procedures in the above floor-by-floor iteration 
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design method. One is the preliminary determination of the resistance of 
the columns for any story after the column resistances of all those stories 
below this story have been determined. The other is the adjustment of the 
preliminary resistance thus determined for any story after the column rfe- 
sistances of all those stories above have been determined. The fundamental 
principles in carrying out these two basic procedures are discussed in 
paragraphs 5-21 and 5-22 respectively. The complete design procedxire for 
the design of multi-story buildings is given in detail in paragraph 5-23* 
me assxmrptions used in the floor-by-floor design procedure and the ac¬ 
curacy of this design method are summarized in paragraph 5-25- 

In paragraph 5-20, the absolute displacements, relative displace¬ 
ments and resistance functions of a three-story building are plotted versus 
time. These ciirves are used (l) to illustrate the general characteristics 
of multi-story buildings, and (2) to explain why a somewhat complicated 
procedure has to be used in design. 

5-20 DISPIACEMEETS PJSfD RESISTANCE FUUCTIONS OF A THREE-STORY FRAME 
BUILDING. A three-story frame building designed for blast resistance is 
illustrated in figure 5.34. The floor-by-floor resistances and displace¬ 
ments of this building under the design blast loading are plotted in 
figure 5.35. This building was designed based on the specification that 
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Figure 5,35. External loads, displacements and resistances of threfstory 
frame building as functions of time 


the maximum allowable relative displacement for any stoi’y should not exceed 
0.33 ft. Subscripts 1, 2, and 3 refer to the first, second, and third 
floors above the ground respectively. Hence, the third "floor" is the 

roof. 

Due to the differences in the floor masses and in the external loads 
on the masses, the relative displacements between and = 

Xg - x^ , and between and m^ , = x^ - x^ , are negative shortly 
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after the initiation of the load. The relative displacement, Xg , is zero 

at t^ and the relative displacement, , is zero at t^ . In the time 

interval from 0 to ^2 * resistance, Rg , is negative; and in the time 

interval from 0 to t’ , R^ is negative. The maximum elastic deflections 

for the first, second, and third stories are reached at time, t . t 

le ^ 2e * 

and t^^ respectively. 

The curves in figure 5-35 show three important characteristics which 
are generally true for any multi-story frame building. First, the ratio of 
the absolute displacements of any two floor masses varies over a rather 
large range. Hence, in design, it cannot be assumed that the absolute dis¬ 
placements of any two floor masses always bear a constant ratio. Secondly, 
the resistances, Rg , or R* are negative in a short interval of time 
after the initiation of load. The subsequent displacement of mg or m^ 
is considerably affected by the negative portion of the resistances. Hence, 
in design, it Is necessary to take the negative portion of the resistance 
function into account. Thirdly, the msocimum absolute displacements, , 

Xg^, and of m^, mg, and m^ do not occur at the same time. As 

shown in figure 5* 35b, x^, Xg^, and x^^^^ occur at t^, tg^, and t^^^ 
respectively, and t, < t_ < t_ . When m„ reaches its maximum dis- 

Im 2m 3® 2 

placement, x ^ , the floor mass, m^ , is moving backward due to elastic 
rebound. Hence, in the evaluation of the maximum relative displacement be¬ 
tween the first and second floors, it is necessary to consider the elastic 
rebound of floor mass, m^ . 

The floor-by-floor iteration design procedure presented in these para¬ 
graphs is intended to take these important characteristics into account. 

In the preliminary design of the columns for any story, the resistance as 
a function of time for the story must be estimated. The energy method of 
design given in paragraph 5-12 for single-degree dynamic systems, with 
slight modification can be applied to the determination of the required 
plastic resistance for the story. In the evaluation of the maoclmum allow¬ 
able strain energies of the columns, the maximum absolute displacements for 
masses, m^, nig, and m^ given in equation (5.78) must be used. 

(5.78.) 
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^2m 

(*lm ■ *le) 

+ X^ = X, + x^ 
2m Im 2m 

“^le 

( 5 . 78 b) 

bm = - 

’'ae) " ’'3m 

= X, + X^ + X, 

Im 2m 31*1 

- ^le - ^2e 

( 5 . 78 c) 


where ^3111 prescribed maximum relative displace¬ 
ments for the first, second, and third stories respectively; and Xg^ 

are respectively the maximum elastic relative deflections for the first sind 
second story columns. 

In equation (5-78), the maximum absolute displacements x^^, 

are expressed in terms of the maximum allowable relative displacements 
by assianing that the rebound of the lower stories equals their maximum 
elastic relative deflection. The maximum absolute displacements thus de¬ 
termined are used for the computation of the energy absorption by eqtiations 
given in paragraph 5 - 21 c. 

5-21 THE PRELIMINARY DESIGM OF COIIJMNS FOR AlifY STORY, a. Effective Load 
and the Equations of Motion in Terms of Relative Displacements. In para¬ 
graph 5 - 19 > the general procedure of iteration method of design is given. 
The basic principle of the preliminary design of the column for any stoiry 
after all the stories below this story have been designed is given in this 
paragraph. 



Hie equivalent dynamic 
system of a multi-story frame 
building is given in figure 
5* 33- Any three consecutive 
floors, g - 1, g, and g + 1, 
are shown in figure 5 - 36 . 

The equations of motion of 
these three floor masses are: 


m 


g + 1 


8+1 _ f (t) - 

dt B + J- 


R 


+ 1 g 2 


Figure 5»36. Three consecutive floors of a 
multi^story frame building 


( 5 . 79 a) 
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= f (t) - R + R 


g + 1 


(5.79b) 


U T 

g - 1 


= ,(t) - R . + R 

g - 1 '- ' g - 1 g 


(5.79c) 


In the preliminary design of the columns for story g , it is assumed that 

all the stories above this story are absent. Neglecting R ^ in equa- 

g + 1 

tion ( 5 - 79 b)> the eq'uation of motion for m becomes: 


n -1 = f (t) - R 

8 g' g 


15.80) 


This is a typical equation of motion for a single-degree dynamic system. 

But the design method described in paragraph 5-12 is not applicable to the 

determination of R for one reason. In paragraph 5-12^ the resistance 

g 

function depends on the absolute displacement of the mass supported by the 

columns. But in equation ( 5 . 8 O), although x is the absolute dlsplace- 

g 

ment, the resistance function R depends on the relative displacement 

g 

between the floors g and g - 1 . Hence, in the preliminary design, it 
is necessary to consider the relative displacement between floors as well 
as the absolute displacement of each floor. 


The relative displacements, X^, Xg. 
following equations: 


.X are defined by the 
n 


\ - ^2 


X = X - X - 
g g g - 1 


(5.81) 


X = X - X , 
n n n - 1 


Substituting equation ( 5 . 81 ) into equations (5-79h), and (5.79c)» 

the exact equations of motion in terms of relative displacements are ob¬ 


tained. By neglecting R 


g + 1 


in the preliminary design of columns for 
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5-2Ib 


story s , ’the equation of motion for mass, , in terms of relative dis¬ 
placement, X , is given by: 

O 


"g -1" ^ m 




(5-8e) 


where 


fgx(t) - fg(t) - fg . i(t) + ^ Rg . 1 

end defined as the effective load for story g . The relative 

displaceiaent, X , depends on this effective load which is a completely 
known function since f ^(t) and f (t) are the given loads and the 


resistance function, R , , is known from the preliminary design of the 

columns for story g - 1 . Once f ^C'^) known, the relative displace- 

gx. 

ment, X , and the resistance, E , at any time can be evaluated. 

-j- Cl o 

*" A typical shape of the effec¬ 

tive load, f ^(t) , is shown in 
^ gx 

figure 5*37a» The corresponding rela- 

•I 

5 tive displacement and resistance func¬ 

tion are given in figure 5 * 37 b and 
figure 5-370 respectively. For pre- 
limlnsry design purposes, the actual 
resistance fimction shown in figure 
® 5-370 is replaced by an approximate 

resistance function shown in figure 

5 5.37d. The method of evaluating the 

ti 

parameters of the approximate resist- 

“ ance function, namely R , t*, t , 

g g ge' 

and R is discussed in the follow- 
gm 

i, ing pareigraph. 

• S 

S2 b. Linear Effective Load. In 

•« . . 

• A « 

i£ the preliminary design of two, three, 

- or four story frame buildings, the 

resistance functions for a linear effec* 

tive load, no reverse yielding initial portion of the effective load. 



resistance functions for a linear effec¬ 
tive load, no reverse yielding 
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for any story above the first can be approximated, with sufficient accuracy, 
by a linear effective load as shown in figure 5.37a. The initial portion 
of the actual^effective load is generally always negative because the initif 
motion of the lower floor masses is greater than the upper floor masses. 
•Ihis situation is reversed as the deflections increase. The general expres¬ 
sion of the linear effective load is: 




B 


gX 


(‘-i) 


( 5 . 83 ) 


Note that is a negative quantity for reasons mentioned above. Approx: 

, are 


mate formulas for the evaluation of the parameters, B „ and T ^ 

gx gx 

given in paragraph 5-23- 

Substituting eqviatlon (5-83) into equation (5.82), the following 
equation is obtained: 

d^x 

B_^(l - ~ \ - f 1 + S U-X (5.8U) 

O O 


m 

g 


dt^ 


gx 




R 


The resistance function, 
equation (5.84). 

Prom equation (5-84), the time. 


g , in equation (5-82) is replaced by kgX^ in 

is zero, can be eval- 


t* , when 
g 


g 


uated. 
0 to 


The average value of , that is R , in the time interval from 

g 

dX /dt at t* also can be evaluated 
g g 


g 

t* , and the relative velocity 
g 


These quantities, evaluated on the assumption that the columns in story g 

remain elastic when X is negative, are plotted in non-dimensional form 

g 

in figure 5*38. 

In actml design, a trial column section is selected from which the 
spring constant. 


the ultimate resistance. 


k , the maximum elastic relative deflection, X , and 
g a® 

R_ , can be computed. The value of T_v/T’ is 


gm 


‘gX' '■gn 


determined where 


T* 

gn 


is the period of a fictitious dynamic system ex¬ 


pressed by equation (5*84), and given by equation (5.85). 


T‘ 

gn 


®*«/f— 

y[i. 


m 


3- 


("A -1)]' 


(5.85) 


g 


Using 




in figure 5.38, the non-dimensional quantities; 
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c. Equivalent System for Use with Design Cliarts. The design charts 
of paragraph 5-10 were constructed for the design of single-degree-of- 
freedom dynamic systems. These charts can he used to simplify the design 
of multi-story frame buildings by introducing an equivalent single-degree- 
of-freedom system for the g-th story. 

In the previous paragraph the approximate resistance, R , at any 

O 

time for the g-th story was developed. Suppose the external load, f (t), 
is a triangular load given by: 



(5.87) 


The load and resistance functions acting on m are shown in figure 5*39®’* 
The resistance function does not resemble the resistance function for a 
single-degree-of-freedom system with zero Initial displacement and velocity. 
Therefore it must he replaced by the resistance function, [^®g(^^eq ^ 

equivalent single-degree dynamic system shown in figure 5*39b. This is ac¬ 
complished by revising the external load in a manner consistent with the 
change in the resistance function. For the equivalent system, the maximum 



Ti me Time 


(a) Actual System 


(b) Equivalent Single-Degree 
System 


Figure 5.39. Equivalent single-degree dynamic system for the g-th story of a frame building 
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5-21c 


resistance is chosen to he equal to R_ and the time when the maximum 

gm 

elastic deflection occurs is equal to t^^ . The relationship between 


[V*>] 


and R_(t) is given by: 
eq B 




= R (t) - R* + R" 
eq ^ 


( 5 . 88 ) 


or 


K^(t) - 


eq 


+ R» - R" 


where R* and R" are the curves of the approximate and equivalent re¬ 
sistance functions respectively for times less than t . R" is assumed 

ge 

to be a parabolic curve. 

Substituting equation (5.88) into equation (5.8o), the following 
equation is obtained: 




m 


■g 


dt 


(5.89) 


Equation ( 5 . 89 ) indicates that the movement of the floor mass, m , can be 

O 

represented by the movement of the mass of an equivalent single-degree dy¬ 
namic system. The external load on the equivalent system is: 

[^fg(t) - R« + R"j 

instead of f (t) . The work done on the equivalent system is; 

s 

t 

\ ' /[V‘> ■ "'■] i “X ^ ^ 

The maximum allowable strain energy of the equivalent system is: 


/[i>,(t)] ^ ax B" S dt . - B_ 5? dt 




ge 


"gm ^ 
gm dt 
ge 


R" ^ 'i't + R fx - X 
dt gm I gm ge 


(^gm “ ^ge) 


( 5 . 9 a 


where • x is the absolute displacement of m at time, t 

ge g gc 
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The first term in equation (5-90) is identical with the first term in 
equation ( 5 . 91 ). Hence, in comparing the maximum work done and allowable 
energy absorption of the equivalent system, only the second terms in equa¬ 
tions ( 5 * 90 ) and (5- 91 ) need be considered. The method described in para¬ 
graph 5* 12a can be used for the evaluation of the maximum work done. The 
absolute maximum work done by the load ™ist be first eval¬ 

uated. Next, figure 5*27 can be used to obtained the work done ratio, and 
the actml work done and allowable energy absorption determine the suit¬ 
ability of the selected trial section. 

It is seen from the above discussion that, by introducing an equiva¬ 
lent single-degree dynamic system, the design charts given in paragraph 
5-10 for single-degree dynamic systems can also be used for the design of 
multi-story frame buildings. The methods of evaluating the fictitious maxi¬ 
mum work done, work done ratio, and allowable energy absorption for the 
equivalent system of the g-th story are described in the following three 
paragraphs. 

Ficitious Maximum Work Done, W . !]3ie fictitious maximum work done 

_ gp 


by the load is given by; 


U -& 

gp Sfflg 




fg(t) - R*(t)j 




( 5 . 92 ) 


For a given value of T „/t* , the values of R„ and t * cam be obtained 

gA' gn ' 6 g 

from figure 5 . 38 . Within the range of 

0.15 < T „/T‘ <1-0 

gx' gn 


the value of R-t* varies from 

O 


0.6 


1 + 


S’ g - 1) 


to 0.85 




^(V^g -1) 

For the purpose of preliminary design, the following empirical equation may 
be used: 

0.75 ■ (5.93) 


g g 


^ ^ (V-g -1) 
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Substituting the value of R t’ obtained from ecLiiation (5-93) into eq.ua- 

s s 

tion ( 5 - 92 ), the fictitious maximum work done, W , can be determined. 

SP 


Work Done Ratio, C 

^ r 


The work done ratio, = Was shown in 


figure 5*27 depends on and . For the eq.uivalent system of the 
g-th story, C „ can be approximated from the relation: 


‘^gR “ 


I R t* 

1 +_e_s 


( 5 . 9 »^) 


but the ratio C _ = T /T^ 
gT g' gn 


has to be evaluated before 


can be de¬ 


termined from figure 5-27- 

A single-degree dynamic system involves three non-dimensional param¬ 
eters. For the equivalent system of the g-th story, the parameters are 


chosen to be C t /t , and C T , where 
gP ge' g' g ' 


is the time when the 


maximum elastic deflection. 



Figure 5.40. ^ge^'^g curves for elasto-plastic systems, 
triangular load 


X , occurs. The relation- 
ge * 

ship between these parameters 
is given in figure 5**+0 for 
the case of a triangular load. 
This figure can be used for 
the evaluation of any param¬ 
eter when the other two j)a- 
rameters are known. For the 
equivalent system of a g-th 
story, is given by 

equation (5*9^)> t is 
given by equation ( 5 . 86 ), and 

T„ is the given load dura- 
g 

tion. Hence, C „ and 
gP 

known. The C 
ge' g gT 

ratio can be obtained from 
figure 5-^0. The value thus 
obtained is defined as the 
equivalent value and denoted 
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[^gl] 


Using C 


(v) 


, the work done ratio, C „ 
eq ^ gW 


then be determined from figure 5-27 from which the maximum work done, W 


can be determined. Prom figure 5-29, the time ratio, is obtained, 

and the time, t , when X occ-urs can be determined. 


Maximum Allowable Energy Absorption. The maximum allowable energy- 
absorption as represented by the second term in equation ( 5 . 91 ) is: 


E_ = R 


gm 1 gm 


(5.95) 


The value, x , obtained from equation ( 5 . 80 ) by integration, is 


given approximately by: 


^ge “ 


(5.96) 


Comparison of Maximum Work Done and Maximum Allowable Energy Absorp¬ 


tion. The value-of W is given by equation (5 •92). The work done ratio, 

GJr 

C , can be evaluated by the method described in the preceding paragraph, 

gw 

The maximum allowable strain energy, E , is determined from equation 

o 

(5»95)* Comparison of the maximum work done and maximum allowable energy 
absorption determines the suitability of the selected trial section. 

Prom the above discussion, it is seen that, by introducing the equiva¬ 
lent system for the g-th story, and by using an equivalent value, > 

the design charts, figures 5 >27, 5 * 28 , and 5 . 29 , which are prepared for 
single-degree dynamic systems, can be used for the preliminary design of 
the columns for a multi-story frame building. The detail procedixre for de¬ 
sign is given in paragraph 5 - 23 . 

5-22 ADJUSTMENT OF PRELIMINARY CQUJMN RESISTANCE. In the design procedure 
outlined in paragraph 5 - 19 ^ the preliminary design for each story is made 
assuming all hi^er stories are absent. That is while making a preliminary 
design of the g-th story, the presence of stories g + 1, g + 2, etc. is 
Ignored. Actually, there is a force, + q ^ acting on mass, m^ , from 
the columns of story g + 1 (see figure 5 . 38 ). The effect of this force 
may be taken into account by an adjustment of the preliminary resistance of 
the colxmms for the g-th story. 
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5-23 


The equation of motion for , when hoth R 


sidered, is given by equation ( 5 - 97 ) 

"g ^\ ®g +1' \ 1 1] 


g + 1 


are con- 


( 5 . 97 ) 


According to equation (5•97)# the resistance function, R 


can be 


considered as am additional load on m 


g + 1 

The approximate shapes of R 


and Rg as functions of time are shown in figure 5 >41. When R^ 


neglected in the preliminary design of R , the maximum displacement, x- 

g g 


occurs at 


^(g + l)m 


The problem now is 


to determine by what amount R 


(g +l)e 


should be increased so that when the 
additional load, R 1 is acting 


—- - 

Rg +1 

1 


Rgm t 

t 


^gm 

Rg 



g + 1 ^ 

, the maximum allowable dis- 


.placement of m is still equal to 

O 


The additional load, R 


g + 1 ' 


changes both the velocity and dis- 


Figure 5.41. Approximate resistance functions 
Rg+ j(t)and R^(t) 


placement of m at t . Sup- 
pose Rg^ is changed by an amoxmt 
equal to AR as shown in figure 


5.1+1; the magnitude of AR should 


be of such value as to cancel the ef¬ 


fect of R„ ^ , 
g + 1 

equation ( 5 - 98 )- 


The modified value of R is given approximately by 


Modified R 


R + AR = R 
gm gm gm 


r-^ 


JL±1 

I 


( 5 . 98 ) 


where 


Rg^ = value obtained in the preliminary design 

I=:R It -t\-Rt* 
g gm \ gm ge) g g 


g + 1 > g 


= R 

(t - 

t 1 

gm 

\ gm 

gey 

= R, 

(e 

+ l)m 

ft 

L gm 


'(g + l)e 


R ,t» 

g + 1 g + 1 


5-23 COMPLETE DESIGN PROCEDURE FOR MULTI-STORY FRAME BUILDINGS. A complett 
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design procedinre for the design of multi-story frame buildings is given in 
this paragraph. 

The procedure is set up on the ass-umptions that (l) the external 
loads can be approximated by triangular load curves, and (2) the actual ef¬ 
fective loads can be approximated by linear effective loads. The procedure 
is applicable only to the design of two, three, or four story buildings. 

The extension of this procedure to the design of higher buildings is dis¬ 
cussed in paragraph 5-24. Numerical examples illustrating the application 
of this procedure to the* design of two and three story buildings are given 
in EM 1110-345-418. 

a. Preliminary Design of the First Story Col-umns. In the prelimi¬ 
nary design, the first floor is considered as a single-degree dynamic 
system with mass, , subjected to external load, f^(t) , and all the 
other floors are assumed to be absent. The energy method of design given 
in paragraph 5-12a is used. The maximum elastic deflection, , the 

plastic resistance, , and the spring constant, , of the selected 

column section are determined. The time, t^ , is determined from figure 
5.29 the time, t^^ , is determined from figure 5-40. 

h. Preliminary Design of Col\mins Above the First Story. The proce¬ 
dure given in this section is applicable to the design of the columns of 
any story above the first story. Minor differences in the detail computa¬ 
tion for different stories are noted whenever it is necessary. 

Step 1. Evaluate the parameters, and T^ , of the linear ef¬ 

fective load according to the following approximate formulas: 

For the second story: 


®2X “ ®1 ” 

^^2 " \e”l®2X 

For any story above the second story: 

V - (V"g - 1) - 1 ■ - 1 


(5-99) 

( 5 . 100 ) 


( 5 . 101 ) 
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5-23b 



( 5 - 102 ) 


where B,, 3„, B are the peak values of the external loads, 

1' 2-' g 'I' ■" 

f2(t), and and, T^, T^, are the d\jrations of the external 

loads. In determining ®g^ Tg, and T^, the given external 

loads should be approximated by triangular load curves. Depending on which 
story is under consideration, the approximation should be close in the fol¬ 
lowing time interval: 

for the second story, 0 < t < 2t^^ , 

for the third story, 0-< t < St^^ , and 

for the fourth story, 0 < t < ^t^^ . 

In the design of buildings that are higher than two stories the value of 
B^ may be zero or negative in which case the procedxrre of paragraph 5-24 
should be used. 

Step 2. Compute R t* according to equation (5-93)- 

S 6 



step 3* Compute the fictitious maximum work done, , according 


to equation (5.92). 


W =7“ 

gP 2m„ 


-S + R t* 


Before doing this, the given load, f (t) , is approximated by a re- 

O 


vised triangular load curve. Depending on which-story is under considera¬ 
tion, the appi'oximation should be close in the following time interval: 
for the second story, 0 < t < , 

for the third story, 0 < t < 1.3t» , 


for the fourth story, 0 < t < 1.3t2jjj f and generally 


for the g-th story, 0 < t < l-3t^g _ 


The values of B, 


Tg obtained for the revised approximate tri¬ 


angular load will in general differ from those obtained in step 1 because 
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they simulate the given load over a different interval of time. 

Step i4-. Find , As a first approximation, the following equa¬ 

tions may he used. 


for any intermediate story, R = W /x 

m gp' gm 

for the top story, R = I.3W /x 

Sm gp' gm 


(5.103) 

(5.104) 


where x^ is the maximum allowable absolute displacement for the g-th 
story and is given by equations similar to equation (5.78), 

^Im ~ ^Im 

^2m ^Im ^2m ' ^le 
^3m “ ^Im ^2m ^3m " ^le " ^2e 


Step 5" Estimate the maximum elastic deflection, X . If the 

ge 

heights of the first and g-th stories are equal, X may be assumed equal 

gG 

to Xj^g . Determine spring constant, k , from the relation. 


K - 

g gm' ge 

Step 6. Determine according to equation (5.85)» 


(5.105) 


T' * 2« 
gn 4 


- -B--— 

1 + (m /m T 1 k 

L \ g' g - V J e 


step 7. Determine T^/T^ where is given by equation (5.100), 


or equation (5.102). 

Step 8. Using T ^/T' in figure 5.38^ find 


R / m \ 

^ \ g - 1/ 


and-^(l +- 


ti A k T 

-1/ V 


Then compute t' , R , and dX /dt . 

00 O 

Step 9. Compute t according to equation (5.86), 

g6 


^ge “ ■’'g SdXg/db 


Step 10. Compute t_ /T^ and C ^ from equation (5-9^), 

g0 g gA 
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5-23C 


■ 


R /B 
gm^ g 


1 + 


R t’ 
g g 

g g 


Step 11. 

5.40. 

Using 

t /T , and C - . determine 
ge' g ^ gR ' 

(V)eq 

from figure 

Step 12. 

Using 

C _ and 1 C „ 1 , determine 

gR \ gT/eq ’ 

C and 

gW 


figures 5-27 and 5-29 

respectively. Then compute 

and t from the 

gm 


relations. 


W = C ,,W 
gm gW gp 


t = 
gm 


(5.106) 

( 5 .107) 


The value of t 


gm 


= /t /t \t 
^ gm^ gj g 

is used to check whether the triangular load approxima¬ 


tion in step 3 is acceptable. That is^ if is much different from 

1.3t(g l)m assumed in step 3> a revision of the triangular load ap¬ 
proximation may be necessary. 

according to equation ( 5 * 96 )> 

B 

i: 

g 


Step 13 . Compute x 


ge 


ivi: 


X = 

ge 2m„ 


B t 
g 

3T 


^+R 


Step l4. Determine the maxim-um allowable strain energy, 
cording to equation (5-95)• 


Eg , ac- 


E = R 


gm 




- X 1 

gm ge) 

Step 15 * Compaj'ison of E and W determines the suitability of 

g g^ 

R . If the agreement between E and W is not within lOjto. the pro- 
gra g gm / 7 X' 

cedure from step 3 to step 1^+ is repeated by using a new value of 
given by^ 


W 


R = 
gm 


gm 


2[x - X \ 

\ gm gej 


(5.103) 


and W agree within lOfo, the value of R obtained 
gm ° " gm 


However, if . E 

from equation (5-108) is approximately the maximum resistance required, 
c. Adjustment of Preliminary Column Resistance R 


gm 


the g + 1 story columns on the movement of floor mass, 


m 


g 


The effect of 
, is taken into 


93 



5-23<i 


EM IIIO-345-IH5 
15 Mar 57 


account by modifying the value of . The modified value of R is 

igm gm 


given by equation (5*98)> as 


where 


Modified R = R 

gm gm 


{ 


+ I 

g g + 1 , g 

^g 


) 


R g^ = value obtained in the preliminary design 


It - t \ 

g gm ^ gm gel 


R t' 
g g 


g + 


1 , g ^(g + l)m [^gm ■ "^(g + l)e 


- R -,t' , 

g + 1 g + 1 


Starting from the top story, the plastic resistances of different stories 

should be modified one after the other. In the modification of R by 

gm 

equation (5*98 )j "the value of 2.)m modified value 

rather than the value obtained in the preliminary design. In the computa¬ 
tions of I^ for the bottom story, the value of R^ is zero and the value 
of tg^ is equal to 2t^^/3 . 

d. Selection of Final Coltimn Sections. The maximum horizontal re¬ 
sistance of the selected columns for the g-th story, including the effect 
of vertical loads acting on the columns as outlined below, should be close 

to the modified value of R obtained in the previous paragraph. In de- 

gm 

termining the effect of vertical load, its average value in the time 
interval, 

t < t < t 
ge gm 

should be used. The expression of the ultimate resistance of the selected 
columns is given by: 


- a 


(5.109) 


where 


= plastic resisting moment of the selected section after the ef- 
^ feet of the vertical load is taken into account. 

ZP =■ the time average of the total vertical load in the time interval, 

^ t < t < t 

ge gm 


n = number of columns in the g-th story. 
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h = height of the g-th story. 

s 

After the column section has been selected, the spring constant, k 

6 

and the limiting elastic deflection, X , are determined. In the evalua- 

g6 

tion of these parameters, the effect of girder flexibility may be approxi-, 

mately accounted fpr (see EM 1110-345-417 and EM 1110-345-418). 

e. Nvunerical Analysis of the Entire Building. A numerical analysis 

for the entire buildix^ should be carried out to determine the suitability 

of the selected columns. If the result of the numerical ainalysis Indicates 

that the maximum absolute displacement, x > is greatly different from 

the design specification, all the columns in eind below the g-th story 

should be revised. The revised value of R can be determined by a trial 

gm 

and error method or from the relation. 


t 





(x - x' \ 

1 

( x' 

- X \ 

1 gm gej 

g\ gm/ 

1 gm 

gm| 


X - X’ 

gm ge 


(5.110) 


where 


R* 

gm 


and R^g ^ are the average values of the maximum resistance 

used in the numerical analysis over the time interval, t < t < t ; and 

^ ge gm . 

x^ and Xg^ are the displacements obtained from the numerical analysis. 

R/ ^ V is the revised maximum resistance of story g 4- 1 : x 

(S IM , j 


prescribed maximum allowable absolute displacement; and f 


g 


l^Sm) 


gm 


magnitude of the external load, f (t) , at the time when x* 

g gm 


is the 
is the 
occurs. 

New column sections should be selected according to the revised 
values of the plastic resistances. The spring constants and the maximum 
elastic deflection should then be determined. The numerical analysis 
should be repeated using the revised column sections. 


5-24 PROCEDURE FOR BUILDINGS OF MORE THAN TWO STORIES. The design pro¬ 
cedure given in paragraph 5-23 is for the case of a linear effective load 
which is applicable only to buildings two stories or less in height. In 
the design of higher buildings (higher than two stories) the value of B 
computed from eq.uation (5-99) '^7 t>e zero or negative. The effective load, 
f2 jj(t) , and the relative displacement, Xg , are then positive immediately 
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after time, t = 0 . Thus when computed from equation ( 5 . 99 ) is zero 

or negative, the time, t^^ , can be taken to be 5t^g/3 . Similarly, if 
computed from equation (5.101) is zero or negative, the value of t 


can be taken to be 


."^20 ^ 


For buildings of more than two stories. 


if B „ is close to zero or negative. 


\j\j 4tiv-.jL w Vw/JL o J. V j.\^ should be considered eqioal to 

zero and the procedxires from step 5 to step 9 in paragraph 5 - 23 b should be 

omitted in the preliminary design of the second or third story colxnnns.,, 

5-25 ACCURACY OF FLOOR-BY-FLOOR DESIGN METHOD. A number of assumptions 

have been used in this design method- First is the replacement of actual 

effective load by a linear load in which the pareimeters, B ,, and T „ 

& gX 

are determined by approximate formulas- Second, approximate formulas are 
used in the determination of t„ and C „ . Third, in the evaluation of 
strain energies and work done, the portion of R"(t) shown in figure 5.39b 
is neglected. And fourth, in the determination of the plastic resistance, 

, of the chosen colxunn section from equation ( 5 .IO 9 ), the time average 
of the total vertical load on the columns is used. 


There is error involved in each assumption. The accumulative error 
for am actual design problem may be quite large. However, based on a 
number of examples which have been worked covering a wide range of build¬ 
ings under various loadings, the maximum resistance of the columns for any 
story obtained by this floor-by-floor iteration method has been found to 
be generally within 20 ^ of the desired value, for cases in which the pre¬ 
scribed maximum allowable relative displacement was at least three times 
larger than the maximum relative elastic deflection. 

For multi-story frame buildings, a small change in the plastic re¬ 
sistance generally results in a large change in the relative displacement. 
Although the result of the numerical analysis may indicate a large error 
in the relative displacement, this usvially represents small error in the 
plastic resistance. Since the plastic resistance cannot be evaluated 
exactly, accuracy greater than 20^ in the design of coltmms on the basis of 
relative displacement is not warranted. An error of 205t in the relative 
displacement may correspond to an error of only 5^ in the plastic 


resistance. 


101 


EM 1110 - 345 - 14-15 
15 Max 57 


5-26 


DESIGN OF MQLTI-STORY BUILDINGS FOR ELASTIC DEFORMATION 

5-26 INTRODUCTION. In the design of a multi-story building for elastic 
deformation the maximum allowable lateral deflection of the columns, under 
dynamic load, should not exceed that which causes yielding. In the equiva¬ 
lent system for frame buildings as shown in figure 5*33» the moment in the 

columns of the g-th story is proportional to the relative displacement X . 

6 

Direct stresses are developed in the col\amns due to over-all bending but 
these may be ignored in buildings of three or four stories. In frame 
buildings, where the direct stress due to over-all bending is appreciable, 
the preliminary design could still be made without including the effect of 
the direct stress. This would be considered, however, in selecting the 
column sections. Hence, the elastic design of frame buildings for elastic 
deformation can be based on the requirement that the maximum allowable 
relative displacement for any story must be smaller than the relative 
elastic deflection at the yield point stress. 

However, in the equivalent system of a shear wall building as shown 
in figure 5 - 32 , the stresses in the shear wall at any story axe a function 
of the over-all bending and shear in the wall and hence depend on the rela¬ 
tive displacements of all stories rather than just the relative displacement 
of adjacent stories. The use of absolute displacements is desirable in the 
design of shear walls because it is the simplest method to express the rela¬ 
tive displacement between all floors in a shear wall building, whether the 
floors are adjacent or non-adjacent. 

The characteristics of a multi-story building designed for plastic 
deformation discussed in paragraph 5-20 are also generally true when de¬ 
signed for elastic deformation. These characteristics are: (l) the rela¬ 
tive displacements in the upper stories are negative shortly after the 
initiation of the load; and (2) the elastic rebound of floor masses must 
be considered in the evaluation of the relative displacements. In design 
for plaatic deformation, the resistance, R , is assumed to be a constant 

O 

after the elastic limit is reached. The design of the g-th story colxomns 
is reduced to the determination of the shape R as shown in figure 5'37d. 

s 

Moreover, the elastic rebound of floor masses is considerably smaller than 
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the allowable relative displacement. Hence, the use of approximate formula 
in estimating the magnitude of elastic rebo\ind is permissible. But in an 
elastic system, the resistance, R , always varies with time. The elastic 

O 

rebound is about of the same order of magnitude as the relative displace¬ 
ment, and approximate evaluation of the elastic rebound is not permissible. 
Because of these differences, the elastic design of multi-story buildings 
is more involved than the design for plastic deformation. In most cases, 
a trial and error design procedure may have to be used. The elastic de¬ 
sign of multi-story frame and shear wall buildings is briefly discussed 
in paragraphs 5-27 and 5-28 respectively. The design of shear* wail build¬ 
ings is more fully discussed in EM 1110-3^5-^19* 

5-27 DESIGN OF MULTI-STORY FRAME BUILDINGS FOR EIASTIC DEFORMATION. One 
possible approach -to the elastic design of multi-story buildings is the 
preparation of non-dimensional graphs. The feasibility of this method is 
illustrated by the characteristics of a two-story frame building. The 
equivalent dynamic system of the 
building is shown in figure 
5 .U 2 . The external loads, f^^Ct) 
euad f 2 ('*-) assumed to be 

triangular with load durations 
T^ and Tg respectively. The 

relative displacements, B| (l-t/Tj ) = f|{t) 


Bad-t/Tg) 


and X_ when the peak external 
loads are applied statically, 
are given by 



///////77y////77 7 7 


( = —- 

'•as 2k. 


■5.111) 


Figure 5,42. Equivalent dynamic system for 
tuio-story building 


"Is 


+ Bg 

2e; 


( 5 . 112 ) 


X, 


and 


X denote the maximum relative displacements under dy- 
Im 


Suppose _ 

namic load; the dynamic load factors for the second and first stories are 

defined by the following relations; 


(D.L.F. = Xg^X 


2s 


(5.113) 
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(D.L.P, 


(5.1li^) 


These dynamic load factors depend on the non-dimensional parsuneters ^ 2 /^ * 
Tg/^ip \/^in * ^]ji equal to 2«‘\/m^/k^ . The gen¬ 

eral shape of the dynamic load factor curves is shown in figure 5*^3 for 

specific values of Tg/T^^, 

and Bg/B^ . Before these curves 

11 used for practical design, a 

11 ^ series of similar graphs must be pre- 

i “ '- pared for different combinations of 

Si -- , / / 

O V—kj/k, incr»fl»l(>g ^2'^1 COVeT 

'- the probable range of practical 

T,/T., 

buildings. For practical cases, it 
I g. . is believed the range of m^/m^, 

^ within narrow limits. It is even 

f*" XyT / X^ki/k, inerfoting possiblc that they might be ass\jmed 

___X ^^onstant for certain types of eon- 

struction. Suppose foinr different 
values eire used for each non- 


k^/ki Incrtatino 


Figure 5.43. Dynamic load factors for two-story . 

L ij- t ! j ■ / - .... dimensional parameter, then the 

building for elastic design, m 2 /mj - constant, * 

T 2 /TJ- constant, B 2 /Bi = constant total number of graphs needed is 6h- 

By similar considerations for a 

three-stdry building, the total number of graphs required is 2h01. Since 
such a Istrge number of graphs are needed, it appears Impractical in general 
to prepare non-dimensional graphs for use in the design of multi-story dy¬ 
namic systems. 

For practical design, a trial and error procedure may have to be 
used. A set of columns for the entire building is selected. The suit¬ 
ability of the selected columns is determined by the result of a numerical 
analysis. In the preliminary selection of columns for the g-th story, an 
equivalent single-degree dynamic system, shown in figure 5-^^> may be used. 
In the equivalent system, all the floor masses and external loads above 
the g-th story are assmed to be concentrated at the g-th floor. The 
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Mg = (Dg +mg+| + mg^.2+...m„ 

FgCt) 

= fg(t)+fg + |(t) 

+-fn(t) 


Figure 5.44. Equivalent single-degree system for the g-th story 

design method on the basis of dynamic load factor given in paragraph 5-12c 
is used for the determination of the columns for the g-th story. The 
system shown in figure 5.44 is used for the preliminary design of columns 
in the first, second, up to (n - l)th story. For the top story, the ex¬ 
ternal load is assumed to be 2f^(t) instead of fj^Ct) . 

Because of higher modes of oscillation in multi-story buildings, a 
large number of maximum and minimum values are presented when the relative 
displacements are plotted versus time. Hence, in the n-umerical analysis, 
the integration process should be carried out far enough to insvire that 
the absolute maximum relative displacements for all stories are obtained. 
5-28 DESIGN OF MULTI-STORY SHEAR WALL BUILDINGS FOR EIASTIC DEFORMATION. 

In the elastic design of shear wall buildings, the controlling factor may 
be shear and/or bending stresses. For the evaluation of both stresses, the 
building is replaced by a multi-degree dynamic system discussed in para¬ 
graph 5-l6a. The method of evaluating the parameters of the dynamic system, 
the detailed preliminary design procedure, a discussion of design tech¬ 
niques for both elastic and plastic action, and numerical examples of 
shear wall building design are given in EM 1110-345-419. 
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EFFECT OF VARIATION OF PARAMETERS ON THE RESPONSE 
OF SINGLE-DEGREE DYNAMIC SYSTEMS 

A-Ol INTOODUCTION. In paragraph 5-10, the response of single-degree dy¬ 
namic systems subjected to simplified loads is presented in non-dimensional 
graphs. For a given problem, when the peak load B , the load duration T , 

the mass m , the spring constant k and the plastic resistance R are 

m 

all known, the maximum work done, , the maximum displacement, x . and 

in ' m * 

the time, t^ , when occurs can be determined from these non- 

dimensional graphs. For the convenience of discussion, the first five var¬ 
iables (B, T, m, k, and R^) are considered as independent' variables in 
this appendix. 

For design applications, in addition to these non-dimensional graphs, 
it would be helpful to know how the response of a dynamic system varies 
with the shape, the peak value, and the duration of the load pulse, and how 
the response is affected by the mass, the spring constant, and the plastic 
resistance. In other words, it is of interest to know the effects of small 
Increments of the independent variables on the response. Two practical de¬ 
sign problems, in which this is of importance, are the following; 

The first problem is concerned with uncertainties. Since all the in¬ 
dependent variables are subject to uncertainties in their determination, it 
would be of interest to know what is the effect of these \incertainties on 
the response and what is the necessary accuracy in the determination of the 
variables such that the response of a dynamics system is within a given 
tolerable limit. 

The second problem occurs when it is intended to answer the following 
two questions which often arise at the completion of a trial design compu¬ 
tation. Oliese questions are; (l) What is the permissible difference be¬ 
tween the specified maximum displacement and the corresponding computed 
value of a given design; and (2) if a trial design is judged to be .unsatis¬ 
factory, by what amounts should the variables be changed in the next trial 
design. For numerical examples, suppose for a given system a 10^ change 
in the mass results in change in the maximum displacement, then an 


107 



EM 1110-345-^15 
15 Mar 57 


A-01 


accurate determination of the mass would not be necessary. On the other 
hand; if a 10^ change in the plastic resistance results in a 50 ^ change 
in , a close agreement between the specified maximum displacement and 
the corresponding computed value would not be warranted since the deter¬ 
mination of the plastic resistance may be subject to an uncertainty greater 
than 10^. 

The effects of small increments of the independent variables on the 
response are discussed in this appendix. These effects are expressed in 
terms of percentage-increment ratios which are non-dimensional quantities 
and are defined in paragraph A-02. In paragraph A- 03 , the technique of 
non-dimensional analysis is applied in deriving certain functional rela¬ 
tionships among the percentage-increment ratios. Because of these rela¬ 
tionships, the effect of small increments of all of the variables (B, T, 
m, k, and can be condensed into two sets of curves. The percentage- 

increment ratio for linearly elastic, completely plastic, and elasto- 
plastic systems subjected to both rectangular and triangular loads are dis¬ 
cussed in paragraphs A-04, A-05, and A -06 respectively. The effect of load, 
shape on the response is discussed in paragraph A-OT- The expressions of 
various percentage-increment ratios for linearly elastic, completely 
plastic, and elasto-plastic systems subjected to rectangular and triangular, 
loads are summarized in paragraph A- 08 . 

It is to be noted that all of the material in this appendix is de¬ 
rived from the corresponding material in EM 1110-345-415 through a process 

lent ratios are related to the 

and the graphs in this ap¬ 
pendix are related to the slopes of the corresponding graphs in paragraph 
5-10. Although these two sets of graphs are derivable from each other, 
they are intended for different applications. Generally speaking, for 
routine design computation, it would be easier and the result more accu¬ 
rate when the graphs in paragraph 5-10 are used, while for those special 
applications discussed earlier in this paragraph, it would be better to 
use the graphs in this appendix. Beside these special applications, the 
graphs of percentage-increment ratios also illustrate many important prop¬ 
erties of structures under dynamic loads. For this reason, although 


of differentiation. The percentage^rincre 
partial derivatives, (for example, , 
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quantitative data are presented in figures A.l ttirough A. 8, the discussion 
in paragraphs A-04, A-05, and A -06 is more or less qualitative in nature. 
This quantitative discussion should then serve as a rule of thumb for a 
designer to make proper Judgement on the many practical problems which are 
not necessarily related to any routine design computation. 


A-02 SENSITTVrrr AND EERCENTAGE-INCKEIMEHT RATIO. 


Sensitivity. The 


maximum displacement, , generally depends on all the variables, B, 

T, m, k, and R^ . Suppose for the time being only the load dura.tion, 
T , is varying. Then is a function of T . If the load duration is 
changed from T to T + AT , the corresponding maximum displacement is 


changed from x to x + Ax 
m m m 

related by equation (A.l): 


For small increments. Ax and Zd? are 

m 


Ax = 

HI 


(A.l) 


Ax Bx_ 
q_m m 

xT " AT ~ “5T 


(A.2) 


where 

AT « the increment of the load dviration. 

Zflc = the corresponding increment of the maximum displacement. 

3x 

—^ = a partial derivative. - 

The quantity, S^^ , being a partial derivative, is defined as the Sensi¬ 
tivity of the maximum displacement with respect to the load duration. Thus 
the sensitivity S^^ is equal to the ratio of two increments Ax^ and AT 
and is a dimensional quantity: For example, when S^ is equal to 0.01 
inch/sec, this indicates that if the load diiration is increased by 0.1 
seconds, the maximum displacement would be increased by 0.001 inches pro¬ 
vided the increments are small in comparison with their original values. 

Althou^ the dimensional quantity, sensitivity, is very useful in 
numerical computation, it does not express clearly the physical significance 
of "sensitivity." From the numerical magnitude of the sensitivity, it is- 

109 


EM 1110-3*^5-^+15 
15 Mar 57 


A-02b 


difficult to perceive whether the maximum displacement is "very sensitive" 

or "not sensitive" to any variation of the load duration. This is because 

the original values x and T to which Zflc and AT are to be added 

m m 

are not included in the definition of the sensitivity. Hence in this ap¬ 
pendix, instead of using increments and sensitivities, the non-dimensional 
quantities, percentage-increments and percentage-increment ratios are used. 
These non-dimensional quantities are defined in paragraph A-02b. 

b. Percentage-Increment Ratios. If the left hand side of equation 

- Xjjj 

(A.l) is divided by x^^^ and the right hand side is divided by T , the 
result is given by equation (A.3). 



= I 


xT T 


(A.3) 


where 



percentage-increment in x^ 


T 


percentage-increment in 


T 


and 

T - — 
xT ~ x^ St 
m 

= The ratio of the percentage change in x to the percentage 
change in T while all the remaining independent variables are 
held constant. 

The quantity I^^ , with two subscripts, is defined as the percentage- 
increment ratio. 

In the discussion so far, x and T are considered respectively 

in 

as the dependent and independent variables. By a similar procedure, if the 

variables are the maximum work done and the peak load B , the percentage- 

increment ratio I^^ would give the ratio of the percentage-increment of 

W . to the corresponding percentage-increment of the peak load B . 

M 

There are a number of percentage-increment ratios for a single-degree 
dynamic system subjected to simplified loads. These ratios are described 
in paragraph A-02c. 

no 
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Percentage- Increment Ratios for Single-Degree Dynamic 
®- single-degree dynamic system subjected to simplified loads 
there are five independent variables, B, T, m, k, and . When any 
one of these variables is changed, the maximum displacement, in general, is 
affected. The total effect on x^^^ when all the variables are changed can 
be obtained approximately by the method of superposition. Thus the total 
increment Ax , for the general case instead of being given by equation 
(A.I), is given by: 


Ax = 
m 


AT + 


Ak + 


AB + 


(A.4) 


According to equation (A.2), the different partial derivatives in equation 
(A.4) can be defined as different sensitivities. 

Applying the technique used in deriving equation (A. 3 ), the general 
expression for the percentage-increment in the maximum displacement can be 
derived and this is given by eqmtion (A. 5 ). 


xT T 


+ T ^ + T — + T ^4-T _! 

xm m ^xk k -^xB B -^xR R 


(A. 5 ) 


same method of derivation, if W instead of x is consider^ as 

in m j 

the independent variable, the general expression of the percentage-injcrement 
of the niaxinium work done is given by equation (A.6): 




(a\6 ) 


The quantities, I^, I^, I^, I^, and in equation (A. 5 ) are the 

percentage-increment ratios for the maximum displacement, while the quanti-\ 
^ W W ^ ^ ^ equation (A.6) are the 

percentage-increment •ratios for the maximum work done. 

The percentage-increment ratio gives quantitatively the effect of a 
given veuriable on the response of the dynamic system. For example if the 
numerical values of for different structures are different, the 

larger the value of is, the more the maximum displacement will be af¬ 

fected by any change in the load duration. Similarly, for a given structure 
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if Ip is nmerically larger than I , the maximum displacement will he 

Xii xm 

affected more by the plastic resistance than by the mass, m . On the other 
handj if is numerically small, say less than 0 . 2 , the maximum dis¬ 

placement is practically a constant when the spriiog constant k varies. 

Equations (A. 5 ) and (A.6) are derived by partial differentiation ant^ 
superposition. Hence, in the evaluation of the percentage change in x^ 
or W^. for a given structTore dynamic system, these simplified equations, 
in general, are applicable only when the increments in B, T, m, k, or 

R are not more than 10 ^. 
in 

In equations (A. 5 ) and (A. 6 ), there are ten percentage-increment 
ratios. These ratios are not all independent. The relationships among 
them are derived in paragraph A- 03 . 

A -03 REIATIOHSHIPS AMONG TEE mCEMAOE-IHCTlEMElilT RATIOS, a. Increment 
Ratio for Non-Dimensional Parameters. The discussion in paragraph 5-10 
shows that only three non-dimensional parameters are needed to express the 
functional relationships of all the variables of an elasto-plastic system 

subjected to simplified loads. For example in figures 5*25 and 5*28, the 

X m ^ 

in T iTL 

three parameters are, “ ^ f which are denoted by C^, 0 ^, 

^e "^n 

and C„ respectively. That is 
K 


■'R 


_m 

”” X 

(A.7) 

e 


ii 

(A.8) 

R 


m 

~ B 

(A.9) 


In terms of these parameters, the dependent variable of the dynamic system 

is C instead of x , and the independent variables are 0^ and C„ 

X m ' ^ T R 

instead of B, T, m, k and R . For small increments in C-, and 0^, , 
' ' m T R ' 

the corresponding increment is given by equation (A.IO). 

8c ac 
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ac be 

where and are small increments, and are par¬ 
tial derivatives. ^ ^ 


Applying the same technique used in the derivation of equation (A. 3 ) 
and (A. 5 ), the percentage-increment of , can be derived and is given-by 
equation (A,11) 


'fl 


\ 

.''x^ 


T (A-11) 


where 



= the percentage-increment ratio of with respect to • 

Tliere are only two percentage-increment ratios in equation (A.ll). 
It can be concluded that of all the percentage-increment ratios of an 
elasto-plastic system, only two ratios are independent. Hence the ten 
percentage-increment ratios in equations (A.5) and (A.6) can be expressed 
in terms of and • 

This conclusion can also be deducted by physical reasoning. is 

a function of Cj^ and . If for a given case is a constant, 

would remain a constant provided that is unchanged although T , m 

and/or k may be varying. This indicates that I^# 

be expressed in terms of • Similarly, it can be concluded that ^ 
and I „ can be expressed in terms of ■ The equations expressing 
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those relationships are given in paragraph A-03'b. 

b. Percentage-Increment Ratios in Terms of and Substi¬ 


tuting for X in equation (A. 7 )^ the following equation is obtained: 

K, 0 



X k 

m 

R 

m 


(A. 12) 


Taking logarithms of both sides of this equation, then differentiating, and 
then replacing the differential notation "d" by the increment notation 
"a," equation (A.I 3 ) is obtained. 






R 


m 


(A. 13 ) 


Similarly, from equations (A.8) and (A. 9 ), the following two equations are 
obtained respectively: 


- T ^ 2 k 
^m 



substituting equations (A.I 3 ), (A.l4) and 
rearranging: 


2 m 


(A.14) 


B 


(A. 15 ) 


(A. 15 ) into equation (A.ll) and 



^ i ^ (1 i I ^ 

IT ■ 2 m 2/k 



(A. 16) 


Comparing equation (A.I 6 ) with equation (A. 5 ), the following relationships 
among the percentage-increment ratios are obtained: 


1 ■ xT 


^2 " " ^xB “ (^xR 




(A. 17 ) 
(A.18) 


So far in the derivation, x^ is considered as the dependent variable- By 
the same method, the percentage-increment ratio when is the dependent 

111:- 



A-03b 
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variable can be derived. For an elasto-plastic system, referring to equa¬ 


tion (5- 38 )# the expression for is given by 


m in\ 


(A.19 ) 


Taking logarithms of both sides, differentiating and rearranging: 


2 C Ax 2 /C - 1 ) AR 


V ~ 2C - 1 X 
m X m 


1 R 2C - 1 k . 
m X 


(A. 20 ) 


Substituting equation (A.I 6 ) into (A.20), the expression of in terms 

vni 

of the percentage changes of the dependent variables is obtained; 


W 2C - 1 -^1 T 
m X 


2 / C 

jc T ^ 1-1 X 

- 1 •*•1 m " " 2C^ - 

I ^ + (r 

2C^ - 1 -^2 B r 


2C^ - 1 ^ij k 


: \ AR 

-X_ I \ 
-1-^21 R 

/ ® 


(A.21) 


Conqjaring equation (A.21) with equation (A.6), the relationships among the 
percentage-increment ratios of are obtained and these are given by 
equations (A.22) and (A. 23 ). 



(A.22) 
(A.23) 


Equations (A.20) through (A.23) are derived on the basis that is 
given by equation (A.I 9 ), which is true for an elasto-plastic system. It 
C€Ui be shown that these equations are also applicable for both linearly 
elastic and completely plastic systems provided that is replaced by 

unity for the former while the limiting value with approaching in¬ 

finity is used for the latter. 

Equations (A.l?) and (A.22) show that t t y t 

and can be expressed in terms of t while equations (A.I 8 ) ai^ (A. 23 ) 
show that and can be expressed in terms of 1^ . 

Therefore, the graphical presentation of the percentage-increment ratios 
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is greatly sin^ilified by these functional relationships. When any one of 
the percentage-increment ratios in the first group is computed, all the re¬ 
maining I*s are also known. The same is time for the second group of 
I's . Thus only two sets of cxarves are needed for the con 5 )lete graphical 
presentation of all the I's of an elasto-plastic system subjected to 
sinqplified loads. In this appendix, one set of curves is for , which 
is also equal to I^ip . The other set is for I^ which is equal to 
^1- + Ig^ . These curves are described in paragraph A- 06 . 

c. Other Percentage-Increment Ratios. The discussion so far covers 
the cases when the independent variables are B, T, m, k, and and 

the dependent variables ahe and . For special application, the 

percentage-increment ratio for variables other than those mentioned above 
may be needed. Because of the non-dimensional natvire of the percentage- 
increment ratio, the expressions for additional percentage-increment ratios 
can be derived from those already presented by a method of superposition. 

A few examples are given for illustration. 

As a first example, suppose the effect of the limiting elastic dis¬ 
placement x^ on Xjj^ is to be considered. Since 

R 

m 


therefore 


Ax ZiR 
e m 



k 


(A.2l^) 


The percentage-increment ratio of x^^ with respect to 
by; 


I = I t, - I 1. 
xe xR xk 


X 


e 


is then given 
(A. 25 ) 


Where is ratio of the percentage change in x^^^ to the corresponding 

percentage change in x^ . 

As a second exanple, the load jjulse H , which is the area under the 
load time curve as defined by equation (5'5)> is considered as an inde¬ 
pendent variable. For rectangular or triangular load, the load pulse is 
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proportional to the product of B and T . Then equation (A. 26) is 


obtained: 


^xH " ^xT 


{A.26) 


where ^xH is the percentage-increment ratio of the maximum displacement 


Xj^ with respect to the load pulse H 


As a third example, the effect of ^ on the maximum displacement is 
considered. For a given trieuigular or rectangular load, this -variable is 
the ratio of the two sides of the triangle or rectangle. Hence in a very 
restricted sense, it represents the "shape" or "proportion" of the load 
curve with a given geometry. Let this variable be denoted by J , that is: 


The percentage-increment ratio is given by: 


(A. 27) 


^xJ "" ^xB " ^xT 


(A.28) 


Where I , is the ratio of the percentage change in x to the correspond- 
xj _ m 

ing percentage change in ^ . 

A-04 FERCTamGE-INCTai>^ RATIOS FOR LIHEARI2C ELASTIC SYSTEMS. The 

percentage-increment ratios for the maximum displacements of linearly 

elastic systems subjected to triangular or rectangular loads are given in 

this paragraph. In deriving the expressions for different percen-tage- 

increment ratios, the methods described in paragraphs 5-07b, A-02 and A-03 

T 

are used. For exaiiq)le, the case of a triangular load with > 0*371 is 
considered. This corresponds to t^^^ < T . The expression for for 

this case is given by equation (A.29)* 


X = r 1 
m k I 


-1/2 JtT' 


(A.29) 


The expressions for I derived from this equation are re- 

XJd XI 

spectively given by the following two equations: 


"xB = 


(A.30) 
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^xT 





2jt 

T 


1 + 


1 

T J 



^ _ 1 


2rt - 


tan 2«Cn 


T 


(A.31) 


Applying equation (A.I 7 ), the corresponding expression for and 

can he derived from equation (A. 3 l)* 

Similar expressions for all the percentage-increment ratios for tri- 
T 

angular loads with ^ < O. 37 I , and also for rectangular loads can be ob¬ 
tained. These e^qaressions are all summarized in paragraph A-03a. 

The percentage-increment ratios for lineeirly elastic systems are 
plotted in figures A.l and A.2. The first figure is for rectangular load 
while the second is for triangular load. From these two figures;, the fol¬ 
lowing conclusions can be drawn: 

( 1 ) For both types of load, is a constant and equal to unity. 

The maximum displacement of a linearly elastic system is always directly 
proportional to the peak load B . 

( 2 ) When is less than 0.2, is approximately constant and 

equals unity for both types of loads. Thus for short duration load, the 
maximum displacement of a linearly elastic system is directly proportional 
to the load dvtration T . On the other hand, for a rectangular load, I „ 
is zero when > 0.5 an<i for a triangular load, is less than 0.1 

when > 3'0- Therefore the maximum displacement is practically inde¬ 
pendent of the loeui diiration for long dioration load. 

( 3 ) For both types of loads, when C_ < 0.2, both I and I are 

approximately equal to -0.5- Thus for rectangular or triangular loads of 
short duration the maximum displacement is inversely proportional to both 
the square root of the mass m and the square root of the spring constant 
k . On the other hand, for a rectangular load when Ct> 0. 5> equali 

zero and equals -1.0. These numerical values are also approximately 

correct for a triangular load when > 3 . 0 . Thus for long duration tri¬ 
angular’or rectangular load, the maximum displacement of a linear elastic 
system is- practically independent of the mass hut is inversely proportiona] 
to the spring constant k . 


no 


0.02 


0.2 


0.3 0.4 


0.6 0.8 l.( 


0.03 0.04 0.06 0.08 0.1 


Figure A A. Percentage increment ratios for a linearly elastic system 
subjected to rectangular loads 



Figure A.2. Percentage increment ratios for a linearly elastic system 
subjected to triangular loads 
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(4) The numerical values of 1^ and are always 

less than unity, llius the reliability of the maximum displacement de¬ 
termined by a design computation is better than or as good as the uncer¬ 
tainties in the determination of any variable. 


A -05 PERCEHTAGE-CTCREMEMT RATIOS FOR COMPLETELy PLASTIC SYSTEMS. Applying 
the method described in paragraph 5 - 07 c, the expression of x^^ for a com¬ 
pletely plastic system subjected to rectangular lo8d,s can be obtained and 
is given by equation (A.32 ). 



(A. 32 ) 


The percentage-increment ratios and derived from equation 

(A. 32 ) are given by: 


= 2 {A. 33 ) 



The corresponding expressions for ^ 

can be obtained from equations (A.IT)# (A.I 8 ), (A. 22 ) and (A,23). Similarly 
the corresponding expressions for triangular loads can be derived. These 
cjqpressions are summarized in paragraph A-08b. 

Th e percentage-increment ratios for both and of a completely 

plastic system subjected to rectangular and triangular loads are given in 
figures A. 3 and A.4 respectively. From these two figures, the following 
conclusions can be drawn: 

( 1 ) For both types of loads, and are eqxxal to 2, and 

and are equal to \niity. Hence, for a completely plastic system sub¬ 
jected to simplified loads the maximum displacement and the maximum 

work done are: (l) proportional to the square of load duration T 

and ( 2 ) inversely proportional to the mass m . 

( 2 ) The numerical values of or are always greater than 2 

and the nvimerical value of is always greater than unity. When 

< 0 . 3 , and are approximately equal to 2 and -1 respectively. 

:^c 
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^Fi^are A^3. Percentage increment ratios fora 
completely plastic system subjected to 
rectangular loads 


Figure A.4, Percentage increment ratios for a 
completely plastic system subjected to 
triangular loads 


Hence, for such cases the inaximuin displacement is proportional to the 
sq.uare of the peak load, B , and inversely proportional to the plastic re¬ 
sistance, R . 

(3) When Cp > 0 . 5 , for rectangular loads, > 3 and < -2; 

and for triangular loads, > 4 and < - 3 . Thus the maximum dis¬ 

placement is extremely sensitive to variations in B and R . For such 

m 

cases it would he meaningless, in a trial and error design procedure, to 
seek an agreement better than + 30 ^ between the maximum displacement obtained 
by a design cbmputation and that of the design specification. 

A -06 FERCENTA-GE-IHCREMEMT RATIOS FOR EIASTO-PIASTIC SYSTEMS. The expres¬ 
sions for the percentage-increment ratio of an elasto-plastic system sub¬ 
jected to simplified loads can be derived by the methods described in para¬ 
graphs 5-07d» A -02 and A-03- An example is used to illustrate the procedure 
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in the derivation. The independent variables of the dynamic system chosen 

as example satisfy the following conditions: (l) the load is a triangular 

load] ( 2 ) the displacement reaches the limiting elastic displacement 

at time t where >.t < T; and ( 3 ) the maximum displacement is reached at 
0 0 

t where t < T . 
in m 

When t < t , the equation of motion is given by equation (5-33a). 

6 

When t < t < t .the equation of motion is given by equation (5-33b). 
e m 

The expression for x^ derived from these two differential eq\jations is 
given by equation (A.35)* 



(A.35) 


where 


V = the velocity of the mass at time t . 

• e ® 

This equation shows that x^ is a fxinction of B, T, m, k, t^, 

t anH V . The last three variables: t^, t^, and v^, in turn, de¬ 

pend on B, T, m, k, and • Considering the effect of load duration 
T on x^ , the general expression of , instead of being given by 

equation (A.3), is given by equation (A. 36 ). 


xT 


T 


m 


-sf*3r-5f"5r-sf*sr'^ 

e m 0 


(A. 36 ) 


After the partial derivatives in equation (B. 36 ) are evaluated, 
Sion of is given by: 



the expres- 


(A. 37 ) 


By a similar procedure, the expressions for all the percentage-increment 
ratios of an elasto-plastic system subjected to triangular or rectangular 
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loads can be derived. These expressions are summarized in para^aph A-08c. 

The percentage-increment ratios for elasto-plastio-systems subjected 
to rectangular aind triangular loads are given in figures A.5 through A.8. 

Ihe first two figures are for rectangular loads and the last two, for tri¬ 
angular loads. In figures A. 5 and A.7, the percentage-increment ratios 
with respect to T, m and k are plotted against with as the 

running parameter. In figures A.6 and A.8, the independent variables are 
B and . IVom these fovir figures, the following conclusions can be 
drawn; 

(1) The effect of load dijration T . The curves in figires A.5 and 

A. 7 are divided into two groups by the ciirve corresponding to = 1 . 

When < 1, the value of I^^ , with a few exceptions, lies between 1 
and 2. Hence for a given C_ , the maximum displacement is proportional to 
T at small values of and is proportional to the square of T at large 

values of . The smaller the values of , the smaller will be the 

value of when I^ approaches 2. On the other hand when > 1.0, 

the value of I^^ lies approximately between zero and unity. Hence, at 
small values of , the mpyimum displacement is approximately propor¬ 
tional to the load duration T , while at large values of C^, , the maxi¬ 
mum displacement is practically independent of the load duration T . 

( 2 ) The effect of the mass m and the spring constant k on the 

maximum displacement can be approximately summarized into the following two 
statements, (l) When Cp < 1, the values of lie approximately be¬ 
tween -O .5 and -1.0, and the values of li® between -0.5 and 0. Hence, 

at small values of , the maximum displacement is inversely proportional 
to the square root of both the mass, m , and the spring constant, k 

At large values of C„ , x is inversely proportional to m but is inde- 

pendent of k . (2) When 1.0. At small values of \ ® 

versely proportional to square root of both m and k . But at large 

values of C„ , x is inversely proportional to k but is independe 
T m 

( 3 ) The effect of the peak load B and the plastic resistance 

As shown in figures A.6 euid A.8, the numerical values of xR 

scatter over a much wider range than the corresponding value 


123 




igure A.7. Percentage-increment ratios and for Figure A, 8. Percentage^increment ratios 7^^ ^xB 

on elasto^plastic system subjected to triangular loads elasto-plastic system subjected to triangular IocmIs 
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convenience of discussion, the curves in these two 
figures can be very roughly divided into three regions, (l) The first 
region covers all cases when < 0.6. Practically for all values of Cp , 
the values of lie between 1 and 2, and the values of lie be¬ 
tween 0 and -1. Thus for low values of is directly proportional 

to B but independent of while for high values of 0^ , x^^ is pro¬ 

portional to the square of B and inversely proportional to R . 

m 

parameters of one story buildings are generally located in this region. 

Since the ninuerlcal values of both I _ and I _ ere not very large, a 
close agreement between the maximum displacement obtained by a design com¬ 
putation and that of the design specification cem generally be achieved. 

(2) The second region corresponds to cases when > 1-0 ^ ^ 1.4 or 

Cp < 0.4. In this region both and are practically constant.. 

The numerical values of neither one is very large. The cases when Cj^ < 0.4 
are not of practical importance because the coiresponding values of are 
always greater than 20. On the other hand, the cases when > 1.4 gen¬ 
erally correspond td the region in which the design of structtffal members 
such as beams, walls or columns are located. ( 3 ) The third region cor¬ 
responds to cases when > 1.0 and 0.4 < < 1.4. In this region I^ 

is generally greater than 2.5 and I^ is less than -I. 5 . Therefore, the 
maximum displacement x^ is very sensitive to any variation in B or R^^^ 

especially when C_ lies between 0.8 and 1.2. Ih design problems, both B 

R 

and R are subject to tincertainties in their determinations. Therefore 

the value of x obtained by a design conqputation would not be very re- 
in 

liable. Hence it would be a good design practice, in general, to select 
the parameter of a design such that the combination of and Cp Ilea 

outside of the third region. 

A -07 THE EFFECT OF LOAD SHAPE ON THE RESFOHSE OF DIHAMIC SYSTEMS. The 
effect of load shape is partially investigated in EM 1110-345-415 where the 
responses of a dynamic system subjected to different types of load, as shown 
in figure 5 . 8 , are analyzed and presented separately. For each type of load 
with a given geometry, l.e., a triangle or a rectangle, its "shape" or "pro¬ 
portion" is completely specified by two variables. These are the peak load 
B and the load duration T . When either B or T is varying, the 
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response is also varying continuously. Therefore, in paragraph A-02, par¬ 
tial derivatives of x with respect to B or T can be evaluated and 

m 

the effect of B and T can be expressed q.uantitatively in terms of 
percentage-increment ratios. Here again the percentage-increment ratios for 
rectangular and triangular loads are considered and presented separately. 

If the concept of percentage-increment ratio is to be extended to the 
analysis of the effect of load shape on the response, it would not be siaffi- 
cient to consider Just discrete types of loeid, for example rectangular and 
triangular loeuis. Instead, another variable has to be added to the load 
pulse so as to form a family of load curves, and the geometrical shape-of 
the load curve will then be varying in a continuous fashion. For exan^ile, 
in figure A. 9 , when T^ varies from 0 to T , the geometrical load shape 

, . varies continuously from a triangular 

f(t) f(t) »(t) f(t) ^ 

B B— V B -- B - 1 load to a rectangular load. A 

N. 1 \ 1 \ percentage-increment ratio I » 

I——t i y -t I—-4—* 

^ . ' could then be defined which expresses 

Trtongulor T| Increoting Raetongulor ^ ... ^ 

LoodT, *0 LooaT,'T the relative magnitude of the per- 

j -.IT centage changes in x„ and T, . 

Figure A.9. A set of load curves with 1 j mi 

as a parameter And the effect of changing the geo¬ 

metrical shape of the load on the 
response could then be analyzed quantitatively. 

It is clear that a rigorous investigation of the effect of the load 
shape cannot be carried out without introducing another variable. For an 
elasto-plastic system, there would be 4 non-dimensional variables Involved 
and the result could not be conveniently presented in graphical form. Al¬ 
though for linearly elastic or completely plastic systems there would be no 
difficulty in presenting I _ for a given family of load shapes. Its 

xii 

practical application would be limited to the case where the actual load 
shape variation follows the same pattern as the variation in the given 
family. Because of these difficulties, as far as the effect of load shape 
is concerned, a qualitative discussion, supplemented by non-dimensional 
grains whenever convenient, is presented in this paragraph. Specifically, 
two cases are considered. In the first case, the response of an elasto- 
plastic system to rectangular and triangular loads is compeared with each 


Figure A.9. A set of load curves with Tj 
05 a parameter 
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other. This will indicate in what region^ if any, the load shape is not 
important. In the second case, the percentage-increment ratio I t ^ given 
hy equation (A. 28 ), for rectangular and triangular loads is presented in 
graphical form. These curves will show how the response is affected by 
changing the "shape" or "proportion" of the load curve without changing its 
geometry. 

Comparison of the Responses of an Elasto-Plastic System Subjected 
to Rectangular and Triecngular Loads. In figure A. 10, the ratios of the 
maximum displacement for an elasto-plastic system subjected to a rectangular 
load to the corresponding value for a triangular load are plotted. There 
are three regions in this figure where the shape of the load pulse has no 



Figure A.10. Ratio ofnrnximum displacements for rectangular to the corresponding values 
in triangular load, elasto-plastic system 
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Importairt effect on the response, (l) In iTirst region, where > 0.5 

and C_ < 0.4, the ratio of x^j^’s is eq.\xal to 4. The cor- 

Xt .1 

responding values of and are a greater than unity. In this 

region, the maximum strain energy is appi ox jLBnQ_.j.g 2 y proportional to x but 

. HI 

the maximum work done is approximately propo^r-t ional to the square of the 
area under the load-tin^ curve* WiereforCj^ *fctte maximum displacement is ap¬ 
proximately proportional to the square of area under the load curve 

while the detail geometry of the load curv^ in^K^rtsoit. (2) In the 

second region where < 0^5 ’t.he ratio of x^’s is approxi¬ 

mately equal to 2. In this regiotii x ■< SDc and t > T . Althou^ W 
is still approximately proportional to' the square of the area of the load 
pulse, the nBximum strain energy is also pr'ojjortional to the square of . 
There the maximum displacement is proport iorxeul to the area under the load— 
time curve and is not much affected by the <3Le-tail shape of the load. ( 3 ) 

In the third region vhere ^ 2 at^ >- 14.2, the ratio of xi^'s is ap¬ 

proximately equal to vinity. This corrcspoxacLs to the cases where t^ is 
mTiA.n In comparison with T . The faeactBiuiEi displacement is proportional to 
the peak load B while the load shape 1« rko-fc important. 

Outside of the above three regions, t;lxe! ratios of x^'s for rectangu¬ 
lar and triangular loads vary considerably wi-th both and Cj^ . No 

sinrple conclusion can be drawn. Generally ■fctie response is expected to be 
affected by the detail geometry or the loafed curve. 

b. The Percentage-Increment Ratio • The percentage-increment 

ratio I - given by equation {A. 28) it? plo'fc'ted in figures A .11 and A.12 
for rectangular and triangular loads respesotii-vely. As discussed in para¬ 
graph A- 03 c, the value of I^ Indicates effect of varying the "shape" 

or "proportion" of a load curve without rtueuciSing its basic geometry. In 
figures A. 11 and A. 12, the numerical valu«s of' ax-e small when 

< 0 . 3 . In this region the response wt XX not be much affected when one 
triangular load is replaced by another tr* ia-xxgular load, or one rectangle 
by another rectangle, or any load curve %rln.Xo;h resembles a triangle by th.e 
corresponding triangular load. 

A -08 EXPRESSIONS FOR MAXIMUM PERCTIIWAGE-INCKEMEI^ RATIOS 

The expressions of t , x , I anit linearly elastic. 


1,2c 


Trtongulor Rttlstonct 
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completely plastic and elasto-plastic systems subjected to rectangular and 
triangular loads are given in this paragraph. From the expressions of 
these fo-ur Quantities, the expressions of many other quantities can be de¬ 
rived. For example, in terms of x , the expression for W is: 

' m ■' m 

when X < x 
m e 


and when x > x 
m e 


W = ^ k x^ 
in 2 m 



(A.38 


(A.39 


In terms of I ™ and I „ , the expressions of many other percentage- 
XT XK 

increment ratios can be obtained. These relationships are given by equa¬ 
tions (A. ho) and (A.hi). 



a. Expressions for Linearly Elastic Systems. 

Rectangular Loads 

(i) For t < T 
m 

t ^ 
m _ 1. 

T - 2 

k X 

D.L.F. =V^ = 2 

where 

D.L.F. = dynamic load factor 


(A.ho) 


(A.hi) 


(A.h2a) 

(A.h2b) 


xT 


= 0 


(A.h2c) 


130 


A-08a 


(ii) For t > T 
m 




cos 2n - 1 


\ 1 -if 2« C, 

[isrs, Cj 

D.L.P. = 2 sin C^j 
^xT “ 


D.L.P. = 2 sin 
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(A.42d; 


(A.43a) 


(A. 43b) 


(A.43c) 




(A.43d) 


Triangular loads 


(i) For t < T 


m 1 . -1 


n p 
).L.F. =21 


(a- s) 


tan 2it C„ 


(A.44a) 


(A.44b) 


r 1 , 2 „ 2 1 r 

I . .1 . ■ °T _ 

xT ^ P -"I- T " --. ■ “1 

1 + 4rt Cj 2« - tan”^ ^2n 


(ii) For > T 


D.L.F- - 




■■ T ” 2it 
n 


^xR=°' 


1 - cos 2jt 
sin /2it C^\ - 25r 


2 + 4«^Cj^ - 2 COB 


^2jt - 4« sin ^2« C^j 2 


I m —X ■f' 
xT 


1 - COS ^2k 

D.L.F. 


(A.44c) 


(A.44d) 


(A.45a) 


(A. 46b) 


(A.46c) 


I _ = 0 
xE 


(A.46d) 
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Expressions for Cong)letely Plastic Systems. 
Rectangular Load 


" 1 

m 1 


Triangular Load 


"R 


"R 


T 


X 

m 


2 m 

■ S 

^xT 

» 2 

T « • 

1 

^xR 


Sind t < T 
m 


t 

-£ = 2 1 
T ^ 1 

(-S) 

*m 4 

(l C 

12 0 ^ *3 

2 m • 

\ ?, 

^xT 

- 2 

T _ _ 


: 

^ “ ^R 

and t > T 
m 


t 

m 

1 

■■"au"' IBS 

T 

2 Cr 
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c* Expressions for Elaato-Plastic Systems. 

Rectangular Load 

( 1 ) For and < T . .The expressions for this 

the same as those given by equation (A.I 2 ). 

(il) For x^ < Xg and t^^ > T . The expressions for this 

the same as those given by equation (A.43). 

(ill) For x„ > X , t > T and t > T . 
m e' m e 
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(iv) For X > X , t < T and t < T . 

HI e m 


(v) For 
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Itiajigiilar Load 

(i) For X < X , and t < T . 'Hie expressions are the same a 

' m e ^ m 

those given by equation (A.1+1+)- 

(ii) For X < X , and t > T . The expi’essions are the same 
those given by equation (A. 1+5). 

(iii) For x > x , t > T, and t > T . 
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(iv) For X > X , t < T, and t < T 
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(v) For x^>x^, t^< T and > T 
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